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20.  Abstract  (Continued) 

u  %  j*  Tr®cluencief  °f  significance  for  long  range  surveillance  systems  (5  to  100 
Hz)  data  for  model  building  is  often  insufficient  or  too  poorly  resolved.  In  spite 
of  the  vast  quantity  of  marine  geological  and  geophysical  data  generally  available 
to  the  acoustician,  a  majority  of  the  data  describe  sediments  and  crustal  struc¬ 
ture  at  depths  greater  than  a  few  hundred  meters  or  less  than  10  meters  Cor¬ 
respondingly,  these  depths  are  important  at  frequencies  below  a  few  Hertz  or 
greater  than  a  few  hundred  Hertz. 

The  purpose  of  the  study  documented  in  ttiis  report  was  to  develop  geoacous- 
t^c  models  capable  of  predicting  the  essential  bottom  reflectivity  behavior  relevant 
to  long-range,  low-frequency  surveillance  systems  from  acoustic  survey  measure- 
ments  and  available  marine  geological  data. 

It  was  found  that  at  five  thickly  sedimented  bottom  sites,  simple  geoacoustic 
models  could  be  derived  from  measurements  of  bottom  reflectivity  data  taken  at 
grazing  angles  from  5°  to  80°.  Deconvolution  processing  and  analysis  of  the 
resultant  bottom  impulse  responses  led  to  the  form  of  the  models  and  to  quanti¬ 
tative  values  of  model  parameters.  Specifically,  amplitude  and  arrival  time  analy¬ 
sis  ot  dominant  peaks  in  the  impulse  responses  was  employed.  Model-predicted 
?.rknva) l.Jf,mes  differed  from  measured  arrival  times  by  3.1  msec  on  the  average. 
This  difference  is  small  compared  to  the  processing  pulse  duration  of  14.2  msec 
Accurate  arrival  time  prediction  such  as  exhibited  here  is  one  requirement  of  a 
geoacoustic  model  if  it  is  to  be  incorporated  into  a  propagation  model  for  detailed 
predictions  of  acoustic  field  intensity. 

Geological  analyses  of  available  environmental  data  were  conducted  to  develop 
estimates  of  sediment  sound  speed,  density  and  attenuation  profiles.  Geoloqical 
data  were  obtained  from  (1)  cores,  (2)  seismic  profiler  records,  (3)  3.5  kHz  pro¬ 
filer  records,  (4)  wide  angle  reflection  profiling  (sonobuoy  data),  and  (5)  Deep 
Sea  Drilling  Project  (DSD)  holes  drilled  in  locations  with  sediment  analogous  to 
the  site  of  interest.  Information  relevant  to  acoustic  properties  are  summarized 
including:  water  depth,  bottom  and  water  sound  speed,  sediment  thickness 
and  reflectors  within  the  sediments.  General  discussions  of  geological  features 
dixD m9  aCOUSt'C  ProPerties  at  several  of  the  stations  are  presented  in  Appen- 

Geologically  inferred  predictions  of  sediment  surface  sound  speed  and  density 
were  in  agreement  with  the  acoustical  estimates  but  the  total  sediment  sound- 
speed  model  was  inadequate,  as  judged  by  errors  in  the  predicted  arrival  times. 
Bom  the  acoustical  and  geological  estimates  of  attenuation  were  subject  to  large 
uncertainties,  and  the  geological  estimate  was  consistenly  higher  than  the  acous¬ 
tical  estimate. 

The  acoustical  model  results  of  this  study  suggest  generalization  to  other 
deep  ocean  abyssal-plain  and  abyssal-hills  bottom  sites.  However,  since  the  geo¬ 
logical  data  base  was  of  the  type  that  is  generally  available,  the  difficulties  with 
geological  model  prediction  may  also  be  encountered  at  other  sites. 

Recommendations  for  improving  accuracy  of  geologically  inferred  models 
include  (1)  recording  and  analyzing  full  dynamic  range  of  amplitude  information 
in  normal-incidence  profiler  data;  (2)  pulse  compression  and  dereverberation 
processing  (deconvolution)  of  profiler  and  wide-angle  reflectivity  data;  and 
(3)  improved  sampling  in  DSDP  cores  of  the  top  few  hundred  meters. 
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Foreword 


This  report  documents  geoacoustic  models  designed  to 
predict  the  nature  of  acoustic  boundary  interaction. 
For  five  sites  involving  thickly  sedimented  bottoms, 
simple  models  are  derived  from  measurements  of  reflec¬ 
tivity.  Also,  recommendations  are  provided  for  improv¬ 
ing  the  accuracy  of  geoacoustic  models  in  general. 

>0.  O.  P-faOt**. 

G.T.  Phelps,  Captain,  DSN 

Commanding  Officer,  NORDA 
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Executive  Summary 


The  performance  of  surveillance  systems  mounted,  sus¬ 
pended,  or  towed  near  the  bottom  is  governed  in  part 
by  the  nature  of  the  acoustic  bottom  interaction.  Bot¬ 
tom  properties  affecting  this  interaction  include, 
among  others,  the  spatial  variation  of  sound  speed, 
density,  and  attenuation.  Specification  of  these  bot¬ 
tom  properties  is  generally  referred  to  as  a  geoacous¬ 
tic  model.  Research  to  date  has,  to  a  large  degree, 
involved  the  quest  for  appropriate  geoacoustic  models 
and  their  validation  by  comparison  of  measured  and 
model-predicted  bottom-loss  functions . 

At  frequencies  of  significance  for  long  range  surveil¬ 
lance  systems  (5  to  100  Hz)  data  for  model  building  is 
often  insufficient  or  too  poorly  resolved.  In  spite  of 
the  vast  quantity  of  marine  geological  and  geophysical 
data  generally  available  to  the  acoustician,  a  major¬ 
ity  of  the  data  describe  sediments  and  crustal  struc¬ 
ture  at  depths  greater  than  a  few  hundred  meters  or 
less  than  10  meters.  Correspondingly,  these  depths  are 
important  at  frequencies  below  a  few  Hertz  or  greater 
than  a  few  hundred  Hertz. 

The  purpose  of  the  study  documented  in  this  report  was 
to  develop  geoacoustic  models  capable  of  predicting 
the  essential  bottom  reflectivity  behavior  relevant  to 
long-range,  low-frequency  surveillance  systems  from 
acoustic  survey  measurements  and  available  marine  geo¬ 
logical  data. 

It  was  found  that  at  five  thickly  sedimented  bottom 
sites,  simple  geoacoustic  models  could  oe  derived  from 
measurements  of  bottom  reflectivity  data  taken  at 
grazing  angles  from  5°  to  80°.  Deconvolution  pro¬ 
cessing  and  analysis  of  the  resultant  bottom  impulse 
responses  lead  to  the  form  of  the  models  and  to  quan¬ 
titative  values  of  model  parameters.  Specifically, 
amplitude  and  arrival  time  analysis  of  dominant  peaks 
in  the  impulse  responses  was  employed.  Model-predicted 
arrival  times  differed  from  measured  arrival  times  by 
3.1  msec  on  the  average.  This  difference  is  small  com¬ 
pared  to  the  processing  pulse  duration  of  14.2  msec. 
Accurate  arrival  time  prediction  such  as  exhibited 
here  is  one  requirement  of  a  geoacoustic  model  if  it 
is  to  be  incorporated  into  a  propagation  model  for 
detailed  predictions  of  acoustic  field  intensity. 


Executive  Summary 


Geological  analyses  of  available  environmental  data 
were  conducted  to  develop  estimates  of  sediment  sound 
speed,  density  and  attenuation  profiles.  Geological 
data  were  obtained  from  (1)  cores,  (2)  seismic  pro¬ 
filer  records,  (3)  3.5  kHz  profiler  records,  (4)  wide 
angle  reflection  profiling  (sonobuoy  data)  and  (5) 
Deep  Sea  Drilling  Project  (DSDP)  holes  drilled  in 
locations  with  sediment  analogous  to  the  site  of  in¬ 
terest.  Information  relevant  to  acoustic  properties 
are  summarized  including:  water  depth,  bottom  and 
water  sound  speed,  sediment  thickness  and  reflectors 
within  the  sediments.  General  discussions  of  geologi¬ 
cal  features  affecting  acoustic  properties  at  several 
of  the  stations  are  presented  in  Appendix  D. 

Geologically  inferred  predictions  of  sediment  surface 
sound  speed  and  density  were  in  agreement  with  the 
acoustical  estimates  but  the  total  sediment  sound- 
speed  model  was  inadequate,  as  judged  by  errors  in  the 
predicted  arrival  times.  Both  the  acoustical  and  geo¬ 
logical  estimates  of  attenuation  were  subject  to  large 
uncertainties,  and  the  geological  estimate  was  consis¬ 
tently  higher  than  the  acoustical  estimate. 

The  acoustical  model  results  of  this  study  suggest 
generalization  to  other  deep  ocean  abyssal-plain  and 
abyssal-hills  bottom  sites.  However,  since  the  geolog¬ 
ical  data  base  was  of  the  type  that  is  generally 
available,  the  difficulties  with  geological  model  pre¬ 
diction  may  also  be  encountered  at  other  sites. 

Recommendations  for  improving  accuracy  of  geologically 
inferred  models  include:  (1)  recording  and  analyzing 
full  dynamic  range  of  amplitude  information  in  normal- 
incidence  profiler  data;  (2)  pulse  compression  and 
dereverberation  processing  (deconvolution)  of  profiler 
and  wide-angle  reflectivity  data;  and  (3)  improved 
sampling  in  DSDP  cores  of  the  top  few  hundred  meters. 
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in  Figure  53.  Heavy  dashed  lines  define 
sediment  thickness  of  1.0,  1.5,  and  2.0 
seconds  two-way  travel  time  observed  on 
profiler  recored  (modified  from  Ref. 

88).  DSDP  Site  118  in  the  Biscay  Abys¬ 
sal  Plain  is  also  indicated.  Cores  in¬ 
dicated  were  taken  by  R/V  VEMA  (V)  and 
USNS  KANE  (K). 

Figure  53.  Seismic  profile  records  across  the  76 

Northern  Iberian  Abyssal  Plain.  Track 
is  indicated  in  Figure  52.  Profile  was 
obtained  on  USNS  KANE.  Sections  A  and 
B  and  Units  1  and  2  are  discussed  in 
text . 
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at  Station  2.  Heavy  line  is  based  on 
gradient  derived  by  the  CRIC  method  at 
DSDP  Site  118.  The  Houtz-Ewing  curve 
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Hamilton  curve  (45)  is  from  sonobuoy 
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Figure  55.  Area  map  of  Stations  4  and  5  in  the  78 

Tagus  Abyssal  Plain  with  bathymetry. 

Base  map  is  by  Rud diman  and  Glover 
(48).  Cores  were  taken  by  Ruddiman  on 
WILKES.  Profiles  C-C 1  and  B-BT  are 
shown  in  Figure  59.  Lines  with  arrows 
show  bomb  runs  for  the  two  stations. 
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records  taken  along  bomb  run  at  Station 

4,  Tagus  Abyssal  Plain. 

Figure  57.  3.5  kHz  subbottom  depth  recorder  80 

records  taken  along  bomb  run  at  Station 

5,  Tagus  Abyssal  Plain. 

Figure  58.  Seismic  profiller  records  near  Stations  81 
4  and  5.  Tracks  B-B 1  and  C-C T  are  shown 
in  Figure  55.  Depths  In  figure  are  un¬ 
corrected  . 


Illustrations 


Figure  59.  Sound  speed  in  cores  from  the  Tagus  82 

Abyssal  Plain.  Locations  of  cores  are 
shown  in  Figure  55.  Figure  is  from 
Ruddiman  and  Glover  (48). 


Figure  60.  Area  map  of  Stations  7  and  8  at  Dragon  83 
Seamount  with  bathymetric  contours. 

Contours  are  based  on  Hayes  et  al .  (55) 
and  are  modified  with  new  data  obtained 
along  tracks  (indicated  by  dots  and 
heavy  solid  lines).  Profiler  records 
and  3.5  kHz  records  along  heavy  solid 
lines  are  shown  in  Figures  61,  62,  and 
63.  The  peak  of  Dragon  Seamount  is 
about  7  nm  south  of  Station  7.  Cores 
V30-193  and  V27-162  and  DSDP  Site  136 
are  indicated. 


Figure  61.  3.5  kHz  subbottom  depth  recorder  record  84 
taken  along  bomb  run  at  Station  B. 

Figure  62.  Seismic  profiler  records  near  Station  85 
8.  Tracks  and  lettered  locations  W,  Wf, 

K,  and  Kf  are  indicated  in  Figure  60, 

Figure  63.  3.5  kHz  subbottora  depth  recorder  record  86 
taken  along  bomb  run  at  Station  7, 
seismic  profiler  record  crossing  bomb 
runs  for  Stations  7  and  8,  and  line 
drawing  of  sea  floor  along  Station  7 
bomb  run  with  no  vertical  exaggeration. 
Track  locations  are  shown  in  Figure  60. 

The  arrow  connecting  the  two  profiles 
are  the  crossing  points. 

Figure  64.  Sound  speed  in  cores  taken  at  Station  7  87 

(SB3)  and  Station  8  (SB4).  Sound  speed 
values  of  the  core  and  bottom  water 
have  been  adjusted  to  laboratory  condi¬ 
tions  at  20°C. 


Figure  65.  Predicted  sediment  sound  speed  profiles  88 
at  Station  8.  The  heavy  line  with  a 
gradient  of  2.33  sec“l  is  the  average 
gradient  determined  at  DSDP 
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Illustrations 


Figure  66.  Area  map  for  Station  9  on  the  lower  89 

continental  rise  south  of  Madeira.  Base 
map  is  from  Uchipi  [59].  Station  8  and 
DSDP  Site  136  are  indicated.  Tracks  of 
the  WILKES  and  KANE  from  Station  8  to 
the  end  of  the  cruise  at  Las  Palmas  are  . 
shown.  Sonobuoy  stations  and  cores  were 
taken  by  Lamont-Doherty  Geological  Ob¬ 
servatory. 

Figure  67.  3.5  kHz  subbottom  depth  recorder  and  90 
seismic  profiler  records  taken  along 
bomb  run  at  Station  9.  Tracks  are  shown 
in  Figure  66. 

Figure  68.  Sound  speed  gradient  determination  at  91 
DSDP  Site  118  using  the  CRIC  method. 

The  Mid-R  reflector,  observed  at  0.48 
sec  two-way  travel  time  on  the  profiler 
record  shown  on  the  left,  corresponds 
to  sandstone  layers  drilled  at  400m. 

Figure  69.  Results  of  Lamont-Doherty  sonobuoy  sta-  92 
tions  off  Northwest  Africa.  Location  of 
stations  are  shown  in  Figure  66. 

Figure  70.  Sediment  sound  speed  profile  determina-  93 
tions  from  wide-angle  reflection  data 
for  Group  W  data.  Calculated  interval 
sound  speed  is  plotted  as  a  function  of 
depth  to  the  midpoint  of  the  layer. 

Sound  speeds  and  depths  correspond  to 
Group  W  data  presented  In  Figure  69. 

Dashed  curve  is  the  least-squares  fit, 
second  degree  curve  through  the  Group  W 
data  set.  Data  near  Station  9  are  indi¬ 
cated  by  squares. 

Figure  71.  Sound  speed  profile  determina tions  from  94 
wide-angle  reflection  data  for  Group  E 
data.  See  explanation  for  Figure  70. 
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Illustrations 


Figure  72.  Sediment  porosity  values  at  DSDP  Site  95 
136.  Gross  descripton  of  lithology  and 
age  of  the  sediment  profile  at  DSDP 
site  is  indicated  on  the  left.  Core 
sampling  was  not  continuous.  Slanted 
solid  and  dashed  straight  line  segments 
represent  the  range  of  porosity  meas¬ 
ured  over  the  sample  length.  Solid 
lines  are  GRAPE  measurements;  dashed 
are  conventional  weight  measurement 
method.  The  percent  rebound  correspond¬ 
ing  to  the  sample  depth  is  indicated. 

Curves  are  from  Hamilton  (83).  See 
text . 

Figure  73.  Water  sound  speed  profiles  at  acoustic  96 
stations.  Source,  receivers,  and  water 
depth  is  indicated  for  the  bomb  run  at 
each  station.  Location  map  is  modified 
from  Fenner  and  Bucca  (87). 
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Symbols  and  Abbreviations 


Co,i 

Sound  speed  in  water  a  or  sediment  layer 
i,  m/sec. 

CR 

Sound  speed  ratio  at  first  interface,  used 
without  subscript  when  it  is  clear  that  it 
derives  from  Rayleigh  or  arrival-time  anal¬ 
ysis.  CR  -  Cse(j^ment/Cwater . 

CRi 

Sound  speed  ratio  at  interface  i,  from 
acoustic  arrival  time  analysis. 

CRi  “  Ci/Ci_i. 

crg,r,t 

CR  determined  from  geological  data  (G), 
from  Rayleigh  coefficient  analysis  of  acous¬ 
tic  data  (R),  or  from  time  of  arrival  analy¬ 
sis  of  acoustic  data  (T). 

CRIC 

Surface  sediment  sound  speed  and  interval 
sound  speed  method  of  determining  g. 

DR 

Density  ratio  at  first  interface,  used  with¬ 
out  subscript  when  it  is  clear  that  the 
acoustically  inferred  value  from  Rayleigh- 
coefficient  analysis  is  intended. 

DR  =  p  sediment /p  water. 

drG,R 

DR  determined  from  geological  data  (G)  or 
from  the  Rayleigh-coefficient  analysis  of 
the  acoustic  data  (R). 

DSDP 

Deep  Sea  Drilling  Project. 

f 

Frequency. 

g 

Sound  speed  gradient  in  a  layer  (positive 
for  sound  speed  increasing  with  depth) 
sec”l . 

h 

Thickness  of  a  layer,  m. 

HF 

High  Frequency. 

IC 

Interval  or  average  sound  speed  of  a  layer. 

i 

Layer  index  increasing  with  depth  starting 
with  the  water  as  layer  zero.  Also  used  for 
the  unit  imaginary  symbol  i  =  (-l)*'*. 
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Symbols  and  Abbreviations 


Wave  number ,  2  7 r/X  • 


Low  Frequency 

Amplitude  of  an  arrival  in  the  unit  impulse 
response  of  the  bottom. 

Slant  range,  i.e.,  total  path  length  from 
source  to  receiver  along  bottom  interacting 
ray  path. 

Plane  wave  reflection  coefficient  at  the  top 
of  a  layer. 


Root  Mean  Square. 

Total  two-way  ray  path-length  through  the 
sediments 


Two-way  ray  path-length  of  nth  arrival 
propagation  path  through  layer  i. 

Plane  wave  transmission  coefficient  at  top 
of  a  layer 

Measured  two-way  travel  time  difference; 
i.e.,  the  time  difference  between  an  obser¬ 
ved  arrival  and  the  first  bottom  reflected 
arrival. 

Two-way  travel  time  elapsed  in  propagation 
of  nth  arrival  through  layer  ±. 

Two-way  travel  time  differences  between 
first  and  nth  arrival. 

Root-mean-square  difference  between  meas¬ 
ured  and  computed  arrival  times  for  a  given 
station  and  acoustically  inferred  model  (A) 
or  geologically  inferred  model  (G). 

Depth  below  water-sediment  interface. 

Plane-wave  attenuation  coefficient,  dB/m. 

Grazing  angle  at  the  water-sediment  in¬ 
terface  (o),  or  at  bottom  of  ith  layer  (i). 


arrival  times  to  obtain  estimates  of 
sediment  sound  speed,  density,  and  at¬ 
tenuation. 

Results  of  the  analyses  indicated  the 
utility  of  the  experimental  design  and 
deconvolution  processing  for  the  inves¬ 
tigation  of  bottom-interaction  phenom¬ 
ena.  Wide-angle  seismic  reflection  tech¬ 
niques  are  usually  limited  to  layers 
more  than  approximately  200  m  thick  in 
deep  water  [39].  The  deconvolution  tech¬ 
niques  employed  here  for  sediment  sound- 
speed  analysis  extend  the  capability  to 
layers  as  thin  as  10  m. 

The  low-frequency  impulse  response  be¬ 
havior  (at  all  but  the  sediment-free 
bottom  site,  station  7)  was  modeled  by 
simple  refracting-  and/or  reflecting- 
layer  models.  The  average  root-mean- 
square  (RMS)  arrival-time  error  (pre¬ 
dicted  arrival  time  versus  observed)  for 
all  the  models  was  3.1  msec,  with  a 
maximum  RMS  error  of  8.6  msec  (Station 
2,  model  6).  These  RMS  errors  are  smal¬ 
ler  than  the  14.2  msec,  3  dB,  full-width 
pulse  duration  for  the  low-frequency 
impulse  responses.  The  sediment  sound- 
speed  structure  at  shallow  depths,  as 
determined  from  two-layer  modeling  of 
high-frequency  impulse  responses,  was 
suggestive  of  higher  sound  speeds  than 
would  be  expected  on  the  basis  of  geo¬ 
logical  information. 

Geological  analyses  of  available  envi¬ 
ronmental  data  were  conducted  to  develop 
estimates  of  sediment  sound-speed,  den¬ 
sity  and  attenuation  profiles.  Geologi¬ 
cal  data  were  obtained  from  (1)  cores, 
(2)  seismic-profiler  records,  (3)  3.5 

kHz  profiler  records,  (4)  wide-angle  re¬ 
flection  profiling  (sonobuoy  data)  and 
(5)  Deep  Sea  Drilling  Project  (DSDP) 
holes  drilled  in  locations  with  sedi¬ 
ments  analogous  to  the  site  of  interest. 
Whenever  possible  the  sound-speed  pro¬ 
file  was  specified  in  terms  of  an  effec¬ 
tive  constant  gradient  and  determined 
from  DSDP  sites  with  a  sedimentary 
structure  similar  to  that  of  the  station 
as  determined  from  other  geological  con¬ 
siderations.  Other  methods  for  deriving 


the  gradient  are  also  discussed  in 
Appendix  E. 

Sediment  density  and  surficial  sound- 
speed  were  generally  measured  from  core 
samples.  Attenuation  estimates  were 
based  on  core  samples  and  relationships 
presented  by  Hamilton  [45].  In  addition, 
other  information  relevant  to  the  acous¬ 
tic  properties  are  summarized  including: 
water  depth,  bottom-water  sound-speed, 
sediment  thickness,  and  reflectors  with¬ 
in  the  sediments.  General  discussions  of 
geological  features  affecting  acoustic 
properties  at  several  of  the  stations 
are  presented  in  Appendix  D. 

Geologically  inferred  predictions  of 
sediment  surface  sound-speed  and  density 
were  in  agreement  with  the  acoustical 
estimates  (DRR,  and  CRR) ,  but  the 
overall  sound-speed  model  was  inade¬ 
quate,  as  judged  by  errors  in  the  pre¬ 
dicted  arrival  times  (RMS  t  q  =  50 

msec).  Both  the  acoustical  and  geologi¬ 
cal  estimates  of  attenuation  were  sub¬ 
ject  to  large  uncertainties,  and  the 
geological  estimate  was  consistently 
higher  than  the  acoustical  estimate. 

The  acoustical  model  results  of  this 
study  suggest  generalization  to  other 
deep  ocean  abyssal-plain  and  abyssal- 
hills  bottom  sites.  However,  since  the 
geological  data  base  was  of  the  type 
that  is  generally  available,  the  diffi¬ 
culties  with  geological  model  prediction 
may  also  be  encountered  at  other  sites. 

Recommendations  for  improving  accuracy 
of  geologically  inferred  models  include: 
(1)  recording  and  analyzing  full  dynamic 
range  of  amplitude  information  in  nor¬ 
mal-incidence  profiler  data;  (2)  pulse 
compression  and  dereverberation  pro¬ 
cessing  (deconvolution)  of  profiler  and 
wide-angle-reflectivity  data;  and,  (3) 
improved  sampling  in  DSDP  cores  of  the 
top  few  hundred  meters. 

In  the  organization  of  the  report,  many 
details  have  been  included  which  the 
reader  may  wish  to  skip  on  first  read¬ 
ing.  Section  II  is  a  brief  historical 


2 


review  of  progress  in  the  study  of  the 
acoustic  bottom-interaction  problem. 
Section  III  describes  the  usual  two-ship 
experiments  often  used  for  collecting 
bottom-reflectivity  data.  Section  IV. A 
may  be  skipped  by  the  reader  familiar 
with  deconvolution  processing.  In  Sec¬ 
tion  IV. B,  the  geoacoustic  models  are 
described  together  with  the  analysis 
procedures  for  determining  model  param¬ 
eters  from  the  measured  data.  Section 
IV.  C  gives  the  detailed  acoustic  re¬ 
sults,  and  Section  V  gives  the  detailed 
geological  results  for  each  station.  A 
table  summarizing  the  most  significant 
results  is  located  at  the  end  of  these 
last  two  sections.  A  detailed  comparison 
of  the  acoustical  and  geological  results 
is  discussed  in  Section  VI,  and  conden¬ 
sed  statements  of  conclusions  and  recom¬ 
mendations  are  given  in  Sections  VII  and 
VIII.  References  and  illustrations  con¬ 
stitute  Sections  IX  and  X  followed  by 
several  appendices  with  more  detailed 
information.  Appendix  A  begins  with  a 
summary  of  formulas  used  in  the  acousti¬ 
cal  analyses.  A  fuller  discussion  of  the 
characteristics  of  the  acoustic  data  is 
reserved  for  Appendix  B.  The  available 
geological  and  geophysical  information 
relevant  to  acoustic  properties  at  each 
site  is  discussed  in  Appendix  C.  Addi¬ 
tional  information  explaining  the  ef¬ 
fects  of  marine  geological  processes  on 
the  bottom  acoustic  response  and  de¬ 
scribing  presently  available  techniques 
for  measurement  and  prediction  of  geo¬ 
acoustic  parameters  is  included  in  Ap¬ 
pendices  D  and  E.  The  paper  concludes 
with  Appendix  F,  which  gives  a  brief 
discussion  of  the  structure  of  the 
sound-speed  profiles  at  each  of  the  sta¬ 
tions.  It  is  suggested  that  a  rapid 
reading  of  the  paper  include  Sections 
IV. B,  IV. C,  V,  VI,  VII,  and  VIII. 

II.  Background 

One  of  the  earliest  theories  to  be  ap¬ 
plied  to  the  bottom-interaction  problem 
was  that  of  Lord  Rayleigh  [1]  during  the 
1890* s.  Rayleigh  derived  the  reflection 
coefficient  for  plane  waves  incident  at 
arbitrary  angles  on  the  boundary  separa¬ 
ting  two  homogeneous  fluid  half-spaces. 


According  to  Sergeev  [2],  continuous 
ultrasonic  echo-sounding  measurements 
were  made  from  a  moving  ship  during  the 
period  1918-1920  to  evaluate  bottom  to¬ 
pography  and  to  determine  whether  the 
bottom  was  hard  or  soft.  Liebermann  [3] 
recorded  24  kHz  signals  that  were  recei¬ 
ved  via  direct  and  bottom-reflection 
paths  at  a  9°  grazing  angle.  By  meas¬ 
uring  the  maximum  and  minimum  intensi¬ 
ties  as  the  source  height  was  varied,  he 
was  able  to  determine  the  bottom-reflec¬ 
tion  coefficient.  He  concluded  that  the 
reflectivity  depended  on  sediment  grain 
size  and  sorting.  By  1949  Sergeev  also 
began  to  utilize  echo  pulse  amplitude 
and  waveform  information  to  correlate 
with  bottom  type. 

During  seismic  refraction  studies,  Hill 
[4]  observed  arrivals  which  he  attrib¬ 
uted  to  refraction  within  a  sound-speed 
gradient  in  the  bottom.  Several  of  the 
later  studies  of  the  sound  speed  gradi¬ 
ent  are  discussed  in  Appendix  E. 

Hamilton  et  al.  [5]  measured  in  situ  and 
laboratory  values  of  sediment  density, 
porosity,  grain  size,  sound-speed,  and 
attenuation  at  100  kHz.  They  discussed 
the  honeycomb  and  honeycomb  flocculant 
sediment  structures  formed  by  deposition 
of  small  particles,  the  bridging  effect 
of  platey  minerals,  and  the  attendant 
relationship  between  sediment  porosity 
and  grain  size.  From  measurements  of 
sediment  density  and  sound-speed,  nor¬ 
mal-incidence  bottom-loss  values  were 
calculated  and  found  to  be  in  good 
agreement  with  bottom-loss  measurements 
taken  by  Mackenzie  [6].  Breslau’s  [7] 
work  on  12  kHz  normal-incidence  bottom- 
reflectivity  data  statistically  demon¬ 
strated  highly  significant  correlations 
between  bottom  loss  and:  (1)  porosity; 

(2)  amount  of  sand,  silt  or  clay;  and 

(3)  sorting. 

In  the  work  of  Mackenzie  [6],  the  fluid 
half-space  model  of  the  bottom  was 
brought  to  full  quantitative  develop¬ 
ment.  In  his  measurements  of  bottom 
reflectivity  at  grazing  angles  from  12 
to  90°  and  at  frequencies  of  4,  7.5, 

and  16  kHz,  Mackenzie  found  that  the 
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Rayleigh-Morse  theory  (in  which  the 
Rayleigh  equation  is  modified  for  bottom 
absorption  by  inclusion  of  a  complex 
sound-speed)  provided  a  good  theoretical 
description  of  the  data*  He  thus  demon¬ 
strated  that  an  older  theory  of  normally 
reacting  surfaces  [8]  was  inappropriate 
for  bottom  sediments.  In  addition,  the 
theory  fit  the  data  best  when  attenua¬ 
tion  coefficients  were  extrapolated  from 
measurements  at  30  kHz  according  to  a 
linear  or  first— power  dependence  on  fre- 
que  ncy . 

A  flurry  of  research  activity  (Jones  et 
al.  [9],  Barnard  et  al.  [10]  in  1964, 
Brown  and  Ricard  [11],  Bucker  et  al. 
[12],  Menotti  et  al.  [13],  Cole  [14]  in 
1965,  Smith  and  Li  [15]  in  1966,  and 
Winokur  and  Bohn  [16]  in  1968)  was  con¬ 
cerned  with  the  next  logical  step  for 
improving  bottom-loss  models:  multi¬ 
layered  fluid  and  viscoelastic  models. 
These  models  permitted  more  accurate 
bottom-loss  curve  fitting  and  explained 
some  of  the  variability  observed  as  a 
function  of  grazing  angle.  An  important 
step  was  taken  by  Morris  [17],  who  dem¬ 
onstrated  that  better  data  fits  were 
obtained  after  incorporating  sound-speed 
gradients.  In  particular,  they  provided 
an  explanation  for  low  losses  observed 
at  low  grazing  angles. 

In  1970,  Hastrup  [18]  and  Hovem  [19] 
introduced  the  technique  of  deconvolu¬ 
tion  for  processing  bottom— ref lectivity 
data.  Deconvolution  removes  bubble-pulse 
reflections  from  bottom-reflected  traces 
generated  by  explosive  charges,  and 
increases  the  resolution  or  spikiness  of 
primary  reflected  pulses.  Sediment- 
reflecting  horizons  are  more  easily 
identified,  thereby  aiding  in  the  syn¬ 
thesis  of  layered  models.  Hastrup  com¬ 
puted  theoretical  impulse  responses  from 
sediment  core  measurements  and  found 
that  they  compared  favorably  with  the 
experimentally  measured  impulse  respon¬ 
ses.  Frequency  domain  deconvolution 
techniques  were  also  employed  by  Van 
der  Veen  and  Santaniello  [20]  and  refin¬ 
ed  by  Lackoff  and  LeBlanc  [21],  Dicus 
[22],  and  Herstein  [23].  Homomorphic 


deconvolution  was  demonstrated  by  Stoffa 
et  al.  [24,  25]. 

Interest  in  the  effects  of  a  sound-speed 
gradient  continued.  Hanna  [26]  showed 
how  the  bottom  gradient  could  be  esti¬ 
mated  from  the  obeservation  of  inter¬ 
ference  nulls  due  to  superposition  of 
the  reflected  and  refracted  ray  at  a 
bottom-mounted  hydrophone  using  a  CW 
source.  Christensen  et  al.  [27]  ident¬ 
ified  the  refracted  arrival  associated 
with  a  bottom  sound-speed  gradient  by 
comparing  high-pass  and  low-pass  filter¬ 
ed  bottom-reflected  traces.  In  addition, 
they  showed  that  the  refracted  pulse 
arrival  times  for  varying  source-to- 
receiver  ranges  were  predicted  by  ray 
trace  computations  with  a  suitably 
chosen  gradient  value.  Dicus  [28]  meas¬ 
ured  bottom  impulse  responses  for  a 
gradient  bottom  and  demonstrated  precise 
quantitative  modeling  of  the  impulse 
responses  and  associated  bottom-loss 
functions.  Santaniello  et  al .  [29]  dis¬ 

cussed  the  focusing  of  sound  rays  in  a 
gradient  bottom,  the  contribution  of 
lateral  (head)  waves,  and  the  effect  of 
different  measurement  geometries  in  con¬ 
nection  with  bottom-loss  data  that 
showed  significant  negative  losses  over 
a  restricted  grazing  angle  range. 
Stickler  [30,  31]  pointed  out  the  inade¬ 
quacy  of  analyzing  data  with  the  as¬ 
sumption  of  incident  plane  waves,  and 
showed  numerical  examples  of  the  shift 
in  apparent  critical  angle  and  occur¬ 
rence  of  negative  losses  associated  with 
different  measurement  geometries.  For  an 
isospeed  bottom,  he  presented  a  new 
method  for  measuring  bottom  sound  speed 
with  bottom-mounted  sensors  in  the  case 
of  a  point  source  [32]. 

III.  Description  of  Experiment 

The  acoustic  data  were  obtained  during  a 
two-ship  operation  conducted  by  the 
Naval  Oceanographic  Office  (NAVOCEANO) 
in  the  northeast  Atlantic  from  25  Sep¬ 
tember  to  17  October  1973  between  Barry, 
Wales,  and  Las  Palmas,  Canary  Islands. 

The  locations  of  the  acoustic  stations 
are  shown  in  Figure  1  and  lisced  in 
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Table  1.  Station  locations  and  times 


Station  Province 


Time/Date  (1973)  Latitude/Longitude 


2 

Northern  Iberian 

Abyssal  Plain 

0110/29  Sep 
0503/29  Sep 

42°49.0'N/13°43.5,W 
42°55 .3 'N/13°01 .0 'W 

4 

Southeast  Tagus 

Abyssal  Plain 

2255/1  Oct 
0250/2  Oct 

37°10.7'N/llo15.0'W 

37°30.0,N/11°38.9,W 

5 

Southwest  Tagus 

Abyssal  Plain 

2246/3  Oct 
0224/4  Oct 

37°7.4'N/12°25.4,W 

37°24.2'N/11°59.1'W 

7 

Dragon  Seamount 

1900/7  Oct 
2118/7  Oct 

34°59.9,N/16°28.5'W 

34°40.0,N/16°27.8'W 

8 

Dragon  Seamount  Plateau 

1111/8  Oct 
1349/8  Oct 

34°55.4'N/16°37.2,W 

34°39.5'N/16°52.5'W 

9 

Lower  Continental  Rise 
South  of  Madeira 

0635/11  Oct 
0930/11  Oct 

29°0.0,N/20°49.5,W 

29°11.4,N/20°27.5,W 

Table  1.  Of  the  six  stations,  four  are 
in  deep  ocean  basins  (2,4,5,  and  9),  one 
is  on  an  extremely  rough  basement  out¬ 
crop,  a  seamount  (Station  7),  and  one  is 
on  a  shallow,  thickly  sedimented  feature 
distant  from  shore  (Station  8). 

Bottom-reflectivity  measurements  were 
obtained  using  the  experimental  geometry 
shown  in  Figure  2.  USNS  WILKES  was  an¬ 
chored  during  the  experiment  in  order  to 
receive  signals  at  a  deep  hydrophone 
from  SUS  charges  launched  from  USNS  KANE 
while  traveling  away  from  WILKES  along  a 
track  referred  to  as  the  bomb  run.  Ex¬ 
plosive  source  charges  (MK  61  SUS,  0.8 
kgm  TNT)  were  detonated  at  244  m  depth. 
Signals  were  received  on  an  omnidirec¬ 
tional  hydrophone  at  approximate  depths 
of  2500  m  for  stations  2,  4,  and  5;  1700 
m  for  stations  8  and  9;  and  244  m  for 
station  7.  The  depth  of  the  source  was 
chosen  as  a  compromise  between  greater 
low-frequency  energy  at  shallow  depths 
and  greater  multipath  resolution  at 
greater  depths.  Receiver  depths  were 
chosen  to  optimize  reception  of  direct- 
water-path  signals  except  at  station  7. 
At  station  7  a  shallow  receiver  was 
chosen  to  inhibit  multipath  separation 
and  allow  integration  over  all  surface- 
reflected  multipaths.  Multipath  resolu¬ 
tion  at  station  7  was  impossible  for  any 


experimental  geometry  because  of  the 
long  response  time  of  the  reverberant 
bottom. 

Charges  were  detonated  at  spacings  rang¬ 
ing  from  2  to  10  minutes;  the  interval 
was  chosen  to  give  a  grazing-angle 
increment  of  about  3°  for  each  geo¬ 
metry.  The  initial  charge  was  detonated 
at  a  distance  from  the  receiver,  which 
provided  a  bottom-reflected  ray  path  at 
near-normal  grazing  angle.  The  minimum 
grazing  angle  was  obtained  at  the  maxi¬ 
mum  range  near  the  first  convergence 
zone.  The  actual  grazing  angles  were 
computed  during  data  analysis  by  using  a 
ray-trace  program  [33]  incorporating 
sound-speed  profiles  obtained  from  SVSTD 
(sound  velocity/ salinity/ temperature/ 
depth)  measurements  taken  either  immedi¬ 
ately  before  or  after  the  experiment. 

The  hydrophone  used  was  a  Clevite  Model 
CH-1A.  The  system  response  was  flat 
within  1  dB  below  2.2  kHz  and  4  dB  down 
at  4  kHz.  All  measurements  utilized  the 
self-calibration  property  of  the  syn¬ 
thetic-replica  deconvolution  technique 
[22]. 

The  received  signals  were  fed  into  a 
high  impedance  (500  megohm  in  parallel 
with  100  pfd)  Ithaco  453  amplifier. 
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The  amplified  signal  was  monitored  and 
adjusted  in  5  dB  increments  to  keep  the 
input  signal  to  the  recorder  at  about  -6 
dB  re/1  v.  Two  identical  Lockheed  7- 
track  recorders  with  a  dynamic  range  of 
35  dB  were  used  with  a  10  dB  gain  sep¬ 
aration*  The  data  were  recorded  in  both 
FM  and  direct  modes  at  7-1/2  ips.  The  FM 
response  was  0-2500  Hz;  the  direct  res¬ 
ponse  was  100-25,000  Hz.  Two  other  chan¬ 
nels  were  used:  one  for  verbal  annota¬ 
tions  and  the  other  to  record  a  Systron 
Donner  8150  time-code  generator  signal. 
The  remaining  channel  was  unrecorded. 

IV.  Acoustical  Data  Analysis 

A.  Digital  Processing 
1.  Deconvolution 

The  data  were  processed  digitally  to 
obtain  the  ocean-bottom  impulse  res¬ 
ponse.  Once  accomplished,  little  addi¬ 
tional  effort  was  required  to  obtain 
bottom-loss  results  digitally.  Prior  to 
digitization,  the  analog  data  tapes  were 
pre-edited  by  listening  to  the  data  and 
voice  tracks  simultaneously  and  decoding 
the  recorded  time  code.  In  this  way  ab¬ 
solute  start  and  stop  times  were  deter¬ 
mined,  which  bracketed  the  desired  data 
for  digitizing.  The  analog  tapes  were 
then  digitized  at  the  computer  center  of 
the  Naval  Ship  Research  and  Development 
Center,  Carderock,  Maryland.  The  start 
times  determined  during  the  pre-editing 
were  preset  in  a  time  decoder  for  ini¬ 
tiating  digitization.  Following  digiti¬ 
zation,  the  newly  created  digital  tapes 
were  processed  on  a  Control  Data  Corpor¬ 
ation  Model  3800  computer  at  the  Naval 
Research  Laboratory,  Washington,  D.C. 

Impulse  responses  were  calculated  using 
a  synthetic-deconvolution  technique 
[22}.  With  the  synthetic  method,  a  rep¬ 
lica  of  the  acoustic  signal  incident  at 
the  ocean  bottom  is  generated  by  a  time 
expansion/contraction  transformation  of 
a  directly  measured  source  signal.  The 
direct  signal  is  measured  only  once, 
while  the  amount  of  expansion  or  con¬ 
traction  is  varied  on  a  shot-by-shot 
basis.  The  expansion  factor  adjusts  the 


bubble-pulse  period  of  the  direct  signal 
to  that  of  the  particular  shot  being 
processed.  The  bubble-pulse  period  is 
measured  as  the  time  between  the  peak  of 
the  shock  pulse  and  the  peak  of  the 
first  bubble  pulse. 

Bubble-pulse  periods  were  generally 
measured  to  an  accuracy  of  +0.1%  using 
the  inverse  filter  technique  on  measured 
direct-path  signals.  At  the  longer  rang¬ 
es  corresponding  to  grazing  angles  of 
less  than  approximately  20°  9  suitable 
direct-path  signals  were  not  available 
because  of  refractive  shadowing,  distor¬ 
tion,  or  excessive  noise.  In  most  of 
these  instances,  the  bottom-reflected 
signal  was  used  for  bubble-pulse-period 
determination.  Primary  reflected  pulses 
were  separated  in  time  from  the  bubble 
pulses  and  hence  did  not  interfere  with 
them.  At  Station  7,  bottom-reflected 
signals  could  not  be  used  with  the  in- 
verse-filter  technique  because  of  inter¬ 
ference  between  bubble  pulses  and  pri¬ 
mary  reflected  pulses.  In  this  case,  the 
cepstrum  technique  was  applied. 

The  deconvolutions  were  accomplished  in 
the  frequency  domain  via  discrete  finite 
Fourier  transforms.  The  effects  of  noise 
were  minimized  and  the  impulse  responses 
stabilized  by  linear,  minimum-mean- 
square-error  estimation.  Geometric 
spreading  and  absorption  losses  incurred 
along  the  water  paths  were  removed  in 
order  that  the  impulse  responses  repre¬ 
sent  only  the  contribution  of  the  bottom 
interaction. 

An  eigen-ray  computation  was  made  for 
each  station  using  the  Naval  Research 
Laboratory  program,  Trimain  [33],  with 
modification  for  printing  bottom  graz¬ 
ing  angles  and  ray  slant  ranges.  For 
each  station,  Trimain  provided  the 
transmission-loss  and  slant-range  infor¬ 
mation  needed  to  make  the  geometric¬ 
spreading  and  absorption-loss  correc¬ 
tions  . 

A  zero-phase  filter  with  Gaussian- 
shaped  amplitude  spectrum  was  applied  to 
the  impulse  responses  to  obtain  a  desir¬ 
able  pulse  shape  free  of  side-lobes. 
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The  3  dB  points  of  the  filter  response 
were  set  to  43  and  87  Hz  for  the  low- 
frequency  traces  and  to  2*1  kHz  for  the 
low-pass  filtered  broadband  traces.  For 
the  low  frequency  band,  the  lower  cutoff 
frequency  was  chosen  to  exclude  low- 
frequency  noise  and  the  upper  cutoff 
frequency  was  chosen  to  limit  the  com¬ 
plication  of  additional  layering  resol¬ 
ved  at  greater  band-widths.  The  cutoff 
frequency  for  the  broadband  or  high- 
frequency  band  was  chosen  to  eliminate 
high-frequency  noise  imposed  by  the 
high-frequency  roll-off  of  the  tape- 
recorder  frequency-modulation  circuitry. 
The  impulse-response  amplitudes  were 
normalized  to  the  output  corresponding 
to  a  unit  amplitude  band-limited  input. 

2.  Bottom  Loss  Computations 

Bottom-loss  can  be  defined  as  the  ratio 
of  the  plane  wave  energy  incident  on  the 
ocean  bottom  to  the  reflected  energy 
expressed  in  decibels  as  a  function  of 
frequency  and  grazing  angle.  In  calcula¬ 
ting  bottom-loss  from  the  experimental 
data,  geometric  spreading  and  absorption 
losses  must  be  calculated  from  the  ray 
path  associated  with  specular  reflec¬ 
tion.  In  the  case  of  a  rough  bottom, 
there  are  many  possible  bottom-reflec¬ 
ting  rays  connecting  the  source  and 
receiver  for  any  geometry.  For  Stations 
7  and  8  the  bathymetry  was  variable 
about  a  mean  slope.  The  geometric 
spreading  and  absorption  losses  were 
then  computed  based  on  the  ray  path 
associated  with  specular  reflection  from 
a  plane  inclined  at  the  mean  slope. 

For  all  the  stations  except  7,  the 
single  bottom-interacting  signals  were 
easily  time-separated  from  water  surface 
reflected  multiples.  At  Station  7,  time 
separation  was  impossible  both  before 
and  after  deconvolution.  In  this  case, 
however,  the  four  main  arrivals  associ¬ 
ated  with  a  single  bottom  bounce  could 
be  separated  from  the  next  group  of 
arrivals  associated  with  two  bottom 
bounces.  Therefore,  the  first  group  of 
arrivals  were  isolated  and  a  6  db  cor¬ 
rection  applied  to  remove  the  average 
effect  of  the  four  arrivals. 


Bottom-loss  versus  grazing-angle  func¬ 
tions  were  computed  for  each  station  in 
one-third-octave  frequency  bands.  The 
squared  magnitudes  of  the  Fourier  coef¬ 
ficients  of  direct  and  bottom-reflected 
signals  were  summed  separately  in  one- 
third-octave  frequency  bands.  Each  sum 
is  proportional  to  the  respective  acous¬ 
tic  energy.  After  correction  for  water- 
column  geometric  spreading  and  absorp¬ 
tion  losses,  and  correction  of  the  bot¬ 
tom-reflected  energy  for  error'  due  to 
noise,  bottom-loss  was  computed  as  the 
ratio  of  direct  to  bottom-reflected 
energy  expressed  in  decibels.  No  cor¬ 
rections  were  made  for  noise  in  the 
direct  signal  because  the  signal-to- 
noise  ratio  was  much  higher  than  that 
for  the  bottom-reflected  signals. 

The  noise  was  assumed  to  be  colored 
Gaussian  noise  having  an  approximately 
flat  spectrum  within  each  one-third- 
octave  frequency  band.  The  energy  within 
a  one-third-octave  band  is,  then,  the 
sum  of  squares  of  non-zero  mean  Gaussian 
variates  which  have  equal  variances.  If 
there  were  no  signal  and  only  noise, 
then  the  variates  would  be  zero  mean  and 
the  one-third-octave  band  energy  would 
be  distributed  as  the  familiar  Chi- 
Square  variate.  With  signal  energy,  the 
probability  distribution  function  be¬ 
comes  a  noncentral  Chi-Square  distribu¬ 
tion. 

At  low  frequencies,  a  lengthy  sample  of 
noise  is  required  to  allow  for  spectral 
smoothing.  This  step  was  prohibitive 
because  of  machine  computation  time.  As 
an  alternative,  from  10  to  20  short 
noise  samples  already  taken  for  the  de- 
convolution  noise  filtering  were  ensem¬ 
ble  averaged  to  obtain  a  spectrum.  The 
amplitude  of  the  relative  noise  spectrum 
was  adjusted  for  each  shot  by  measuring 
the  energy  in  a  short  noise  sample  that 
was  taken  ahead  of  each  shot. 

The  error-bar  curves  that  bracket  the 
bottom-loss  data  correspond  to  three 
standard  deviations.  Although  the  non¬ 
normal  distribution  precludes  simple 
identification  of  confidence  limits, 
error  bars  corresponding  to  three 
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standard  deviations  should  give  good  or¬ 
der  of  magnitude  estimates  of  the  error 
with  approximately  the  same  confidence 
limit  (97%)  as  that  for  a  normal  distri¬ 
bution. 

B.  Geoacousticai  Modeling 
L  Models 

Three  single-layer  and  four  double-layer 
sediment  models  with  reflecting  or  re¬ 
fracting  layers,  with  or  without  an 
underlying  half -space,  are  employed  to 
describe  the  observed  impulse  responses. 
In  Figures  3-9  the  ray  diagrams  and  the 
sound-speed  profiles  corresponding  to 
the  models  are  shown.  The  rays  numbered 
from  1  to  3  correspond  to  the  first 
three  arrivals  in  order  of  increasing 
arrival  time  arising  from  reflections 
from  layer  interfaces  or  turning  rays 
within  a  layer.  All  the  hypothesized 
models  are  particular  cases  of  total 
internal  reflection;  that  is,  sound  rays 
are  totally  returned  to  the  water  column 
by  either  refraction  or  reflection. 

The  first  model  is  simply  a  constant 
sound-speed  layer  overlying  a  half¬ 
space  with  reflection  coefficient  R  = 
exp(i7T/2).  The  second  model  is  the  same, 
except  for  an  additional  constant  sound- 
speed  layer.  The  third  model  is  a  layer 
with  a  sound-speed  gradient  and  an 
assumed  thickness  several  wavelengths 
greater  than  the  deepest  turning  point. 
The  sound-speed  profile  in  the  sediment 
is  assumed  to  be  the  well-known,  pseudo- 
linear  profile  [17]  with  interface  grad¬ 
ient  g.  For  the  pseudo-linear  profile 
the  sound  speed  is  given  by 

C(z)  =  C-p  ( 1— 2g (z— zp) /Cp)  (1) 

where  C(z)  is  the  sound  speed  at  depth 
z,  Zp  is  the  depth  at  which  the 
pseudo-linear  profile  begins,  Cp  is 
the  sound-speed  at  depth  Zp,  and  g  is 
the  sound-speed  gradient  at  depth  zp. 

The  fourth  model  is  a  layer  with  a  grad¬ 
ient  in  which  deeper  penetrating  rays 


are  intercepted  by  reflection  from  half¬ 
space  with  reflection  coefficient  R  = 
exp(i7r/2).  Model  5  is  a  constant  sound- 
speed  layer  overlying  a  half-space  with 
a  gradient.  Model  6  consists  of  a  layer 
with  a  gradient  overlying  a  half-space 
with  a  gradient.  Model  7  is  a  layer  with 
a  gradient  above  a  constant  sound-speed 
layer  overlying  a  half-space  with 
reflection  coefficient  R  =  exp(i7r/2). 

In  each  case,  the  model  is  applied  over 
a  restricted  grazing-angle  range  from  a 
low  grazing  angle  consistent  with  a 
high-frequency  asymptotic  or  ray-theory 
approximation,  to  a  high  grazing-angle 
at  which  the  model  fails  to  fit  the 
data. 

When  a  pseudo-linear  sound-speed  half¬ 
space  is  involved,  it  is  assumed  that 
the  applicable  grazing-angle  range 
limits  the  turning  points  to  well  above 
the  critical  depth.  The  critical  depth 
is  the  depth  at  which  the  sound  speed 
approaches  infinity.  From  (1)  the 
critical  depth  is  given  by 

zc  =  zx  +  ^x/2g  (2) 

The  pseudo-linear  half-space  is  a  physi¬ 
cally  unrealistic  model  because  of  the 
sound-speed  singularity  [34].  However, 
the  problem  is  effectively  eliminated 
when  the  critical  depth  occurs  in  the 
region  of  "hidden  depths"  [35],  i.e.,  a 
few  wavelengths  below  the  turning  depth. 
This  condition  is  easily  satisfied  for 
sufficiently  small  grazing  angles. 

It  should  be  noted  that  some  of  the 
models  actually  revert  to  one  of  the 
simpler  models  over  some  portion  of  the 
lower  angular  range.  In  Model  2  at 
grazing-angles  below  the  second-layer 
critical  angle  (if  there  is  one),  the 
reflection  process  is  the  same  as  for 
Model  1.  Model  4  reverts  to  Model  3  when 
the  turning  point  of  the  refracted  ray 
is  above  the  bottom  of  the  layer.  Model 
5  reverts  to  Model  1  at  low  angles  if 
there  is  a  critical  angle  for  the  second 
layer.  Similarly,  Models  6  and  7  become 


8 


Model  4  at  low  angles  if  there  is  a 
critical  angle  for  the  second  layer. 

2.  Parameter  Extraction 

The  models  and  the  analyses  which  follow 
are  based  on  plane  waves;  i.e.,  the  geo¬ 
metry  is  specified  by  the  grazing  angle 
independent  of  the  source  and  receiver 
positions.  The  plane  wave  assumption  is 
an  approximation  that  greatly  simplifies 
the  problem.  However,  care  must  be  exer¬ 
cised  to  limit  its  application  to  geo¬ 
metries  for  which  the  approximation  is 
valid.  Rudnick  [36]  pointed  out  that 
plane  waves  never  really  exist  in  any 
real  measurement.  More  recently,  Stick¬ 
ler  [31]  has  shown  numerical  examples  of 
errors  induced  by  the  plane-wave  assump¬ 
tion  for  the  case  of  reflection  from  a 
half— space  with  a  sound— speed  gradient . 
In  particular,  the  observed  critical 
angle  may  be  shifted  significantly  from 
the  plane-wave  value . 

In  spite  of  these  difficulties,  it  can 
be  shown  that  the  plane— wave  approxima¬ 
tion  may  be  applied  under  certain  re¬ 
strictions  regarding  the  measurement 
geometry  and  the  wavelength.  The  re¬ 
flected  field  of  a  point-source  spher¬ 
ical  wave  interacting  with  a  horizontal¬ 
ly  stratified  bottom  can  be  treated  by 
the  method  of  steepest  descents  in  order 
to  give  an  asymptotic  approximation 
valid  at  high  frequencies  and  at  large 
distances.  The  calculation  shows  that 
the  point-source  reflected  field  is  just 
a  spherical  wave  radiating  from  the 
source-point  image  multiplied  by  the 
plane  wave  reflection  coefficient,  plus 
a  correction  term  which  may  be  kept  as 
small  as  desired  by  choosing  sufficient¬ 
ly  large  values  of  kR  (where  k  =  2n/Aand 
R  is  distance  from  the  receiver  to  the 
source-point  image)  and  by  choosing 
grazing-angles  that  are  not  too  near 
critical.  Details  of  the  calculation  for 
the  case  of  a  constant-sound-speed  half¬ 
space  may  be  found  in  Officer  [37]  and 
Brekovskikh  [38].  The  case  of  a  reflect¬ 
ing  layer  and  a  pseudo-linear-sound- 
speed  halfspace  can  be  treated  in  a 
similar  manner.  The  results  of  these 


calculations  are  summarized  in  Appendix 
A,  Tables  5  and  6. 

In  the  plane-wave-approximation  model, 
parameters  may  be  extracted  from  the 
data  by  analysis  of  arrival  times  and 
arrival-pulse  amplitudes.  Formulas  for 
fitting  arrival  times  are  given  for  each 
the  models  in  Table  7.  The  arrival-time 
analysis  yields  an  estimate  of  the  sedi¬ 
ment  sound-speed  at  the  water-sediment 
interface,  as  well  as  sound-speed  param¬ 
eters  and  layer  thickness  for  succeed¬ 
ing  layers.  First,  arrival-peak  ampli¬ 
tudes  (water-sediment  interface  reflec¬ 
tion  coefficients  R^)  are  used  to  ob¬ 
tain  the  sediment  density  and  a  second 
independent  estimate  of  the  sound-speed 
at  the  water-sediment  interface  by  ap¬ 
plication  of  the  Rayleigh  reflection 
equation.  In  each  case,  certain  func¬ 
tions  of  measured  quantities  are  plotted 
and  a  straight  line  curve  is  fitted. 
Model  parameters  are  determined  from  the 
slope  and  the  intercept  of  the  fitted 
line.  Formulas  giving  the  quantities  to 
be  plotted  and  the  meaning  of  the  slope 
and  the  intercept  parameters  are  found 
in  Table  8. 

Average  attenuation  coefficients  are 
estimated  from  the  peak  amplitudes  of 
the  last  arrival  that  passed  through  a 
turning  point.  When  the  last  arrival 
comes  from  a  reflection,  it  is  assumed 
(as  part  of  the  model  definition)  that 
the  reflection  is  total  with  a  it/ 2  phase 
change  (i.e.,  R  =  exp  iiT/2).  Formulas 
for  the  attenuation  analysis  are  given 
in  Table  9. 

The  parameter-extraction  methods  are 
based  on  assumed  simple  sediment  sound- 
speed  models  that  were  chosen  to  fit 
arrival  times  and  peak  heights  of  the 
processed  impulse  responses.  Justifi¬ 
cation  for  this  model  approach  is  based 
on  the  results  given  in  Reference  28  in 
which  Station  5  data  for  the  present 
study  was  analyzed  in  terms  of  Model  3. 
It  was  shown  that  not  only  did  the  model 
fit  the  arrival  times  observed,  but  the 
entire  set  of  impulse  responses  and 
bottom-loss  curves  were  in  striking 
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agreement  with  theoretical  curves  gener¬ 
ated  by  the  model.  The  agreement  includ¬ 
ed  observed  pulse  shapes  and  some  subtle 
features  in  the  loss  curves.  Further¬ 
more,  the  model  used  in  the  analysis  was 
consistent  with  well-documented  geologi¬ 
cal  evidence  for  the  general  existence 
of  sediment  sound-speed  gradients  [39, 
AO,  41]. 

Some  of  the  data  in  the  present  study 
were  found  to  fit  reflecting-  and  some 
ref ract ing-type  behaviors.  Arrival-time 
behavior  for  pulses  reflected  from  an 
interface  differs  markedly  from  that  of 
refracted  pulses  within  a  layer  with  a 
sound— speed  gradient •  In  those  cases 
where  reflecting  behavior  was  observed, 
the  arrival-time  fit  was  even  better 
than  the  average,  as  judged  by  the  RMS 
travel-time  error  t^.  Thus,  because  of 
the  success  in  modeling  the  refracting 
case  when  a  good  time-of-arrival  fit  was 
obtained,  sediment  parameters  were  also 
derived  from  good  time-of-arrival  fits 
in  reflecting  cases.  The  seven  models 
analyzed  included  combinations  of 
reflecting—  and  ref racting-type  layers 
and  half-spaces., 

C.  Acoustical  Results 

Except  for  Station  7,  the  geometry  of 
the  experiment  permitted  time  separa¬ 
tion  of  the  bottom-interacting  arrival 
from  the  subsequent  water-surface- 
reflected  arrivals.  The  traces  were 
truncated  in  time  to  exclude  the  surface 
arrivals  before  the  deconvolution  pro¬ 
cessing  was  applied.  The  effect  of  the 
truncation  can  be  seen,  for  example,  in 
the  Station  2  traces  (see  Appendix  B) 
beginning  at  approximately  200  msec  in 
the  first  eight  high-angle  traces,  and 
then  decreasing  steadily  to  approximate¬ 
ly  100  msec  in  the  lowest-angle  trace. 
The  effect  is  seen  as  an  apparent  sudden 
decrease  in  the  random  noise  level. 

The  effectiveness  of  the  deconvolution 
3-t  all  the  stations  may  be  judged  by  the 
degree  to  which  bubble-pulse  energy  is 
removed  and  by  the  resultant  noise  back¬ 
ground  as  seen  in  the  figures  in 


Appendix  B.  At  each  of  the  stations  ex¬ 
cept  for  Station  7,  the  reflectors  are 
seen  to  correlate  well  from  trace  to 
trace.  The  complexity  of  reflectors  is 
apparent  in  the  high-frequency  traces; 
however,  at  low  frequencies  a  single 
penetrating  arrival  dominates  over  all 
others . 

The  decay  times  of  the  impulse  responses 
of  Stations  2,  4,  5,  8  and  9  are  on  the 
order  of  0.1  sec  and  are  attributed  to 
multiple  reflections  within  the  sediment 
layers.  In  the  geological  discussions 
that  follow  in  later  sections  it  is 
shown  that  at  these  stations  the  bottom 
is  relatively  smooth  and  thickly  sedi¬ 
mented.  That  is,  the  sediment  cover  is 
thick  enough  that  the  dominant  bottom- 
reflected  energy  interacts  via  bottom 
paths  located  wholly  within  the  sedi¬ 
mentary  cover. 

At  Station  7,  on  the  other  hand,  the 
impulse  responses  reveal  no  coherent 
pattern  of  peaks.  Decay  times  of  the 
impulse  responses  are  on  the  order  of  1 
sec  and  are  attributed  to  random  scat¬ 
tering  from  a  rough,  sediment-free  bot¬ 
tom.  This  hypothesis  is  supported  by  the 
fact  that  the  core  sample  showed  no 
penetration. 

The  intent  of  the  data  processing  and 
modeling  was  to  apply  inversion  tech¬ 
niques  to  obtain  information  about  the 
bottom  environment.  The  coherence  and 
regularity  of  the  impulse  responses  for 
Stations  2,  4,  5,  8,  and  9  permitted 

detailed  modeling  and  parameter  extrac¬ 
tion.  This  was  not  possible  for  Station 
7.  Once  effective  models  are  in  hand, 
the  bottom  response  as  a  function  of 
frequency  and  source-receiver  geometry 
may  be  computed  directly  within  the 
limits  of  the  models.  Questions  regard¬ 
ing  surveillance  system  performance  may 
then  be  addressed. 

The  major  results  of  the  analysis  are 
presented  in  Table  2  in  which  the  infer¬ 
red  sediment  parameters  are  given  for 
the  Rayleigh-reflection-coefficient  an¬ 
alysis  and  the  time-of-arrival  analysis. 
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Table  2.  Summary  of  acoustically  inferred  model  parameters 


Sound-speed  of  the  water  Cw  is  in  meters/second,  gradients  g^,  g2  are  in 
sec“l,  thickness  hi  h2  in  meters,  attenuation  <  ^  in  dB/meter/kHz,  and  time 
t^  in  msec.  The  parameters  DR  and  CR^  are  based  on  Rayleigh-reflection-coeffici¬ 
ent  analysis;  the  remaining  parameters  are  derived  from  time-of-arrival  analysis  for 
each  model. 


Station 

Model 

DR 

CRr 

111 

2 

3 

1549.0 

1.4 

0.97 

280 

2 

5 

1549.0 

1.4 

0.97 

12 

2 

6 

1549.0 

1.4 

0.97 

13 

4 

1 

1548.5 

1.3 

1.04 

75 

4 

4 

1548.5 

1.3 

1.04 

76 

4 

3 

1548.5 

1.3 

1.04 

81 

5 

3 

1548.6 

1.6 

0.99 

200 

5 

6 

1548.6 

1.6 

0.99 

37 

8 

1 

1504.8 

1.4 

0.99 

50 

8 

2 

1504.8 

1.4 

0.99 

12 

8 

4 

1504.8 

1.4 

0.99 

49 

9 

1 

1525.5 

1.5 

0.97 

80 

9 

7 

1525.5 

1.5 

0.97 

29 

The  parameters  specified  include:  (1) 
density  ratio  DR,  (2)  sound-speed  ratio 
CR,  (3)  layer  thickness  h,  (4)  sound- 
speed-gradient  parameter  g,  and  (5)  at¬ 
tenuation  constant  k  The  water 

sound-speed  CW  was  from  a  direct  meas¬ 
urement  and  t^  is  the  root-raean-square 
difference  between  the  measured  and  the 
predicted  arrival  times.  The  detailed 
acoustical  results  are  presented  in  Ap¬ 
pendix  B.  Plots  of  impulse  responses  and 
bottom-loss  curves  are  given  together 
with  a  discussion  of  their  features. 
Plots  of  the  straight-line  curve  fits 
are  also  given  for  each  of  the  model 
analyses . 


CR] 

£1 

hi 

CR? 

Bl 

hA 

La. 

1.05 

0.8 

0.05 

6.1 

1.00 

260 

1.05 

0.8 

0.05 

8.4 

0.96 

10.4 

260 

1.00 

0.8 

0.06 

8.6 

1.06 

0.04 

0.5 

1.04 

0.6 

0.06 

0.9 

1.00 

2.1 

0.07 

0.5 

1.03 

0.8 

0.02 

4.8 

0.99 

2.5 

170 

1.00 

0.8 

0.03 

4.0 

1.08 

0.51 

1.9 

1.03 

36 

1.07 

0.51 

1.0 

1.01 

4.3 

0.27 

0.9 

1.06 

0.27 

1.2 

0.97 

4.4 

42 

1.00 

0.30 

1.4 

Average 

=  3.1 

V.  Geological  Data  Analysis 

For  each  station  geoacoustic  properties 
of  the  sea  floor  were  inferred  from  geo¬ 
logical  data  obtained  from  (1)  cores, 
(2)  seismic-profiler  records,  (3)  3.5 

kHz  profiler  records,  (4)  wide-angle 
reflection  profiling  (sonobuoy  data)  and 
(5)  Deep  Sea  Drilling  Project  (DSDP) 
holes  drilled  in  locations  with  sedi¬ 
ments  analogous  to  the  site  of  interest. 
Although  several  methods  for  estimating 
sound-speed  profiles  were  used,  whenever 
possible  an  effective  constant  gradient 
was  determined  from  DSDP  sites  with  a 
sedimentary  structure  similar  to  that  of 
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the  station  as  determined  from  other 
geological  considerations*  An  average 
interval  sound  speed,  IC,  was  calculated 
from  the  acoustic  travel  time  measured 
on  the  seismic-profiler  record  to  a 
reflecting  horizon  correlated  with  a 
drilled  interface  at  a  known  depth.  The 
gradient  was  easily  obtained  from  the 
relationship  involving  IC  and  the  ratio 
of  sediment  surface  sound  speed  to  that 
of  the  bottom  water,  CR.  This  method  of 
determining  a  gradient  is  referred  to  as 
the  CRIC  method  and  is  used  at  Stations 
2,  4,  5,  and  8.  The  gradient  thus  deter¬ 
mined  is  assumed  constant  in  depth  and 
is  the  same  parameter  as  the  "  linear 

gradient”  of  Hamilton  et  al.  [42]. 

A  distinction  should  be  noted  between 
the  definition  of  the  gradient  employed 
in  the  geological  analysis  and  in  the 
acoustical  analysis.  In  the  geological 
analysis  the  sound-speed  gradient  is 
assumed  constant  independent  of  depth 
and  is  denoted  by  the  symbol  g.  In  the 
acoustical  analysis  the  term  gradient 
and  the  symbol  g  refer  to  a  parameter  in 
the  sound-speed  formula  for  the  pseudo- 
linear  profile  (see  Section  IV*B.l).  In 
particular  g  is  the  sound-speed  gradient 
at  the  top  interface  of  the  pseudo- 

linear  layer.  The  sound-speed  gradient 
continues  to  increase  with  depth,  but 
for  layer  thicknesses  and  turning-point 
depths  much  less  than  the  critical 
depth,  the  gradient  is  approximately 
constant — hence  the  term  pseudo-linear 
sound-speed  profile.  In  the  acoustical 
analysis,  layers  were  thin  in  the  above 
sense;  therefore,  g  values  geologically 
inferred  may  reasonably  be  compared  with 
g  values  acoustically  inferred. 

The  environmental  data  from  the  stations 
are  presented  in  Appendix  C.  Each  sta¬ 
tion,  analysis  contains  a  bathymetric  map 
of  the  vicinity  in  addition  to  just  the 
station  site  in  order  to  provide  infor¬ 
mation  useful  in  extrapolating  the 
acoustic  results  beyond  the  immediate 
station  site.  A  general  discussion  of 
geological  features  that  affect  the 
acoustic  properties  at  several  of  the 
stations  is  deferred  to  Appendix  D. 


Table  3  summarizes  the  acoustical  prop¬ 
erties  of  the  sediments  discussed  in 
the  geological  analysis  in  the  appen¬ 
dices.  In  addition,  other  information 
relevant  to  the  acoustic  properties  are 
summarized:  water  depth,  bottom-water 
sound-speed,  sediment  thickness  and 
reflectors  within  the  sediments. 

The  water  depth  is  listed  where  the 
first  and  last  charges  were  dropped, 
except  for  Station  7  where  the  minimum 
depth  along  the  track  is  also  indicated. 
Uncorrected  depths  are  those  read  di¬ 
rectly  off  the  depth  recorder  and  cali¬ 
brated  for  a  water  sound  speed  of  1500 
m/sec.  The  corrected  values  are  derived 
from  the  tables  of  Gold  and  Audet  [43]. 

The  water  sound-speed  listed  in  Table  3 
is  that  at  the  maximum  depth  reached  by 
the  SVSTD.  Actual  bottom-water  sound- 
speed  used  elsewhere  in  this  report  is 
determined  by  extrapolating  this  value 
of  sound-speed  to  the  ocean  bottom  by 
using  Reference  [44].  The  complete  water 
sound-speed  profiles  are  shown  in  Figure 
73. 

The  sediment  thickness  is  given  in  units 
of  two-way  travel  time  to  the  geologic 
basement  as  determined  from  the  sparker- 
profiler  records.  The  value  is  the  aver¬ 
age  along  the  track  of  the  reflectivity 
experiment  as  best  determined  from  re¬ 
cordings  in  the  vicinity.  The  sparker 
could  not  be  operated  during  the  acous¬ 
tic  experiment  because  the  low-frequency 
signal  (50-150  Hz)  of  the  sparker  source 
would  interfere  with  the  acoustic  meas¬ 
urements.  The  actual  thickness  is  a 
function  of  the  sound-speed  profile; 
however,  for  purposes  of  a  rough  esti¬ 
mate,  the  values  list  approximate  thick¬ 
ness  in  kilometers  (assumed  sound-speed 
of  2000  m/ sec) . 

The  depth  of  the  LF  (low  frequency) 
reflectors  in  two-way  travel  time  are 
read  from  the  sparker-prof iler  records 
at  the  closest  point  to  the  acoustic 
station  that  the  sparker  was  operating. 
For  values  of  less  than  0.5  seconds,  a 
rough  estimate  of  the  depth  can  be  made 
using  an  assumed  sediment  sound  speed  of 
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Table  3,  Geologic  data  and  inferred  acoustic  parameters 


Station 

Province 

Water 

(m) 

uncor. 

depth 

cor . 

Bottom  water 
sound  speed 
(m/ sec) 

Sediment 

thickness 

(sec) 

Depth  of  LF 
reflectors 
(sec) 

No.  of  HF 
reflectors 
per  meter 

Density 

(gm/cm3) 

CR£ 

0.977 

8  i 

(sec  i) 

dB/m^kHz 

2 

Northern 

Iberian 

Abyssal 

Plain 

5250- 

5250 

5306- 

5306 

1551.0 
@  5291  m 

1.2 

1.2 

0.11 

0.38 

0.3-03 

1.51 

0.7 

0.078 

4 

Southeas  t 
Tagus 
Abyssal 
Plain 

5075- 

5065 

5131- 

5121 

1548.5 
@  5115  m 

2.0 

0.09 

0.35 

1.09 

0.45-0.6 

1.53 

0.994 

0.7 

0.076 

5 

Southwest 

Tagus 

Abyssal 

Plain 

5066- 

5065 

5121- 

5120 

1548.6 
@  5120  m 

2.0 

0.05 

0.13 

0.30 

1.00 

0.35-0.6 

1.54 

0.996 

0.7 

0.07  7 

7 

Dragon 

Seamount 

1950- 

1210- 

2700 

1957- 

1218- 

2711 

1503.6 
§  2264  m 

0.001 

2.5- 

2.9 

2.5- 

3.4 

8 

Dragon 

Seamount 

Plateau 

2410- 

2085 

2420- 

2094 

1504.8 
@  2480  m 

0.9 

0.14 

0.51 

0.85 

0.2-0. 3 

1.60 

1.01/3 

2.3 

0.250 

9 

Lower 

Co  ntinen- 

4685- 

4654 

4724- 

4693 

1525.5 
@  4000  m 

1.3 

0.05 

0.11 

0.5-0. 3 

1.52 

0.984 

1.1 

0.100 

tal  Rise  0.23 

South  of 

Madeira 


1500  m/sec.  The  choice  of  the  reflectors 
is  somewhat  subjective,  but  appears  to 
correlate  with  the  layer  interfaces 
determined  from  the  acoustic  data. 

Two  values  for  the  number  of  discrete 
reflectors  per  meter  as  determined  from 
the  3.5  kHz  (HF)  depth  recorder  are 
given.  These  values  are  where  the  first 
and  last  charges  were  dropped,  respect¬ 
ively.  The  maximum  depth  of  penetration 
observed  in  the  records  is  about  50  m 
(assuming  a  sound  speed  of  1500  m/sec  in 
the  sediment).  The  listed  values  are 
determined  from  approximately  the  upper 
25  m.  This  parameter  is  further  discus¬ 
sed  in  Appendix  D.2. 

The  final  four  columns  list  the  geo¬ 
acoustic  parameters  of  the  sediments. 
The  justification  for  choosing  the  par¬ 
ticular  values  is  given  in  Appendix  C. 

A  comparison  of  these  parameters  as 
derived  from  the  geological  and  acousti¬ 
cal  data  will  be  discussed  in  Section  VI 
of  this  report. 

The  wet  density,  P  ,  is  measured  for  core 
samples  or  is  estimated  from  the  known 
geological  properties  of  the  sediments. 

The  ratio  of  the  sound-speed  of  the 
sediment  to  that  of  the  bottom  water 
(CR)  was  generally  estimated  from  core 
data.  The  average  sound  speed  for  the 
upper  5  m  of  the  sea  floor  is  used  for 
this  value.  The  sediment  sound-speed  is, 
in  general,  lower  than  that  of  the  bot¬ 
tom  water.  Exceptions  are  the  sediments 
with  high  carbonate  content  at  Station  8 
and  the  basalt  outcrop  with  a  negligible 
sedimentary  cover  at  Station  7. 

The  sound-speed  gradient,  g,  is  defined 
in  this  report  as  the  difference  in 
sound-speed  at  200  m  and  at  the  sediment 
surface  divided  by  200  m.  Various  meth¬ 
ods  used  to  derive  this  value  are  dis¬ 
cussed  in  detail  in  the  appendices. 
Sonobuoy  data  are  used  in  the  Station  9 
determination.  At  the  other  stations, 
the  two-way  travel  time  to  a  reflector 
of  known  depth  from  a  DSDP  drill  hole 
was  used  to  determine  the  gradient. 


The  last  column  is  the  attenuation  con¬ 
stant  k  •  This  parameter  is  the  propor¬ 
tionality  constant  in  the  equation  that 
assumes  a  linear  dependence  on  attenua¬ 
tion  frequency: a  =  Kf.  The  values  listed 
are  based  on  core  samples  and  relation¬ 
ships  presented  by  Hamilton  [45]. 

VI.  Comparison  of  Acoustical  and  Geological  Results 

Values  of  the  sediment-to- water  density 
ratio  (DR)  and  sound-speed  ratio  (CR)  as 
estimated  from  the  acoustical  and  geo¬ 
logical  results  are  given  in  Table  4. 
Absolute  differences  between  the  acous¬ 
tically  inferred  sediment-to-water  den¬ 
sity  ratio  using  the  Rayleigh-equation 
method  (DRR)  and  the  core-measured 
values  averaged  over  the  top  5  m  of  sed¬ 
iment  (DR q)  averaged  8%.  The  density 
ratio  is  probably  the  least  sensitive 
acoustic  parameter,  and  a  value  of  1.5 
for  all  of  the  stations  is  probably  as 
good  a  choice  as  either  the  DRR  or 
DRq  values. 

Similarly,  the  sediment-to-water  sound- 
speed  ratio  determined  by  application  of 
the  Rayleigh-ref lection  equation  (CRR) 
and  time-of-arrival  methods  (CRT)  may 
be  compared  with  the  geological  value 
(CRq).  Absolute  differences  between 
CRr  and  CRq  averaged  2%,  whereas 
absolute  differences  between  CR*p  and 
CRq  averaged  4%.  The  larger  differ¬ 
ences  for  CR*p  may  be  expected,  since 
CR^p  is  obtained  from  extrapolations  of 
the  behavior  at  depth.  CRR,  on  the 
other  hand,  is  sensitive  to  just  the 
reflection  at  the  water-sediment  inter¬ 
face. 

The  differences  between  CRR  and  CRp 
were  largely  reconciled  at  Stations  2 
and  8  by  the  inclusion  of  a  thin  layer 
at  the  top  of  the  sediments.  The  extra 
layer  provided  a  CRp-  value  more  nearly 
equal  to  CRR,  and,  at  the  same  time, 
gave  correct  arrival-times  for  both  the 
high-frequency  and  low-frequency  impulse 
responses.  Presumably  the  differences  at 
the  other  stations  could  also  be  recon¬ 
ciled  by  more  complicated  layered 
models . 
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Table  4.  Comparison  of  acoustical  and  geological  results 


Acoustic 


Station 

Model 

DRr 

DRg 

CRr 

CRt 

2 

3 

1.4 

1.5 

0.97 

1.05 

2 

5 

1.4 

1.5 

0.97 

1.00 

2 

6 

1.4 

1.5 

0.97 

0.96 

4 

1 

1.3 

1.5 

1.04 

1.06 

4 

4 

1.3 

1.5 

1.04 

1.04 

4 

3 

1.3 

1.5 

1.04 

1.00 

5 

3 

1.6 

1.5 

0.99 

1.03 

5 

6 

1.6 

1.5 

0.99 

0.99 

8 

1 

1.4 

1.6 

0.99 

1.08 

8 

2 

1.4 

1.6 

0.99 

1.03 

8 

4 

1.4 

1.6 

0.99 

1.01 

9 

1 

1.5 

1.5 

0.97 

1.06 

9 

7 

1.5 

1.5 

0.97 

0.97 

CR g 

fcA 

(msec) 

fcG 

(msec) 

KA  <G 

dB/m/kHz 

0.98 

6.1 

103 

0.05 

0.08 

0.98 

8.4 

103 

0.05 

0.08 

0.98 

8.6 

103 

0.06 

0.08 

0.99 

0.5 

33 

0.04 

0.08 

0.99 

0.9 

33 

0.06 

0.08 

0.99 

0.5 

33 

0.07 

0.08 

1.00 

4.8 

50 

0.02 

0.08 

1.00 

4.0 

50 

0.03 

0.08 

1.01 

1.9 

42 

0.51 

0.25 

1.01 

1.0 

42 

0.51 

0.25 

1.01 

0.9 

42 

0.27 

0.25 

0.98 

1.2 

64 

0.27 

0.10 

0.98 

1.4 

64 

0.30 

0.10 

For  a  single-layer  model  at  low  fre¬ 
quencies,  CR*p  should  be  used  rather 
than  CRr.  The  use  of  CR-j.  generates 
small  errors  in  the  amplitude  of  the 
first  reflected  pulse,  but  gives  the 
correct  arrival  times.  Alternatively, 
the  use  of  CR^  gives  the  correct 
amplitude  for  the  first  reflected  pulse, 
but  causes  significant  errors  in  the 
arrival-times . 

The  best  estimate  of  sediment  sound- 
speed  structure  based  solely  on  geologi¬ 
cal  data  at  each  station  is  described  as 
a  single  layer  with  a  gradient  as  defin¬ 
ed  in  Section  V.  The  gradients  are  list¬ 
ed  in  Table  3,  The  corresponding  sound- 
speed  profiles  are  indicated  by  arrows 
in  Figure  10.  The  acoustic  models  are 
more  complex;  the  corresponding  param¬ 
eters  are  listed  in  Table  4  and  pro¬ 
files  illustrated  in  Figure  10.  An  ap¬ 
propriate  means  to  compare  geological 
and  acoustical  models  is  by  the  arrival 
times.  Arrival-time  differences  between 
the  predicted  and  actually  observed 
times  were  calculated  and  the  RMS 


differences  t^  and  tG  presented  in 
Table  4. 

For  each  of  the  stations,  the  acous¬ 
tical  RMS  arrival-time  error  t^  was 
less  than  9  msec.  In  consideration  of 
the  pulse  width  (14.2  msec  correspond¬ 
ing  to  full  width  at  3  dB  points),  the 
theoretical-bottom-loss  pattern  would  be 
expected  to  match  the  observed  pattern 
to  within  half  an  oscillation  period. 
The  RMS  time  difference  was  largest  for 
Stations  2  and  5.  This  might  be  expect¬ 
ed,  since  the  model  at  these  stations 
was  forced  to  describe  a  much  thicker 
sediment  layer.  The  geologically  pre¬ 
dicted  arrival  times  yielded  errors  t  q 
that  were  much  greater  than  the  pulse 
width;  therefore,  there  would  be  no  hope 
in  reproducing  the  bottom-loss  oscilla¬ 
tion  pattern. 

The  uniqueness  of  modeling  tends  to  be 
more  of  a  problem  at  shallow  depths.  For 
example,  at  Stations  2  and  8  where  addi¬ 
tional  modeling  of  the  high-frequency 
impulse  responses  was  attempted,  the 
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time-of -arrival  behavior  did  not  clearly 
separate  into  constant-sound-speed  or 
sound-speed-gradient  behavior  at  shallow 
depths  ( <1 0  m).  This  difficulty,  how¬ 
ever,  is  not  serious  (as  it  would  be  at 
greater  depths)  because  the  absolute 
differences  in  sound  speed  for  the  dif¬ 
ferent  candidate  models  are  small  over  a 
thin  layer.  At  greater  depths,  the 
models  are  clearly  differentiated  by  the 
arrival-time  behavior  as  demonstrated  at 
Stations  2  and  5. 

Details  of  the  shallow  sediment  struc¬ 
ture  (<10  m)  are  relatively  unimportant 
at  low  frequencies  (wavelength  >>10  m) . 
It  primarily  contributes  to  the  effect¬ 
ive  low-frequency  Rayleigh-type  reflec¬ 
tion  at  the  water-sediment  interface.  At 
higher  frequencies  (wavelength  <<10  m) , 
shallow  penetrating  rays  are  well  resol¬ 
ved  from  the  water-sediment  reflected 
pulse,  and  the  details  of  the  sound- 
speed  structure  are  important. 

The  shallow  sediment  structure  inferred 
at  Stations  2  and  8  (Fig.  10)  indicated 
a  rapid  initial  increase  in  sound-speed 
at  the  top  of  the  sediments.  As  discus¬ 
sed  previously,  more  than  one  model  gave 
a  reasonable  time-of-arrival  fit.  The 
shallow  structure  could  be  modeled  by 
either  a  constant  or  a  gradient  sound- 
speed  layer.  The  gradient  layer  permit¬ 
ted  a  lower  value  of  CR^  and  gave 
approximately  the  same  average  sound 
speed  as  the  constant-sound-speed  model. 
Although  the  exact  form  of  the  sound- 
speed  profile  could  not  be  unambiguously 
decided,  the  models  were  consistent  with 
a  sound-speed  which  rose  rapidly  to  val¬ 
ues  on  the  order  of  1600-1700  m/ sec 
within  the  first  10  m  for  Station  2,  and 
within  the  first  50  m  for  Station  8. 

From  the  geological  point  of  view,  the 
shallow  sound-speeds  are  higher  than 
would  be  expected  for  the  sediments 
involved.  In  addition,  the  associated 
sound-speed  gradients  are  much  higher 
than  any  previously  measured.  On  this 
basis,  one  might  drop  the  gradient-layer 
interpretation  and  choose  the  constant- 
sound-speed-layer  model.  In  that  case, 
however,  the  sound-speed  just  takes 


discontinuous  jumps  and  still  arrives  at 
the  high  sound-speeds. 

Although  a  variety  of  bottom  sediment 
structures  could  account  for  the  obser¬ 
ved  high  sound-speeds,  their  occurrence 
is  not  readily  predicted  on  the  basis  of 
the  available  geological  data  alone. 
Nevertheless ,  geological  observations 
are  often  suggestive  or  supportive  of 
acoustic  findings.  For  example,  a  very 
hard  clay  interface  was  sampled  in  the 
core  cutter  from  a  core  taken  south  of 
Station  2.  Likewise,  near  Station  8, 
there  was  an  indication  of  a  highly 
reflective  buried  layer  at  the  start  of 
the  seismic  profile  track.  This  evidence 
of  high-sound-speed  layers  from  the  core 
samples  may  indicate  the  cause  of  the 
acoustically  inferred  high-sound-speeds 
at  shallow  depths  observed  at  Stations  2 
and  8. 

Comparison  of  the  shallow  sediment 
acoustical  results  with  geological  the¬ 
ory  is  not  direct.  Deep-water  sound- 
speed  gradients  determined  from  sono- 
buoy  interval-velocity  measurements  are 
only  available  for  sediment  depths  more 
than  a  few  hundred  meters.  It  would 
hardly  seem  reliable  to  extrapolate 
these  values  for  use  at  0-50  m.  Gravity- 
and  piston-type  cores  cannot  provide  the 
desired  information  because  they  typi¬ 
cally  penetrate  less  than  5  m  of  sedi¬ 
ment  . 

Effective  shallow  sound-speed  gradients 
need  not  be  explained  entirely  in  terms 
of  sediment  compaction,  temperature 
gradients,  or  pressure  effects.  A  dis¬ 
tribution  of  high-sound-speed  sediment 
layers  among  other  low-speed  layers 
could  produce  apparent  high  gradients  at 
shallow  depths,  particularly  when  there 
are  geological  reasons  to  expect  the 
most  recent  sediment  deposition  to  be  of 
the  slower  sound-speed  type.  At  greater 
depths  where  the  comparison  was  more 
direct  (e.g.,  Stations  2  and  5),  the 
acoustically  inferred  sound-speed  grad¬ 
ients  fell  well  within  the  range  of  pre¬ 
viously  reported  values  and  were  com¬ 
pletely  consistent  with  geological 
expectations  for  the  particular  sites. 
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The  geologically  inferred  sound-speed 
gradient  using  an  interval-velocity 
measurement ,  deep-drill  core-measured 
depth,  and  core-measured  CR  value  (CRIC 
method)  was  in  good  agreement  with  the 
acoustically  determined  gradient  for 
acoustic  ray  paths  which  penetrated  200- 
300  m  in  depth  (Stations  2  and  5).  This 
suggests  that  extrapolation  of  deep 
sediment  gradients  can  be  made  from  re¬ 
mote  sites  having  similar  sedimentolog- 
ical  histories.  Although  gradients  were 
accurately  predicted,  the  predicted  ar¬ 
rival  times  t q  were  in  error  by  much 
more  than  a  pulse  width  (see  Table  4). 
In  other  cases  (Stations  8  and  9),  the 
geological  prediction  of  a  deep  gradient 
was  not  applicable  because  either  sig¬ 
nificant  reflection  or  total  reflection 
was  occurring  at  a  much  shallower  layer. 

Estimates  of  acoustic  attenuation  are 
compared  in  Table  4  where  is  the 

acoustical  estimate  and  kq  is  the 
geological  estimate.  For  each  set  of 
values,  Stations  2,  4,  and  5  have  the 
lowest  values  and  Station  8  has  the 
highest,  although  the  acoustic  results 
have  a  greater  range  of  values.  The 
attenuation  parameter  is  the  most  diffi¬ 
cult  parameter  to  estimate  in  either 
approach.  The  geological  estimate  is 
based  on  data  for  which  there  is  much 
scatter  and  the  acoustical  estimate  is 
based  on  so  few  points  it  is  difficult 
to  estimate  its  probable  error. 

Many  sources  of  error  may  contribute  to 
discrepancies  between  acoustically  and 
geologically  inferred  model  parameters. 
The  acoustical  analysis  is  limited  by 
the  imposition  of  a  priori  model  types 
to  which  the  data  is  fitted.  Simple 
models  were  chosen  to  describe  the  beha¬ 
vior  of  dominant  features  in  the  receiv¬ 
ed  signals.  Many  details  of  the  signals 
are  not  included  and,  of  course,  even 
the  dominant  features  are  not  modeled 
exactly. 

Some  sources  of  error  in  the  geological 
analysis  may  stem  from  the  following 
considerations  of  the  way  that  geologi¬ 
cal  survey  data  is  collected:  (1)  On 


seismic-profile  surveys  clipped  signals 
are  recorded  in  order  to  reveal  deep 
reflecting  horizons.  The  accompanying 
disadvantage  is  that  differences  in 
layer  reflectivities  may  not  be  differ¬ 
entiated.  (2)  Surface  reflections  and 
bubble  pulses  in  profiler  and  sonobuoy 
data  limit  depth  resolution  and  preclude 
accurate  measurement  at  depths  of  0  to 
200  m — the  depth  range  that  is  most 
influential  in  low-frequency,  low-graz¬ 
ing-angle  ,  bottom-sediment  interaction. 
(3)  Surficial  sound  speeds  measured  in 
cores  or  extrapolated  from  deeper  meas¬ 
urements  may  not  be  the  same  as  the 
effective  low-frequency,  interface- 
sound-speed  parameter  required  by  the 
models.  (4)  Sound-speed  gradients  pre¬ 
dicted  from  interval-velocity  analyses 
of  layers  deeper  than  200  m  may  not  be 
applicable  at  shallower  depths.  (5) 
Sediment  sampling  in  DSDP  cores  is  often 
insufficient  at  depths  of  0  to  200  m. 

VII.  Conclusions 

A  deep  water  experiment  employing  a  deep 
hydrophone  and  0.8  kg  SUS  charges  deto¬ 
nated  at  244  m  depth  provides  data  use¬ 
ful  for  modeling  low-frequency,  low- 
grazing-angle,  bottom-interaction  phe¬ 
nomena.  Deconvolution  processing  of  bot¬ 
tom  reflectivity  data  eliminates  multi¬ 
ple  bubble  pulses  which  otherwise  inter¬ 
fere  with  the  data  analysis  and  inter¬ 
pretation. 

We  draw  the  following  conclusions  about 
the  specific  sites  analyzed  in  this 
report: 

(1)  Impulse  responses  from  thick-sedi¬ 
ment  bottom  sites  displayed  sequences  of 
pulses  in  an  easily  discerned  pattern 
varying  with  grazing-angle.  Decay  times 
for  the  overall  pattern  were  on  the 
order  of  0.1  sec.  At  low  frequencies  (40 
to  90  Hz)  and  over  some  restricted  range 
of  grazing-angles  bottom  impulse  respon¬ 
ses  could  be  interpreted  in  terms  of 
simple  layered  sediment  models  consist¬ 
ing  of  one  to  three  refracting  and/or 
reflecting  layers.  The  RMS  error  in  the 
model-predicted  arrival  time  along  the 
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dominant  bottom  interacting  path  t^ 
ranged  from  1  to  9  msec — less  than  the 
14.2  msec  acoustic  pulse  width. 

(2)  In  marked  contrast  the  impulse  res¬ 
ponses  from  a  sediment-free  bottom  site 
showed  no  distinctive  pattern  (sugges¬ 
tive  of  roughness  scattering)  and  decay¬ 
ed  in  times  on  the  order  of  1  sec.  No 
modeling  was  attempted. 

(3)  Geological  and  acoustical  estimates 
of  the  sediment-to-water  density  ratio 
(DRr,  DRq)  and  the  sound-speed 
ratios  (CRR,  CRq)  agreed  to  within 
5%. 

(4)  Differences  between  the  geologi¬ 
cally  and  acoustically  inferred  sound- 
speed  profiles  were  on  the  order  of  100 
msec  and  varied  considerably  with  model 
and  depth.  The  geologically  inferred 
sound  speed  was  always  lower,  with  an 
RMS  error  in  predicted  arrival  time 
along  the  dominant  bottom-interacting 
path  t q  of  33  to  103  msec,  i.e.,  from 
1  to  5  times  the  acoustic  pulse  width. 

(5)  In  all  cases  relatively  low  bottom- 
losses  were  measured  over  angular  ranges 
for  which  a  large,  dominant  peak  was 
observed  in  the  impulse  response.  The 
dominant  peak  was  attributed  to  a  total¬ 
ly  refracting  or  reflecting  sediment 
path  that  penetrated  the  sediments  to 
less  than  200  m  depth.  The  only  loss 
mechanism  was  hypothesized  to  be  an  ef¬ 
fective  attenuation  along  the  dominant 
path. 

(6)  Acoustically  inferred,  effective 
attenuation  coefficients,  estimated  to 
be  uncertain  by  at  least  a  factor  of  2, 
ranged  from  0.02  to  0.05  dB/m/kHz. 
Although  the  geological  and  acoustical 
values  showed  the  same  trend  from 
station  to  station,  ratios  of  the  two 
inferred  values  ranged  from  1  to  4. 

Because  of  the  geographical  and  geo¬ 
logical  diversity  of  the  thick-sediment 
sites  analyzed  in  this  study,  it  is 
expected  that  the  same  types  of  low- 
frequency  bottom-interaction  behavior 
would  be  observed  at  many  other  thick 


sediment  sites.  When  this  is  the  case, 
simple  geoacoustic  models  are  possible 
for  predicting  those  aspects  of  surveil¬ 
lance-system  performance  that  depend  on 
bottom  reflectivity.  In  the  absence  of 
detailed  wide-angle  reflectivity  meas¬ 
urements,  such  as  those  analyzed  in  thi/ 
study,  the  derivation  of  accurate  geo¬ 
acoustic  models  on  the  basis  of  general¬ 
ly  available  geological  data  is  limited. 
Nevertheless,  the  utility  of  such  data 
could  be  significantly  enhanced,  as 
suggested  in  the  recommendations  that 
follow. 

VIII.  Recommendations 

The  following  recommendations  are  sug¬ 
gested  for  improving  the  utility  of 
normal-incidence  subbottom-profiler  and 
deep-drill-core  data  for  the  predic¬ 
tion  of  accurate  geoacoustic  models: 

1.  Record  undistorted  signal  amplitude 
of  normal-incidence  profiler  returns, 

1. e.,  without  clipping. 

2.  Apply  signal  processing  techniques  of 
deconvolution  and  deghosting  to  remove 
interference  due  to  bubble  pulses  and 
water-surface  reflections . 

3.  Measure  sound-speed  and  density  in 
the  upper  few  hundred  meters  of  a  deep- 
drill  core  (DSDP)  at  sample  spacings 
small  relative  to  a  wavelength. 

4.  Develop  mathematical  inversion  tech¬ 
niques  of  sound-speed-profile  recon¬ 
struction  that  are  sensitive  at  the  rel¬ 
evant  depths  of  penetration  (0-200  m) 
and  to  the  averaging  effects  of  low  fre¬ 
quencies  and  low  grazing  angles. 
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X.  Illustrations 
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Figure  1  .  Physiographic  map  of  the  Northeast  Atlantic  Ocean  with  locations 
of  bomb  runs  for  stations. 
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Figure  2.  Typical  geometry  of  bomb  run  for  reflectivity  experiment. 


(I)  (2) 


C(Z)  =  Cw  Z<0 
C(Z)  =  C,  0<  Z  <  h 
CR  =  C,  /  Cw 

Figure  3.  Model  1:  Constant-sound-speed  layer  overlying  half-space  with 
reflection  coefficient  R  =  i. 
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(I)  (2) 


(3) 


C  (  Z ) =  C w  Z<0 
C  (Z)  =  C,  0<  Z  <h, 

C(Z)=C2  h,  <  Z  <  h,  +  h2 
CR,=  C,  /  Cw  CR2=  C2/  C, 

Figure  4.  Model  2:  Two  constant-sound-speed  layers  overlying  half-space  with 
reflection  coefficient  R  =  i. 
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Figure  5.  Model  3:  Pseudo-linear-sound-speed  half-space. 
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CR  =  C,  /Cw  c(z)  =  c w  Z<0 

C(Z)=  C,  ( I-  2g  Z/C  , 1/2  0<  Z  <  h 

Figure  6.  Model  4:  Pseudo-linear-sound-speed  layer  overlying  half-space  with 
reflection  coefficient  R  =  i. 
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Figure  7.  Model  5:  Constant-sound-speed  layer  overlying  a  pseudo-linear- 
sound-speed  half-space. 
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(I)  (2)  (3) 


C(Z)  =  Cw  Z<0 

C(Z)=  C,  ( I  -  2g,  Z  /  C ,  ) " 1/2  0<Z<  h, 

C(Z)  =  C2(  I  - 2g2(  Z - h , )  /  C2)  _  1/2  Z  >  h, 

CR,  =  C,  /  Cw 

Figure  8.  Model  6:  Pseudo-linear-sound-speed  layer  overlying  pseudo-linear- 
sound-speed  half-space. 
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Figure  9.  Model  7:  Pseudo-linear-sound-speed  layer  overlying  a  constant-sound- 
speed  layer,  overlying  a  half-space  with  reflection  coefficient  R  =  i. 
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Figure  10.  Predicted  sound-speed  profiles  in  sediments  at  each  station.  Solid  lines 
represent  models  derived  from  geological/geophysical  analyses  and  are 
discussed  in  Appendices  C  and  E  where:  C  =  Cernock  (50);  H  =  Hamilton  et  al . 
(51);  HE  =  Houtz  and  Ewing  (39);  HEB  =  Houtz,  Ewing,  and  Buhl  (60);  L  = 
modified  after  Laughton  (40);  S  =  based  on  surficial  sediment  data  from 
Hamilton  (82);  SB  =  sonobuoy  data  from  this  report;  118  =  based  on  DSDP  Site 
118  data  (49);  and  136  =  based  on  DSDP  Site  136  data  (59).  The  arrow  indicates 
the  choice  of  profile  based  on  geological  considerations.  Dashed  lines  represent 
models  derived  from  acoustical  analyses,  with  numbers  corresponding  to  models 
in  Figures  3-9.  Hatch  marks  indicate  an  impenetrable  interface. 
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Figure  1 1.  Station  2  impulse  responses. 
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Figure  12.  Station  2  impulse  responses  on  expanded  time  axis. 
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Figure  13.  Station  2  one-third-octave  bottom-loss  curves. 
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Figure  14.  Station  2  one-third-octave  bottom-loss  curves. 
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Figure  15.  Station  2  Rayleigh-equation  analysis. 
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Figure  16.  Station  2  model  3  analysis. 
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Figure  18.  Station  2  models  5  and  6  analysis  (deeper  layer). 
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Figure  19.  Station  4  impulse  responses. 
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Figure  20.  Station  4  impulse  responses  on  expanded  time  axis. 
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Figure  21 .  Station  4  one-third-octave  bottom-loss  curves. 
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Figure  22.  Station  4  one-third-octave  bottom-loss  curves. 
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STATION  4 


Figure  25.  Station  4  model  4  analysis  . 
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Figure  26.  Station  5  impulse  responses. 
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Figure  27.  Station  5  impulse  responses  on  expanded  time  axis. 
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Figure  28.  Station  5  one-third-octave  bottom-loss  curves. 
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Figure  29.  Station  5  one-third- octave  bottom-loss  curves. 
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STATION  5 


Figure  30.  Station  5  Rayleigh-equation  analysis. 
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Figure  31  .  Station  5  model  3  analysis. 
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Figure  32.  Station  5  model  6  analysis. 
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Figure  33.  StaHon  8  impulse  responses. 
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Figure  34.  Station  8  impulse  responses  on  expanded  time  axis. 
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Figure  35.  Station  8  one-third-octave  bottom-loss  curve. 
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Figure  36.  Station  8  one-third-octave  bottom-loss  curves. 
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STATION  8 


Figure  37.  Station  8  Rayleigh-equation  analysis. 
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Figure  38.  Station  8  model  1  analysis. 
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Figure  40.  Station  8  model  4  analysis. 
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Figure  41.  Station  9  impulse  responses. 
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Figure  42.  Station  9  impulse  responses  on  expanded  time  axis. 
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Figure  43.  Station  9  one-third-octave  bottom-loss  curves. 
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Figure  44.  Station  9  one-third-octave  bottom-loss  curves. 
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STATION  9 


Figure  45.  Station  9  Rayleigh-equation  analysis. 
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Figure  46.  Station  9  model  1  analysis. 
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Figure  47.  Station  9  model  7  analysis. 
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Figure  48.  Station  7  impulse  responses. 
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Figure  49.  Station  7  impulse  responses  on  expanded  time  axis. 
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Figure  50.  Station  7  one-third-octave  bottom-loss  curves. 
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Figure  52.  Area  map  of  station  2  in  the  Iberian  Abyssal  Plain.  Sediment  thickness 
is  indicated  in  seconds  of  two-way  travel  time. 
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Figure  53.  Seismic-profile  records  near  station  2 
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Figure  54.  Predicted  sediment  sound-speed  profiles  at  station  2. 
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Figure  55.  Area  map  of  stations  4  and  5  in  the  Tagus  Abyssal  Plain. 


— I 
U5 


0000 


TIME/HOURS  (ZULU) 
2330 


STATION  4  SHOT-# I 
2300  \ 


0100 


0030 


0200 


0130 


Figure  56.  3.5  kHz  subbottom  depth-recorder  records  along  bomb  run  at  station  4 
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Figure  57.  3.5  kHz  subbottom  depth-recorder  records  along  bomb  run  at  station  5 
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Figure  58.  Seismic-profiler  records  near  stations  4  and  5 
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Figure  59.  Sound-speed  in  cores  from  the  Tagus  Abyssal  Plain. 
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Figure  60.  Area  map  of  stations  7  and  8  at  Dragon  Seamount. 
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Figure  61  .  3.5  kHz  subbottom  depth-recorder  record  along  bomb  run  at  station  8. 
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Figure  62.  Seismic-profiler  records  near  station  8 
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Figure  63.  3.5  kHz  and  seismic-profiler  records  at  stations  7  and  8. 
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Figure  64.  Sound-speed  in  cores  taken  at  stations  7  and  8. 
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Figure  65.  Predicted  sediment  sound-speed  profiles  at  station  8. 
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Figure  66.  Area  map  for  station  9  on  the  Lower  Continental  Rise  south  of  Madeira. 
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Figure  67.  3.5  kHz  and  seismic-profiler  records  along  bomb  run  at  station  9 
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Figure  68.  Sound-speed  gradient  at  DSDP  Site  1 18  determined  with  CRIC  method. 
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Figure  69.  Results  of  sonobuoy  stations  off  Northwest  Africa. 
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Figure  72.  Sediment  porosity  values  at  DSDP  Site  136. 
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Figure  73.  Water  sound-speed  profiles  at  acoustic  stations. 


Appendix  A. 

Tables  of  Formulas  Used  in  Acoustical  Analysis 


Formulas  given  in  Table  5  specify  condi¬ 
tions  to  be  satisfied  for  the  validity 
of  plane-wave  approximations.  The  condi¬ 
tions  are  expressed  in  terms  of  computed 
parameters  ek,  er^,  and  that  must 
be  small  compared  to  one.  The  formulas 
were  obtained  from  a  calculation  of  re¬ 
flectivity  from  a  point  source  using  the 
method  of  steepest  descents  following 
Brekhovskikh  [38].  The  formula  for  e^. 
was  obtained  simply  from  the  geometry  of 
eigenrays  connecting  source  and  receiv¬ 
er. 

The  reflectivity  calculation  was  applied 
only  to  singly  reflected  or  refracted 
pulse  arrivals  in  the  context  of  a  ray 
theoretical  decomposition,  and  not  to 
the  total  response  from  a  continuous 
wave  (CW)  source.  The  condition  on  k  was 
imposed  in  each  case  to  assume  that 
pulse  arrivals  were  separated  in  time 
for  a  pulse  bandwidth  Am  =  kc.  Otherwise, 
a  separate  restriction  on  k  was  not  re¬ 
quired  for  the  plane  wave  approxmatlon. 
Restrictions  were  found  for  products  of 
the  form  krn  where  n  -  1,  2,  or  3,  but 
these  were  automatically  satisfied  by 
the  separate  conditions  on  k  and  R.  In 
the  case  of  the  pseudolinear  half space, 
the  restriction  on  k  also  satisfies 
requirements  for  the  geometrical  optics 
approximation  (ray  theory)  for  an  obser¬ 
ver  at  the  top  of  the  gradient  half¬ 
space  (Brekhovskikh  [38],  pg.  206,  eq. 
16.39).  When  conditions  for  the  plane- 
wave  approximation  are  fulfilled,  all 
pulse  arrivals  may  be  considered  to 
arise  from  parallel  rays  at  the  same 
angle  of  incidence. 

Table  6  gives  numerical  values  for  the 
quantities  involved  in  the  inequalities 
of  Table  5.  The  ratios  e/R  and  e/k  rep¬ 
resent  the  order  of  magnitude  of  the 
fractional  error  induced  by  neglect  of 
non— plane— wave  terms.  These  ratios  are 


for  worst  case  conditions  at  the  extreme 
of  the  indicated  angular  interval.  The 
error  tends  to  drop  off  rapidly  as  the 
angle  moves  away  from  the  end  points  of 
the  interval;  thus,  much  of  the  data 
corresponds  to  a  better  plane— wave 
approximation  than  indicated  for  the 
worst  case  conditions  of  Table  6. 

The  ratio  £ t/R  is  less  than  0.1  in  all 
cases  and  indicates  a  good  plane— wave 
arrival  time  approximation.  The  other 
three  ratios  (ek>  erl>  erh)  address 
all  aspects  of  the  plane  approximation, 
including  interpretation  of  pulse 
shapes,  amplitudes,  and  arrival  times. 
The  conditions  for  which  ek  erl>  anc* 
erh  are  defined  are  sufficient  for  the 
validity  of  the  plane-wave  approxima¬ 
tion.  The  condition  on  et  is  thus,  in 
some  sense,  a  subset  of  these  conditions 
and  is  far  less  restrictive.  Correspond¬ 
ingly,  the  error  term  (  e  t)  is  much 
less  than  the  more  general  error  terms 
(Ek,  erl,  £rh)* 

The  distinction  between  et  and  e  k, 
£rl»  and  e  ri,  is  important  because 
the  sound— speed  parameters  were  extract¬ 
ed  from  the  data  on  the  basis  of  arrival 
time  analysis.  The  ratio  e  t/R  repre¬ 
sents  the  plane  wave  error  in  arrival 
times;  its  smaller  values  suggest  that 
good  arrival-time  accuracy  is  obtained 
in  the  plane-wave  approximation,  even 
though  the  total  error  is  larger  than 
desired  in  some  cases.  The  larger  values 
for  ek,  er^,  and  erk  may  be  expected 
to  signify  greater  error  in  determina¬ 
tions  based  on  peak  amplitudes  such  as 
in  the  estimation  of  attenuation  coeffi¬ 
cients  . 

Table  7  lists  generalized  formulas  for 
ray  path-lengths  and  travel  times  for 
layers  with  or  without  a  turning  point 
and  for  zero-gradient  layers.  Table  8 
gives  formulas  for  extracting  sound- 
speed  and  layer  thickness  parameters  by 
fitting  straight  lines  to  measured  ar¬ 
rival  time  plots.  Formulas  for  computing 
average  attenuation  are  given  in  Table 
9. 
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Table  5.  Conditions  to  be  satisfied  for  plane  wave  approximation 


The  symbols  used  in  the  table  have  the  following  meanings:  r  is 
the  distance  from  the  source-point  image  to  the  receiver;  9, 
0C,  9r,  are  the  grazing,  critical  grazing,  and  refracted 
grazing  angles,  respectively;  k  is  the  wave  number  (2-Tr/A);  C  is 
the  sound  speed  in  the  water  column,  g  is  the  sound  speed  gradi¬ 
ent  in  the  bottom;  and  h  is  the  layer  thickness.  The  column 
marked  "grazing-angles"  indicates  that  the  inequality  must  be 
satisfied  at  either  the  highest  or  lowest  grazing-angles  used. 
The  first  three  conditions  for  each  model  are  sufficient  to 
establish  validity  of  the  plane-wave  approximation.  The  fourth 
inequality  involving  et  is  sufficient  for  the  accuracy  of 
plane-wave  arrival  times  and  is  less  restrictive  then  the  other 
three  conditions. 


Model 


Grazing- 

Requirements  for  plane  wave  approximation_ angles 


Constant  sound 

k 

» 

ek 

ek  =  1/ (2h  Sin  9r) 

low 

speed  layer 
with  totally 

r 

» 

erl 

Erl  =  2h  Sin20/(Sin3  0r) 

low 

reflecting  bottom 
interface 

r 

» 

e  u 
rh 

£rh  =  h^4 

high 

r 

» 

et 

et  =  2h  Sin  9  Cos  0r/Sin  0r 

low 

Pseudo-linear 

k 

» 

ek 

ek  =  3g/(2C  Sin3  er) 

low 

sound  speed 
half  space 

r 

» 

erl 

erl  =  2c  Sin2  ©/ ( g  Sin  0r) 

low 

r 

» 

erh 

erh  =  4c/g 

high 

r 

» 

£t 

et  =  (C/g)Sin(20r)Sin0 

high 
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Table  6.  Numerical  values  of  the  inequality  conditions  required 
for  plane  wave  approximation 


All  quantities  are  defined  in  Table  5;  angles  are  in  degrees. 
The  columns  labeled  0  give  the  range  of  grazing-angles  over 
which  the  fractional  errors  (e/k,e/R)  may  be  expected  to  be  less 
than  the  values  listed  in  the  adjacent  columns.  The  fractional 
errors  give  an  order  of  magnitude  measure  of  second-order  terms 
in  a  power  series  expansion  of  the  reflectivity  in  which  the 
first-order  term  is  the  plane-wave  contribution. 


Sta. 

Model 

eC 

ek/k 

erl  /R 

e 

£rh/R 

9 

VR 

0 

2 

3 

17.8 

0.2 

0.1 

>22 

0.5 

<35 

0.08 

<35 

4 

1 

19.4 

0.1 

0.1 

>23 

0.002 

<90 

0.02 

>23 

4 

3 

0 

0.4 

0.02 

>15 

0.3 

<35 

0.02 

<35 

5 

3 

13.9 

0.4 

0.1 

>18 

0.6 

<35 

0.08 

<35 

8 

1 

22.2 

0.2 

0.2 

>26 

0.004 

<90 

0.02 

>26 

9 

1 

19.4 

0.1 

0.1 

>23 

0.003 

<90 

0.02 

>23 
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Table  1.  Travel  time  and  path-length  formulas 

Travel-time  and  path-length  differences  are  tn  and  Sn,  re¬ 
spectively,  between  the  first  and  nth  arrivals.  The  incremental 
travel-time  and  path-length  for  arrival  n  in  layer  i  is  indica¬ 
ted  by  tn-£,  Sn£ ,  respectively.  Corresponding  quantities  for  a 
layer  in  which  the  ray  has  a  turning  point  are  given  by  tn and 
SnT  respectively;  i.e.,  T  is  the  index  for  the  layer  containing 
a  turning  point  and  is  dependent  on  the  grazing-angle  0. 


n-1 

tn  “  *  tni  n-1  <  T 

1=1 


tni  =  (2/3gi)[Sin3  0Ri  -  (Sin2  0Ri  -  2gi  hi/Coi)3/2]  i  +  T, 

gi  f  0 

tnt  =  2  sin^  0Rt/(3  vgT) 


tnl  =  2hi  Sin  eRi/Coi  g±  =  0 


n-1 

Sn  =  s  sni  n-l  <  T 

i=l 


sni  =  (Coi/ Bi)  [Sin  eRi  “  (v±  2  -  Cos2  0Ri)1^2/vi 

+  Cos2  0Ri  ln|  (1  +  Sin  6Ri)/[vi-1  +  (v^2  -  cos2  ©Ri)1/2]  | 


i  4  T, 
gi  +  o 


Snt  =  (C0x/gx)  [Sin0RT  +  Cos2  0RT  ln|(l  +  Sin0RT)/Cos  ©RT j ] 


sni  =  2hi/Sin0R1 


gi  =  0 


1/Vi  =  [1  -  Z^hi/C^]1/2 

Cos  0Ri  =  (Coi/CW)  Cos  0O- 

Cq^  =  (CR^  CR2  •  •  •  •  CR^)  (v^  V2  •  •  •  Vi_^)  CW 

0  <  (CoT/28t)  t1  "  (CoT/CW)2  Cos20o]  <  hT  gT  +  0 
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Table  8.  Formulas  for  evaluation  of  model  parameters 


The  variable  y  is  plotted  versus  x  and  a  straight  line  is  fitted. 
Model  parameters  are  computed  from  the  slope  and  x-intercept.  For 
i  ^  T  and  ^  0,  h-^  is  computed  in  terms  of  t^ . 


Function  to  be 

plotted  in  the  Parameter  Extracted 

form:  y  =  a(x  -  xint) 


Rayleigh  Analysis 


i  -  1 


^oi  (xint) 


Arrival  Time  Analysis 


x  =  Cos^  0O 


y  =  (3  t±)  2/3 


i  =  T 


gi  =  0 


y  =  (ti)2 


Coi/CW  =  (xint)-l/2 
h±  =  CW  1  a | L/2 


2 


i  ^  T 

gi  f  o 


Cq£  and  g^  are  first  determined  at  those  angles  for 
which  i  =  T  using  formulas  above  then 


=  coi  [Sin2  eRi  -  (Sin3  eRi  -  ^ti)2/3] 

2 
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Table  9.  Formulas  for  attenuation  analysis 


Formulas  for  attenuation  analysis.  The  symbol  K±  represents  the 
average  attenuation  coefficient  in  dB/m/kHz  encountered  by  the 
ith  arrival  when  i=T.  Formulas  for  Si  are  given  in  Table  7. 


Ki  -  746.2  [1  -  (1  +  0.36  In  (1.2  Pi/Ti))1/2] 

Si 

Ti  =  (1  "  Rl2)  (1  "  R22)  •  •  •  •  (1  -  rJ_!) 

Ri  =  PRi  Sin  ©.,_!  ~  (1/CR-,)  Sin  ©R., 

DRi  Sin  +  (1/CRj^)  Sin  ©Ri 


Cos  ©Ri  =  (Coi/CW)  Cos  0O 


Cos  ©i  =  (Ci/CW)  Cos  e0 


Ci  (CRi  CR2  ....  CRi)  (vi  V2  .  .  .  .  Vi)  CW 


Coi  =  ci/vi 


Pi  is  the  amplitude  of  the  ith  arrival  of  the  unit  impulse 
response,  bandpass-filtered  with  a  Gaussian-shaped  filter 
spectrum  with  3  dB  down  points  at  43  and  87  Hz. 
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Appendix  B. 

Processed  Acoustical  Data  and  Analysis 
for  Each  Station 


Station  2 

In  the  low  frequency  traces  of  Figure 
II,  the  dominant  arrival  seen  at  graz¬ 
ing  angles  of  35°  and  less  is  inter¬ 
preted  as  a  refracted  arrival.  The  same 
type  of  arrival  may  also  be  seen  in  the 
expanded  high-frequency  traces  of  Figure 
12  at  grazing  angles  of  20.7°  and 
less.  Strong  first-arrival  pulses  are 
seen  which  may  contribute  to  low  losses 
at  high  frequencies. 

The  bottom  loss  shown  in  Figures  13  and 
14  is  characterized  at  low  frequencies 
by  a  low-loss  region  and  a  high-loss 
region.  The  angle  at  which  the  transi¬ 
tion  between  the  two  behaviors  occurs  is 
the  highest  angle  for  which  the  strong 
low-frequency  refracted  arrival  is 
observed.  At  frequencies  above  100  Hz 
the  transition  is  not  as  sharp  because 
at  high  frequencies  less  energy  is 
carried  by  the  refracted  arrival  due  to 
sediment  attenuation  along  the  refracted 
path.  Nevertheless,  the  loss  on  the 
average  increases  slightly  with  grazing 
angle;  no  simple  trend  as  a  function  of 
frequency  is  seen. 

At  high  angles  for  which  a  refracted 
arrival  is  not  observed,  the  loss  at  100 
Hz  and  below  is  on  the  order  of  10  db, 
whereas  at  frequencies  above  100  Hz,  the 
loss  Is  on  tne  order  of  5  db  with  some 
significant  fluctuation  in  the  vicinity 
of  2000  Hz.  This  suggests  that  the 
sediments  are  better  reflectors  at  high 
than  at  low  frequencies,  and  is  consis¬ 
tent  with  the  strong  reflectors  observed 
in  the  high  frequency  impulse  responses. 

Strong  fluctuation  peaks  are  observed  at 
low  frequencies  at  angles  around  25°. 
This  is  consistent  with  bottom  loss 
predictions  calculated  for  a  gradient 
layer  model.  Accordingly,  the  loss 
increases  slowly  and  smoothly  with 
grazing  angle  until  an  angle  is  reached 


grazing  angle  until  an  angle  is  reached 
at  which  the  refracted  arrival  Is  time- 
delayed  by  a  half  wavelength.  At  this 
angle,  destructive  interference  produces 
the  first  large  peak  in  the  bottom  loss  . 
The  fluctuation  pattern  then  continues 
to  all  higher  angles.  In  this  data,  the 
refracted  arrival  disappears  at  an  angle 
just  above  the  onset  of  the  fluctuation 
pattern  so  that  only  the  beginning  part 
of  the  pattern  is  seen. 

Further  evidence  for  the  interference 
hypothesis  is  seen  in:  (1)  the  strong 
peak  at  24°,  80  Hz,  and  its  companion 
at  twice  the  frequency  at  24°,  160  Hz; 
and  (2)  the  small  peak  at  20°,  100  Hz, 
and  its  strong  companion  at  20°,  200 

Hz.  Differences  in  the  magnitudes  of  the 
peaks  may  be  due  to  an  undersampling  of 
the  loss  curves  along  the  grazing  angle 
axis . 

The  signal-to-noise  ratio  averaged  over 
all  angles  for  the  100  Hz  one-third- 
octave  frequency  band  was  29  dB. 

The  Rayleigh  analysis  shown  in  Figure  15 
gives  a  density  ratio  DR=1.4  and  a  sound 
speed  ratio  CR=0.97.  The  low  frequency 
arrival  times  were  first  analyzed  as  a 
single  layer  in  Figure  16  giving  g=0.8 
sec”^-,  CR-1.05,  and  h>280  m.  Note  that 
the  data  were  plotted  both  under  the 
assumption  of  a  refracting  layer  and 
under  the  assumption  of  a  reflecting 
layer.  The  refracting  layer  assumption 
provided  a  good  linear  fit,  whereas  the 
reflecting  layer  assumption  did  not. 

Because  the  high  frequency  traces  con¬ 
sisted  primarily  of  two  dominant  arriv¬ 
als  at  low  angles,  the  analysis  was 
extended  to  lower  angles  using  arrival 
times  measured  from  the  well-resolved 
high  frequency  traces.  As  shown  in  Fig¬ 
ure  17,  the  lowest  angle  arrivals  could 
be  interpreted  as  either  a  constant 
sound  speed  layer  (CR^=1.00,  h]^=12 

m)  or  a  refracting  layer  (gi~10 
sec”l,  CR^=0.96).  The  higher  angle 
low-frequency  arrivals  were  reinterpre¬ 
ted  (Fig.  18)  on  the  basis  of  Model  5 
(g2=0.8  sec~l,  CR2=1.05,  h2  > 
260  m)  and  Model  6  (g2=0*8  sec~l, 
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Station  9 

In  the  low-frequency  impulse  responses 
of  Figure  41,  the  dominant  penetrating 
arrival  appears  at  grazing  angles  of 
40.5°  and  less.  At  higher  angles 
numerous  smaller  amplitude  pulses  are 
seen  instead.  The  high-frequency  traces 
of  Figure  41  show  that  the  dominant  ar¬ 
rival  is  made  up  of  a  complex  distribu¬ 
tion  of  reflectors.  The  expanded  high- 
frequency  traces  of  Figure  42  show  that 
the  response  is  too  complicated  for  a 
simple  model  at  high  frequencies. 

The  bottom  loss  at  low  frequencies  (Fig. 
43)  changes  gradually  from  low  losses  at 
low  grazing-angles  to  high  losses  at 
high  angles.  At  frequencies  above  1  kHz, 
the  loss  curves  (Fig.  44)  tend  to  be 
flatter  or  somewhat  concave  downward 
with  grazing-angle.  At  all  angles  and 
frequencies  the  curves  tend  to  have 
superimposed  strong  oscillation  pat¬ 
terns.  The  oscillation  pattern  at  low 
frequencies  and  angles  below  41°  is 
derived  from  the  interference  between 
two  dominant  arrivals,  and  at  higher 
angles  between  the  first  arrival  and 
multiple  arrivals  at  times  greater  than 
100  msec.  Unlike  Station  8,  at  high  fre¬ 
quencies  there  are  numerous  arrivals 
after  the  first.  As  a  result,  the  bottom 
loss  tends  to  smooth  out  at  high  angles 
and  become  more  oscillatory  at  low 
angles.  The  signal-to-noise  ratio  aver¬ 
aged  over  all  angles  for  the  100  Hz  one- 
third-octave  frequency  band  was  22  dB. 

The  Rayleigh  analysis  (DR=1.5,  CR=0.97) 
was  very  precise  as  shown  in  Figure  45 
by  the  close  clustering  of  points  about 
the  straight-line  fit.  The  arrival-time 
analysis  (CR=1.06,  h=80  m)  in  Figure  46 
also  produced  a  good  fit  on  the  basis  of 
a  reflecting  layer.  As  in  the  case  of 
Station  8  there  is  a  large  discrepancy 
between  CRR  and  CRT,  but  again,  this 
may  represent  sound-speed  changes  in  the 
uppermost  part  of  the  sediment,  which  is 
not  resolved  in  the  arrival  times  of  the 
low-frequency  traces.  The  discrepancy 
was  reduced  by  fitting  Model  7  to  the 
data  (Fig.  47:  g=4.4  sec“l ,CRi=0.97 , 

h-^=29  m,  CR2=1*00,  h2“42  m) . 


Station  7 

The  data  from  Station  7  is  strikingly 
different  from  the  rest.  The  impulse 
responses  at  both  high  and  low  fre¬ 

quencies  (Figs.  48  and  49)  are  very 
random,  i.e.,  the  reflected  pulses  do 

not  correlate  from  trace  to  trace.  In 

order  to  show  the  long  duration  of  the 
signals,  the  time  scale  in  the  figures 
has  been  compressed  eight  times  that  of 
the  other  stations.  The  expanded  high- 
frequency  traces  show  that  the  first 
arrivals  are  not  well-defined.  Note  that 
some  of  the  abrupt  signal  decays,  such 
as  seen  in  the  traces  for  gracing  angles 
11.5-12.5°  at  800  msec,  are  really  an 

artifact  of  the  processing.  In  these 
cases,  the  signal  was  not  recorded  or 
digitized  for  a  sufficient  length  of 
time . 

In  general,  the  bottom-loss  curves 
(Figs.  50  and  51)  are  nearly  flat,  but 
are  uneven  functions  of  grazing-angle  at 
all  frequencies  and  angles.  The  losses 
range  primarily  from  15  to  20  dB ,  slop¬ 
ing  gently  to  lower  losses  at  low  angles 
and  for  low  frequencies .  At  frequencies 
below  100  Hz  the  trend  suggests  a  rapid 
decrease  in  the  loss  as  the  grazing- 
angle  decreases  below  10°. 

The  data  are  suggestive  of  scattering 
from  a  very  rough  surface  with  no  de¬ 
tectable  coherent  component.  This  obser¬ 
vation  is  consistent  with  the  fact  that 
this  geophysical  province  was  a  rough, 
unsedimented,  elongated  seamount.  Be¬ 
cause  of  the  lack  of  coherence,  no 
modeling  was  attempted.  It  is  inter¬ 
esting  to  note  that  the  three  standard 
deviation  confidence  limits  on  the  bot¬ 
tom-loss  curves  collapse  to  a  single 
curve,  although  the  average  signal-to- 
noise  ratio  at  100  Hz  (29  dB)  is  not 
much  better  than,  for  example,  Station  8 
(27  dB).  The  reason  is  that  the  number 
of  degrees  of  freedom  is  much  greater 
because  of  signal  duration.  The  greater 
signal  duration  results  in  higher  fre¬ 
quency  resolution  and  hence,  more  inde¬ 
pendent  spectral  estimates  within  a  one- 
third-octave  band.  The  small  improvement 
in  signal-to-noise  ratio  in  comparison 
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with  the  other  stations  is  attributed  to 
the  smaller  dynamic  range  of  the  time 
spread  signal.  The  attenuators  or  ampli¬ 
fiers  were  set  to  a  higher  gain  to  more 
fully  utilize  the  available  dynamic 
range  of  the  tape-recorder.  In  effect 
this  suppressed  the  tape-recorder  noise 
floor  to  a  lower  level  relative  to  the 
average  signal  level. 
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Appendix  C. 

Geological  Analysis  of  Each  Station 


Comprehensive  geological  data  analyses 
for  each  station  site  are  presented  in 
this  appendix.  The  discussion  of  geo- 
logical  considerations  common  to  several 
stations  are  deferred  to  Appendix  D. 
Special  emphasis  is  placed  on  the  deriv¬ 
ation  of  several  possible  sediment 
sound-speed  profiles  as  determined  from 
available  geological  data.  The  compari¬ 
son  with  the  profiles  derived  from  the 
acoustic  data  analysis  is  shown  in  Fig¬ 
ure  10. 

Station  locations  are  defined  as  the 
location  of  the  anchored  ship,  USNS 
WILKES,  during  the  reflectivity  experi¬ 
ments.  USNS  KANE  dropped  SUS  charges 
along  a  radial  path  called  the  bomb  run. 
Seismic  profiler  and  3.5  kHz  records 
taken  along  or  near  the  bomb  runs  are 
reproduced  here.  The  sound-speed  pro¬ 
files  of  the  water  column  at  the  sta¬ 
tions  are  illustrated  in  Appendix  F. 

Station  2 

Station  2  was  located  near  the  northern 
boundary  of  the  Iberian  Abyssal  Plain  as 
shown  in  Figure  52.  The  plain  is  charac¬ 
terized  by  flat-lying,  horizontally 
stratified  sediments  characteristic  of 
abyssal  plains.  The  sediments  are  not  as 
thick  as  those  in  many  abyssal  plains, 
thus  offering  basement  exposures  in 
several  places.  Because  of  the  thinner 
sediments,  it  was  felt  that  this  plain 
would  offer  geoacoustic  properties  dif¬ 
ferent  from  the  Tagus  Abyssal  Plain 
(Stations  4  and  5)  or  the  Madeira  lower 
continental  rise  (Station  9). 

Geological  Information 

The  boundaries  of  the  Iberian  Abyssal 
Plain  are  generally  indicated  in  Figure 
52  by  the  5200  m  contour  on  the  west 
(roughly  along  15°W)  and  by  the 


4400-4800  m  contour  on  the  east  and 
south.  The  plain  has  a  general  slope 
from  northeast  to  southwest  with  an 
average  slope  of  1:3600.  Reynolds  [46] 
compiled  a  geological  analysis  of  the 
whole  abyssal  plain,  emphasizing  the 
area  at  Station  2.  Much  of  the  following 
information  is  taken  from  his  report. 

A  seismic  profiler  record  taken  on  a 
1971  KANE  cruise  near  Station  2  (Fig. 
52)  is  shown  in  Figure  53;  Station  2 
corresponds  to  a  position  on  the  record 
near  the  break  at  2300Z.  The  abyssal 
plain  is  characterized  by  several  fea¬ 
tures  visible  on  this  profile.  Between 
0500  and  0700Z,  a  thick  section  of  sedi¬ 
ment  exists  as  a  transition  zone  between 
the  fault-block  to  the  west  and  the 
rough  flat  block  relief  of  the  continen¬ 
tal  margin  to  the  east.  This  feature  is 
observed  all  along  the  eastern  boundary 
of  the  plain,  varying  between  5  and  40  m 
in  width.  In  all  cross  sections  of  this 
feature,  the  surface  is  nearly  smooth; 
in  some,  reflectors  can  be  traced  from 
the  abyssal  plain  into  this  zone  with  no 
loss  of  continuity,  while  in  other  sec¬ 
tions,  continuity  is  broken. 

Beneath  the  stratified  flat-lying  sedi¬ 
ments,  the  very  rough  acoustic  basement 
can  be  seen  to  the  left  (west)  of  2300Z 
in  Figure  53,  varying  between  0.8  and 
2.0  sec  two-way  travel  time  below  the 
sea  floor.  This  roughness  is  character¬ 
istic  of  the  basement  throughout  the 
abyssal  plain  as  it  rises  to  the  west 
approaching  the  flanks  of  the  Mid- 
Atlantic  Ridge.  At  about  40°15fN  the 
acoustic  basement  descends  to  the  south¬ 
east  with  the  overlying  sediment  forming 
a  wedge  of  over  2.0  sec  of  thickness 
which  increases  to  the  southeast  (see 
Fig.  52).  An  arc  of  sediment  thicker 
than  1.0  sec  extends  from  the  southeast¬ 
ern  corner  of  the  plain  along  the  east¬ 
ern  side  and  across  the  northern  end  of 
the  plain.  Sediment  in  the  southwest 
corner  and  the  west  central  part  of  the 
abyssal  plain  is  generally  less  than  one 
second  thick  and  has  many  basement  expo¬ 
sures.  At  a  latitude  of  about  39°40’N, 
the  basement  emerges  from  the  abyssal 
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plain  to  form  an  east-to-west  topograph¬ 
ic  high  along  the  southern  edge  of  the 
Iberian  Abyssal  Plain. 

Two  main  sedimentary  units  are  apparent 
in  the  Iberian  Abyssal  Plain,  and  are 
separated  by  a  prominent  reflector  that 
is  continuous  throughout  the  basin.  The 
sedimentary  unit  above  this  reflector 
varies  from  a  thickness  of  0.3  sec  in 
the  northwest  to  greater  than  1.0  sec  in 
the  east  and  southeast.  The  prominent 
reflector  is  flat-lying  in  the  north, 
begins  to  dip  to  the  southeast  at  about 
42°N,  and  becomes  less  reflective 
until  it  is  no  longer  discernible  in  the 
south.  In  many  areas  of  the  Iberian 
Abyssal  Plain,  the  upper  sedimentary 
unit  (Unit  1)  lies  unconformably  on  the 
major  reflector  which  marks  the  top  of 
the  underlying  Unit  2  (Fig.  53).  Both 
sedimentary  units  comprise  interbedded 
layers  of  transparent  (pelagic)  and 
reflective  (turbidite)  layers;  however, 
the  upper  unit  is  flat-lying,  but  the 
lower  stratified  unit  is  characterized 
by  faulting  and  enough  nonhorizontal 
orientation  in  the  conformable  reflec¬ 
tors  to  suggest  tectonic  deformation  of 
this  unit  after  deposition. 

Within  the  upper  unit  are  two  distinct 
sections.  The  upper  0.2  sec  (Section  A) 
contains  a  sequence  of  highly  reflective 
layers;  weaker  reflectors  below  give  the 
lower  section  (B)  a  more  transparent 
character.  Section  B  is  observed  only  in 
the  northern  part  of  the  Iberian  Basin. 
To  the  south,  nearing  the  source  of 
sediments  at  Nazare  Canyon,  reflectors 
of  Section  B  gradually  become  more  prom¬ 
inent,  until  at  42°N  no  distinction 
can  be  made  between  Sections  A  and  B. 

Available  3.5  kHz  data  were  analyzed  to 
determine  the  number  of  reflectors  per 
meter,  continuity  of  reflectors,  dip, 
and  depth  of  penetration.  Analysis  of 
3.5  kHz  data  was  hindered  because  the 
depth  in  the  northern  part  of  the  abys¬ 
sal  plain  and  a  portion  of  the  southern 
part  fall  in  the  depth  range  where  the 
return  signal  is  gated  to  eliminate 
interference  by  the  outgoing  pulse. 
Because  of  the  poor  quality  of  the  gated 


records,  no  3.5  kHz  record  is  presented 
for  this  station.  Measurements  made  to 
determine  the  number  of  reflectors  per 
meter  for  the  remaining  areas  in  the 
southern  Iberian  Basin,  however,  give  a 
mean  value  of  0.297  with  a  range  of 
0.267-0.320  (3.37  m  between  reflectors, 
range  3.74-3.12).  As  discussed  in  Appen¬ 
dix  D.2,  this  value  indicates  fewer  re¬ 
flectors  per  meter  than  actually  exist. 

Two  core  samples  taken  north  of  42°  in 
the  plain  consist  of  firm  clay  and  firm 
silty  clay,  both  with  a  carbonate  con¬ 
tent  of  30-40%.  Other  coring  attempts  in 
the  northern  Iberian  Abyssal  Plain  fail¬ 
ed  to  penetrate  the  sea  floor,  but  had 
highly  compacted  clay  samples  in  the 
core  cutters.  The  R/V  VEMA  (Lamont- 
Doherty)  core  taken  near  41°N,  15°W 
shows  several  coarse-grained  turbidite 
layers,  10-15  cm  thick;  a  USNS  KANE  core 
at  the  plain* s  northern  edge  contained 
one  10  cm  turbidite  (Fig.  52).  Core 
studies  show  that  the  pelagic  material 
in  the  Iberian  Basin  consists  of  layers 
of  f o rami ni feral  lutite  with  moderate  to 
moderately  high  calcium  carbonate  con¬ 
tent.  The  turbidites  are  made  up  of  ter¬ 
rigenous  material  and  shallow  water 
detritus  derived  from  the  western  Euro¬ 
pean  land  mass  and  continental  shelf, 
grading  upward  from  coarse-grained  basal 
layers  to  fine-grained  clay  or  silty 
clay.  Contact  of  the  turbidites  with  the 
layer  below  is  very  distinct,  but  the 
upper  contact  is  imperceptible. 

The  sediments  in  the  northern  part  of 
the  Iberian  Abyssal  Plain  are  derived 
predominantly  from  the  Bay  of  Biscay. 
Turbidity  currents  flowing  across  the 
Bay  of  Biscay  are  rejuvenated  as  they 
reach  the  steeper  gradient  of  Theta  Gap 
and  flow  southwestward  into  the  Iberian 
Abyssal  Plain  [47].  Because  of  the  great 
distance  from  its  source,  these  turbid¬ 
ites  would  be  expected  to  have  fine¬ 
grained  basal  layers  (silt  to  sandy 
silt) . 

Nazare  Canyon,  southeast  of  the  plain, 
is  the  major  supplier  of  turbidite 
sediment  to  the  entire  Iberian  Abyssal 
Plain.  Turbidite  layers  nearest  this 
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source  would  be  expected  to  have  coarse¬ 
grained,  poorly  sorted  basal  units.  Pro¬ 
gressing  away  from  the  source,  the  tur- 
bidites  should  become  thinner,  finer 
grained,  and  better  sorted.  The  increase 
in  mean  grain  size  close  to  the  source 
area  suggests  that  greater  impedance 
contrasts  should  be  found  in  the  south¬ 
east  portion  of  the  plain.  This  is  sup¬ 
ported  by  3.5  kHz  and  seismic  records. 

Laughton  [48]  states  that  the  overspill 
of  the  Bay  of  Biscay  did  not  start  until 
the  time  represented  by  a  horizontal 
layer  73  m  below  the  sea  bed.  If  Laugh¬ 
tons  theory  is  correct,  then  almost  all 
of  the  Iberian  Abyssal  Plain  sediment 
other  than  about  the  upper  50-100  m  was 
supplied  by  Nazare  Canyon.  However, 
Reynolds  [46]  notes  evidence  that  a 
thicker  portion  of  the  sedimentary  col¬ 
umn  is  derived  from  the  Biscay  Abyssal 
Plain.  Profiler  records  obtained  west  of 
Theta  Gap  and  in  the  Biscay  Abyssal 
Plain  on  other  NAVOCEANO  cruises  suggest 
a  continuity  of  the  interface  between 
Units  1  and  2  with  a  strong  reflecting 
horizon  in  the  Biscay  Abyssal  Plain. 
This  horizon,  identified  as  the  mid-R 
reflector  [49],  was  drilled  at  DSDP  Site 
118  in  the  Biscay  Abyssal  Plain  (Fig. 
52).  On  the  profiler  records,  the  sedi¬ 
ments  above  the  mid-R  reflector  look 
similar  to  Unit  1  sediments  at  Station 
2.  The  sedimentary  sequence  above  the 
mid-R  reflector  consists  of  a  250  m 
section  of  flat-lying,  highly  reflective 
turbidites  overlying  a  150  m  thick  sec¬ 
tion  of  weaker  reflectors  in  nearly 
transparent  sediment.  These  sections 
appear  to  correspond  with  Sections  A  and 
B,  respectively,  at  Station  2.  If  Laugh¬ 
tons  theory  of  the  time  of  overspill 
into  the  Iberian  Abysssal  Plain  is  in¬ 
correct,  then  possibly  the  age  of  Unit  1 
can  be  correlated  with  sediments  above 
the  mid-R  reflector. 

The  seismic  record  near  Station  2  in  the 
Iberian  Basin  shows  Section  A  to  be  0.2 
sec  thick,  or  about  two-thirds  as  thick 
as  the  upper  section  at  Site  118.  If  the 
lower  interface  of  Section  A  represents 
the  onset  of  glaciation  (2-3  mybp)  when 
turbidity  flows  were  greatly  increased. 


then  the  sedimentation  rate  of  this 
section,  consisting  of  Pleistocene  sedi¬ 
ment,  is  about  5-8  cm/1000  yrs.  This 

rate  is  about  half  that  of  the  corre¬ 
sponding  sediments  at  Site  118  in  the 
Biscay  Abyssal  Plain.  This  lower  rate  is 
reasonable,  since  the  northern  Iberian 
Basin  is  far  from  turbidite  sources  both 
to  the  northeast  and  the  southeast  as 
compared  to  Site  118.  Section  B  corre¬ 
sponds  to  the  underlying,  more  trans¬ 
parent,  sediments  at  Site  118  of  Plio¬ 
cene  Age. 

Sediment  Surface  Sound-Speed  and  Density 

Sound-speed  measurements  made  on  the 
pelagic  sections  of  several  cores  show  a 
sediment  to  bottom  water  sound-speed 
ratio  of  0.974.  Based  on  the  core  and 
3.5  kHz  data,  about  15  cm  of  reflective 
turbidite  layers  (i.e.,  the  coarse¬ 
grained  basal  layers)  with  an  average 
sound  speed  ratio  of  1.100  (1700  m/sec) 
are  estimated  in  the  upper  5  m  of  sedi¬ 
ment.  Using  these  values,  the  average 
sediment  to  bottom  water  sound-speed 
ratio  is  about  0.977,  with  a  sound-speed 
of  1515  m/sec.  Using  the  same  propor¬ 
tions  with  average  densities  of  1.5  and 
1.8  g/cm^  for  the  clay  and  the  turbid¬ 
ites,  respectively,  the  estimated  aver¬ 
age  density  is  1.51  g/cm^. 

Sound-Speed  Profile 

An  estimate  of  a  constant-gradient 
sound-speed  profile  for  Station  2  is 
based  on  the  assumption  that  the  sedi¬ 
ments  are  analogous  to  those  at  DSDP 
Site  118  in  the  Biscay  Abyssal  Plain. 
The  method  used  (CRIC)  and  assumptions 
are  discussed  in  Appendix  E.2.  The  cal¬ 
culated  average  gradient  is  0.67  sec”1 
throughout  the  observed  0.48  sec  thick¬ 
ness  of  sediments  above  the  mid-R  re¬ 
flector.  This  value  for  an  average  grad¬ 
ient  is  the  most  accurate  possible  based 
on  available  geological  and  geophysical 
data . 

Three  other  sound-speed  profiles  esti¬ 
mated  in  different  ways  for  sediments  at 
Station  2  are  also  shown  in  Figure  54. 


Ill 


All  the  profiles  have  been  adjusted  to 
the  surface-sediment  sound  speed  of  1515 
m/sec.  Each  profile  is  drawn  to  a  depth 
corresponding  to  the  two-way  travel  time 
of  0.48  seconds. 

Cernock  [50]  has  measured  the  sound 
speed  of  clay  sediments  during  consoli¬ 
dation  tests  in  the  laboratory  (Appen¬ 
dix  E.l).  The  sample  used  for  the  data 
presented  in  Figure  54  is  from  a  core 
taken  from  an  abyssal  plain  in  the  Gulf 
of  Mexico  near  the  continental  rise. 

The  Houtz  and  Ewing  [39]  curve  (Appendix 
E.3)  is  derived  from  data  from  sediments 
which  include  types  similar  to  those 
found  at  Station  2,  but  are  located  in 
the  western  Atlantic. 

Hamilton  et  al.  [51]  derived  a  sound- 
speed  gradient  for  turbidite  sediments 
using  sonobuoy  wide-angle  reflection 
data  (see  Appendix  E.4).  The  average 
gradient  decreases  with  depth  with  an 
overall  average  value  of  1.05  sec“^. 

In  summary,  four  methods  for  estimating 
the  sound-speed  profile  in  sediments  may 
be  applicable  to  the  northern  Iberian 
Abyssal  Plain.  In  deriving  the  sound- 
speed  profile  for  Site  118  (assuming  a 
similar  profile  at  Station  2),  the  time/ 
depth  parameters  are  the  only  fixed  data 
points.  For  simplicity,  a  constant  grad¬ 
ient  is  assumed.  However,  the  other 
three  methods  suggest  decreasing  grad¬ 
ients  with  depth.  The  two-way  travel 
time  to  400  m  for  Cernock1 s  [50]  profile 
and  the  constant-gradient  profile  are 
within  8  msec,  which  is  within  the  error 
limits  of  the  time  measurement  on  the 
profiler  record  at  Site  118.  Although 
Cernock !s  model  fits  the  observed  param¬ 
eters,  insufficient  work  has  been  done 
in  the  shallow  sediments  to  establish 
how'  to  modify  accurately  the  constant 
gradient  model. 

Attenuation 

Actual  measurements  of  mean  grain  size 
for  cores  in  this  area  are  not  avail¬ 
able.  Based  on  the  core  data  available, 
an  estimated  3%  of  the  core  consists  of 


1700  m/sec  velocity  material.  Referring 
to  Hamilton’s  [52,  45]  sound-speed  and 
mean  grain-size  relationships,  material 
with  this  sound  speed  has  a  grain  size 
of  about  3  In  Figure  3  of  Reference 
[45],  the  value  of  the  attenuation  con¬ 
stant  k  at  i>  =  3  is  0.500  dB/m/kHz.  The 
value  of  0.500  may  be  a  good  average  for 
the  estimated  3%  of  coarse  turbidite 
material  in  the  cores.  This  value,  aver¬ 
aged  with  the  rest  of  the  fine  core 
material  with  K  =  0.065,  gives  a  value 
of  K=  0.078  dB/m/kHz. 

Stations  4  and  5 

Stations  4  and  5  are  located  in  the 
southern  Tagus  Abyssal  Plain  about  130 
km  apart  (Fig.  55).  These  sites  offer  a 
thick  sedimentary  cover  (up  to  2  km) 
over  the  acoustic  basement.  Despite  the 
relatively  small  separation  of  the  two 
stations  within  a  single  geophysical 
province,  sufficient  contrast  in  geo¬ 
acoustic  parameters  exists  to  justify 
bottom-reflectivity  measurements  at  each 
site. 

Geological  Information 

Ruddiman  and  Glover  [53]  studied  core 
data  and  reflection  records  in  the  Tagus 
Abyssal  Plain  in  order  to  define  the 
geoacoustic  parameters  in  the  vicinity 
of  Stations  4  and  5.  The  following  in¬ 
formation  is  taken  from  that  report. 

Sedimentary  Composition  and  Structure 

The  sediment  cores  indicated  in  Figure 
55  were  taken  by  Ruddiman  on  USNS  WILKES 
in  October  and  November  1973. 

The  cores  contained  the  two  basic 
sediment  types  which  have  alternately 
been  deposited  in  the  basin:  (1)  fine¬ 
grained  silty  clays  with  moderate  car¬ 
bonate  content  (pelagic  depositions 
accumulating  at  2-3  cm/1000  yrs),  and 

(2)  well-sorted  fine  sand  or  silt  grad¬ 
ing  upward  through  progressively  finer 
grain  sizes  into  silty  clays.  Duplaix  et 
al.  [54]  have  shown  that  the  heavy  frac¬ 
tion  of  the  mineral  sands  and  silt  in 
Tagus  Abyssal  Plain  cores  are  virtually 
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identical  in  mineralogical  composition 
to  the  sands  and  silts  found  on  the 
continental  shelf  of  Portugal  near  the 
mouth  of  the  Tagus  Fiver.  The  sediments 
on  the  shelf  have  been  periodically  and 
rapidly  transported  through  two  subma¬ 
rine  canyons  to  the  Tagus  Abyssal  Plain 
and  deposited  as  turbidite  layers. 

Figures  56  and  57  show  3.5  kHz  records 
taken  during  the  bomb  runs  for  Stations 
4  and  5.  Two  tracks  with  seismic  pro¬ 
filer  coverage  of  the  Tagus  Abyssal 
Plain  are  shown  in  Figure  55.  Although 
no  tracklines  actually  traverse  the 
length  of  the  bomb  runs,  the  profiles, 
shown  in  Figure  58,  are  representative 
of  the  sediment  profile  along  both 
tracks.  The  two  types  of  sediments  col¬ 
lected  in  the  cores  represent  the  alter¬ 
nating  transparent  and  reflective  se¬ 
quence  necessary  to  produce  the  3.5  kHz 
near-surface  bottom  reflections  and  the 
seismic-profile  deep  reflectors  visible 
over  the  southern  Tagus  Abyssal  Plain. 

Ruddiman  and  Glover  [53]  studied  in 
detail  the  reflectors  in  the  cores  and 
reflection  records  for  comparison  pur¬ 
poses.  Their  discussion  is  summarized  in 
Appendix  D.2.  The  number  of  reflectors 
per  meter  of  sediment  is  about  the  same 
on  the  3.5  kHz  records  throughout  the 
whole  Tagus  Basin;  the  reflectors  are 
often  traceable  for  40-80  km  in  the 
southern  part  of  the  plain,  but  are 
less  continuous  to  the  north.  Groups  of 
reflectors  seen  on  the  seismic  records 
seem  to  be  continuous  throughout  the 
entire  Tagus  Abyssal  Plain.  Character¬ 
istics  of  the  3.5  kHz  records  along  the 
bomb  runs  can  be  described  as  follows. 


Station  4:  Poor  penetration  (averaging 
20  m);  one  reflector  continuously  trace¬ 
able  throughout  the  run;  thickening  of 
layers  near  middle  of  the  track;  and  a 
dip  in  reflectors  up  to  0.4  m/km  toward 
the  southern  end.  Station  5;  Better 
penetration,  averaging  30-35  m;  good 
continuity  of  reflectors;  dip  toward 
middle  of  the  trackline,  averaging  0.1 
m/km;  and  thickness  between  reflectors 
increasing  to  northeast  at  mean  rates  of 
0.3%/km. 

These  dip  and  thickening  measurements 
along  the  bomb  runs  are  local  two-dimen¬ 
sional  trends  superimposed  on  the  gene¬ 
ral  geometry  of  the  abyssal  plain,  which 
has  a  general  dip  and  thickening  of  the 
layers  to  the  southeast.  Thus,  the 
deeper  layers  observed  in  the  seismic 
records  have  greater  dips  (relative  to 
the  horizontal)  than  do  the  near-surface 
layers.  As  in  the  3.5  kHz  records,  the 
deeper  reflectors  tend  to  dip  toward  the 
middle  of  the  run  for  Station  5.  The 
dips  observed  are  summarized  in  Table 
10. 

Two  major  sedimentary  units,  each  great¬ 
er  than  1.0  sec  thickness,  are  visible 
in  the  seismic  profiles.  The  deepest 
unit  comprises  interbedded  reflectors 
and  transparent  layers  in  conformable 
sequences  which  are  highly  deformed. 
Evidently  this  unit  represents  an  anci¬ 
ent  episode  of  turbidite/pelagic  deposi¬ 
tion  which  was  later  deformed  by  tecton¬ 
ic  activity. 

The  upper  layer  shows  flat-laying,  high¬ 
ly  reflective  turbidites  with  great 
regional  continuity.  This  unit  forms  a 


Table  10.  Observed  dips  of  reflectors  at  stations  4  and  5 

Dips  Along  Bomb  Run 


Layer  Depths 

Near  Surface 

0.4-0. 5  sec  (300-375  m) 

0.8-0. 9  sec  (650-745  m) 


Station  4 
0.1  m/km 
0.5-0. 8  m/km 

1.1-1. 6  m/km 


Station  5 
0.1  m/km 

1.1  m/km 

3.2  m/km 


Regional  Dip 

0. 2-0.3  m/km 

1.1- 1. 6  m/km 

2. 2- 3. 2  m/km 
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southward  thickening  and  dipping  wedge 
of  surprising  geometric  regularity.  The 
thickness  of  this  unit  reaches  maximum 
values  of  1.2  sec  (1050  m)  in  the  south¬ 
ern  Tagus  Abyssal  Plain  and  is  roughly 
1.0  sec  (850  m)  beneath  acoustic  Sta¬ 
tions  4  and  5. 

Age  of  Sediments 

Extrapolating  from  core  data,  a  Pleisto¬ 
cene  sedimentation  rate  (combined  pelag¬ 
ic  and  turbidite  input)  is  10  cm/1000 
yrs.  As  Pleistocene  sedimentation  rates 
are  generally  about  twice  that  of  ear¬ 
lier  rates,  an  average  for  the  whole 
upper  unit  is  estimated  to  be  5-10  cm/ 
1000  yrs.  Calculations  show  the  age  at 
the  bottom  of  the  upper  1000  m  wedge  of 
turbidites  to  be  about  10-20  million 
years  (my).  The  unconformity  between  the 
upper  flat-lying  turbidite  wedge  and  the 
deformed  lower  turbidite  unit  is  assumed 
to  be  a  long  hiatus  which  ended  about  15 
million  years  before  the  present  (mybp) 
(based  on  findings  on  such  hiatuses  in  a 
variety  of  DSDP  sites  [55,  56]  in  the 

eastern  North  Atlantic). 

The  deeper  sediment  layers  are  far  more 
difficult  to  date.  Magnetic  anomaly  dat¬ 
ing  shows  that  the  age  of  the  basement 
material  of  the  basin  from  west  to  east 
is  100-120  my  old,  and  115  my  for  the 
region  of  Stations  4  and  5.  Ruddiman  and 
Glover  [53]  speculate  that  the  lower 
deformed  turbidite  unit  probably  dated 
roughly  70-40  mybp  and  overlies  a  200  m 
pelagic  Cretaceous  layer  (aproximately 
70-115  mybp). 

As  no  DSDP  holes  have  been  drilled  in 
the  Tagus  Basin,  data  were  used  from  two 
sites  with  similar  settings  and  sedi¬ 
ments  to  estimate  the  ages  and  composi¬ 
tions  of  the  major  sedimentary  sequences 
(Site  118  [49]  in  the  Biscay  Abyssal 

Plain  [see  Fig.  52]  and  Site  142.  [57]  in 
the  Ceara  Abyssal  Plain  in  the  equator¬ 
ial  Atlantic  off  Brazil).  The  uppermost 
sediment  sequence  at  Site  118  is  a  250  m 
sequence  of  Pleistocene  (0-2  mybp)  tur¬ 
bidites  which  have  basal  silts  and 
sands,  are  interbedded  with  silty  clays, 
and  were  deposited  at  a  combined  rate  of 


12.5  cm/1000  yrs.  From  250  to  400  m,  the 
sediment  sequence  consists  of  Pliocene 
(2-5  mybp)  clays  and  silts  deposited  at 
rates  of  5  cm/1000  yrs.  The  Biscay  Abys¬ 
sal  Plain  sediments  thus  serve  as  a  good 
analogue  to  the  upper  portions  of  the 
Tagus  wedge.  At  Site  142,  the  last  15  my 
of  sediment  comprise  700  m  of  alternat¬ 
ing  turbidite  and  pelagic  material  which 
was  deposited  at  rates  in  excess  of  10 
cm/1000  yrs.  The  sediments  from  400  to 
700  m  at  Site  142  form  a  better  analogue 
to  the  Tagus  record  than  those  at  Site 
118. 

Sediment  Surface  Sound-Speed  and  Density 

Sound-speeds  were  measured  on  cores 
taken  by  USNS  WILKES  (Fig.  59).  The 
pelagic  sediments  have  sound-speeds 
(adjusted  to  laboratory  conditions)  of 
about  1497  m/sec  and  densities  of  about 

1.5  g/cm^.  The  finer-grained  upper 
portions  of  the  turbidites  have  similar 
characteristics;  they  can  be  distin¬ 
guished  from  pelagics  by  color  and 
structure  but  are  acoustically  identi¬ 
cal.  Both  types  are  acoustically  trans¬ 
parent.  The  coarse-grained  basal  layers 
in  the  turbidites  are  acoustically 
reflective;  they  have  velocities  of 
1650-1800  m/sec  and  densities  of  1.8-1. 9 
g/cm^.  The  basal  turbidite  layers  can 
be  considered  simplistically  as  constant 
gradient  layers:  their  bases  form  step¬ 
like,  discontinuous  increases  of  sound 
speed  (from  1500  to  1700  m/sec)  and  den¬ 
sity  (from  1.5  to  1.9  g/cm^)  from  the 
underlying  pelagic  layers,  and  their 
tops  grade  linearly  into  the  sound 
speeds  and  densities  of  the  upper  turbi¬ 
dite  sections,  which  are  identical  to 
the  overlying  pelagic  layers.  Based  on 
core  data,  the  mean  thickness  of  reflec¬ 
tive  basal-turbidite  layers  is  about 
20  cm. 

The  mean  intervening  thickness  of  the 
low-speed  transparent  sediment  is  rough¬ 
ly  100  cm  at  Station  4,  and  80  cm  at 
Station  5.  The  average  speed  and  density 
for  the  upper  5  m  have  been  calculated 
for  both  stations  using  this  information 
as  shown  in  Table  11.  The  average  sound 
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speeds  listed  in  Table  11  have  been  cor¬ 
rected  to  in  situ  conditions. 

Sound-Speed  Profile 

Sound-speed  changes  per  unit  depth 
through  the  sedimentary  column  near  the 
acoustic  stations  in  the  Tagus  Abyssal 
Plain  are  inferred  through  indirect 
means.  This  is  accomplished  by  using 
relevant  DSDP  measurements  of  sound- 
speed  in  sediments  with  ages,  lithol¬ 
ogy,  and  inferred  state  of  compaction 
comparable  to  Tagus  sediments. 

Based  on  the  study  by  Ruddiman  and 
Glover  [53]  the  sound-speed  profile  in 
the  upper  200  m  of  the  Tagus  Abyssal 
Plain  is  assumed  to  be  similar  to  that 

at  Site  118  in  the  Biscay  Abyssal  Plain. 
Using  the  CRIC  method  described  in 
Appendix  E.2,  the  calculated  gradient  is 
0.67  sec  for  a  400  m  section.  As 
discussed  in  the  appendix,  the  gradient 
in  the  upper  200  m  is  probably  slightly 
higher.  The  curve  is  shown  in  Figure  54, 
Two  curves  are  shown  in  Figure  10  for 

Stations  4  and  5.  In  addition  to  the 

Site  118  curve,  Hamilton1  s  [51]  curve 
for  abyssal  plains,  as  discussed  in 
Appendix  E.4  for  turbidite  sediments,  is 
shown.  The  two  curves  commence  at  the 

sediment  interface  with  the  sound-speed 
given  in  Table  11. 


Attenuation 

Grain-size  measurements  were  taken  on 
core  W5  at  Station  4  and  core  W4  at  Sta¬ 
tion  5.  The  sampling  interval  was  suf¬ 
ficient  on  W4  to  provide  a  good  estimate 
of  mean  grain  size  throughout  the  core. 
The  closely  spaced  sound-speed  measure¬ 
ments  confirm  the  thoroughness  of  the 
grain-size  sampling.  On  core  W5,  grain- 
size  was  not  determined  in  the  high 
sound-speed  layers,  so  the  same  average 
grain-size  as  that  in  W4  was  assumed. 
The  percentages  of  basal  turbidite  and 
intervening  layer  thickness  are  the  same 
as  those  used  in  Table  11.  The  attenua¬ 
tion  determinations  k  q  based  on  Hamil¬ 
ton1  s  [45]  relationships  are  shown  in 
Table  12.  Note  that  the  effect  on  the 
average  calculated  attenuation  due  to 
the  lower  percentages  of  basal  turbidi¬ 
ties  at  Station  4  rather  than  at  Station 
5  is  offset  by  the  slightly  finer  grain- 
size  at  Station  5. 

Station  8 

The  site  for  Station  8  offers  a  geologic 
setting  quite  different  from  that  of  the 
other  stations.  The  site  is  located  on  a 
plateau  contiguous  to  Dragon  Seamount 
and  removed  from  the  continental  margin. 
The  locations  of  Station  8  and  the  bomb 
run  are  shown  in  Figure  60.  The  station 


Table  11.  Sound- 

-speed  ratio  and  density  at  stations  4  and  5 

Station 

Sediment  Type 

% 

Sound-Speed  m/sec 
(Ratio) 

Mean  Value  Average 

Density  g/cm^ 

Mean  Value  Averag< 

4 

Basal  Turbidite 

16% 

1600 

(1.051) 

1539 

(0.994) 

1.7 

1.53 

Intervening 

Sediments 

84% 

1497 

(0.984) 

1.5 

5 

Basal  Turbidite 

20% 

1600 

(1.051) 

1542 

(0.996) 

1.7 

1.54 

Intervening 

Sediments 

80% 

1497 

(0.984) 

1.5 
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Table  12.  Sound  attenuation  from  cores  at  stations  4  and  5 


Core  Sediment  Type 

W5  Basal  Turbidite 

Station  4 

Intervening 

W4  Basal  Turbidite 

Station  5 

Intervening 


is  at  a  relatively  shallow  depth  of 
about  2400  m.  The  run  traverses  a  change 
of  depth  from  2460  to  2057  m.  The  sedi¬ 
ments  are  sufficiently  thick  to  provide 
acoustic  propagation  paths  through  the 
ocean  floor,  but  lack  the  bathymetric 
smoothness  found  in  the  abyssal  plains. 
As  can  be  seen  in  the  3.5  kHz  record 
shown  in  Figure  61,  the  topography  has 
slight  relief  of  about  25  m  and  maximum 
slopes  of  3°.  The  site  was  chosen  pri¬ 
marily  because  of  the  contrast  in  sedi¬ 
ments  with  the  deeper  stations  and  also 
because  of  the  slightly  rough  topogra- 
phy. 

Geological  Information 

Sediment  Thickness 

Seismic  records  were  obtained  along  two 
tracks  nearly  parallel  to  the  bomb  run 
(Fig.  60).  The  records  shown  in  Figure 
62  indicate  that  the  sediment  thickness 
is  at  least  0. 7-1.0  sec  two-way  travel 
time  over  the  entire  bomb  run.  The  deep¬ 
est  reflector  observed  along  the  KANE 
track  is  relatively  smooth  and  has  a 
sufficiently  high  acoustic  impedance  to 
mask  any  rough  basement,  whereas  the 
basement  roughness  is  more  apparent  on 
the  WILKES  track  farther  west.  Diffuse, 
but  strong,  reflectors  at  the  seismic 
frequencies  are  observed  about  two- 
thirds  of  the  way  down  the  sedimentary 
column. 


Average 


% 

Grain-Size 

i 

kg 

dB/m/kHz 

16 

6.2 

0.137 

84 

8.8 

0.064 

Average 

0.076 

20 

6.2 

0.137 

80 

9.1 

0.062 

Average  0.077 


The  sediment  thickness  on  the  plateau  is 
well-defined  by  other  NAV0CEAN0  tracks 
traversing  the  area.  The  side  of  the 

plateau  has  a  gentle  slope  to  the  north¬ 
west,  with  increasing  slope  continuing 
counter-clockwise  around  the  periphery 

of  the  plateau.  The  blanket  of  sedi¬ 
ments  has  continuity  of  thickness  and 
appearance  on  traverses  crossing  the 

plateau  to  the  northwest  of  Station  8. 
As  the  slope  increases,  the  sediments 
thin  to  the  extreme  of  no  observable 
sediments  on  the  southern  side  of  the 
plateau,  as  revealed  in  a  seismic  pro¬ 
file  record  (not  shown)  on  a  traverse 
between  Station  8  and  Site  136  to  the 
southeast.  The  ridge  at  point  Kf  in 

Figures  60  and  62  is  the  southern  bound¬ 
ary  of  thick  sediments. 

Sediment  Composition 

Sediments  in  the  general  vicinity  of 
Station  8  are  predominantly  calcareous, 
and  pelagic  in  origin.  Linered  Ewing 
cores,  SB-3  and  SB-4,  were  taken  at  Sta¬ 
tions  7  and  8,  respectively.  Core  SB-3 
penetrated  the  entire  sedimentary  cover 
of  only  82  cm.  Although  the  sedimentary 
cover  at  Station  8  is  thick,  SB-4  only 
penetrated  112  cm  despite  the  use  of  a 
900  kgm  (1  ton)  weightstand,  possibly 
because  of  the  presence  of  at  least  one 
layer  of  coarse  sediment.  The  grain  size 
distribution  was  measured  on  two  samples 
from  each  core.  The  sample  at  22-28  cm 
depth  in  SB-4  reveals  a  coarse  layer 
with  77%  of  the  constituents  in  the  sand 
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size  (>0.063  mm).  The  other  three  sam¬ 
ples  in  the  two  cores  have  sandy-silty 
clay.  Both  cores  consist  of  predominant¬ 
ly  foraminiferal  ooze. 

A  Lamont-Doherty  core,  V30-193,  taken  18 
km  north  of  Station  8  in  a  water  depth 
of  2842  m  (Fig.  60),  consists  predomi¬ 
nantly  of  foraminiferal  marl  ooze  in  its 
448  cm  length.  The  carbonate  content  is 
moderate  to  high.  The  only  detritus 
observed  is  negligible  amounts  of  dark 
minerals  and  metamorphic  rock  fragments. 
Core  V30-192  was  taken  62  km  northwest 
of  Station  8  in  a  depth  of  3808  m  (north 
of  the  area  of  Fig.  60).  At  least  340  cm 
of  sediment  were  obtained  which  consist 
also  of  predominantly  carbonate  sedi¬ 
ments  with  negligible  dark  minerals. 

The  elevation  of  the  site  precludes 
recent  turbidity-current  deposition  of 
terrigenous  sediments  from  the  conti¬ 
nents.  Any  allochthonous  recent  sed¬ 
iments  would  be  derived  from  turbidity 
flows  off  the  Dragon  Seamount  (Station 
7),  from  ice  rafting,  or  from  windblown 
detritus  (including  ash  layers).  The 
seamount  turbidites  would  predominantly 
consist  of  pelagic  sediments  with  some 
fragments  of  basaltic  rock. 

Volcanism 

Uchipi  et  al .  [58]  report  post-Miocene 

volcanism  in  Madeira  and  volcanism  com¬ 
mencing  in  the  Miocene  up  through  the 
present  in  the  Canary  Islands.  Evidence 
of  volcanism  can  be  seen  in  the  3.5  kHz 
and  seismic-profiler  records.  A  well- 
defined  3.5  kHz  reflector  at  14  m  below 
the  sea  floor  can  be  traced  along  the 
entire  record  in  Figure  61.  This  reflec¬ 
tor  is  observed  on  the  tracks  near  Sta¬ 
tion  8  as  far  south  as  32°50fN,  and  as 
far  east  as  16°W.  This  reflector  may 
be  an  ash  layer  from  very  recent  (0.5 
mybp)  volcanism. 

Evidence  of  volcanism  is  also  observed 
in  the  seismic-profiler  record  near  the 
outcrop  (Dragon  Seamount).  Intense  heat¬ 
ing  of  the  sediments  adjacent  to  the 
volcanic  outcrop  during  the  time  of 


volcanism  may  have  resulted  in  an 
increase  in  the  acoustic  impedance  of 
the  sediments  (F.  Bowles,  pers.  comm.). 
Sediments  deposited  after  this  event 
would  have  a  lower  acoustic  impedance . 
Thus,  this  feature  would  appear  as  a 
buried  reflector.  An  example  of  this 
feature  is  the  opaque  reflector,  about 
100  m  below  the  surface  at  point  W  (at 
3.3  sec)  in  Figure  62,  which  continues 
into  Dragon  Seamount  to  the  right.  This 
reflector  is  more  clearly  defined  at 
0540-0600Z  and  again  at  0800Z  on  the 
seismic  record  in  Figure  63.  It  was 
difficult  to  locate  precisely  the  pro¬ 
filer  records  relative  to  the  station 
and  bomb  run.  However,  the  bottom- 
reflectivity  data  indicate  a  layer  at  50 
msec  which  could  possibly  correlate  with 
this  feature. 

DSDP  Site  136 

The  DSDP  Site  136  is  located  in  deep 
water  (4139  m)  95  km  south-southeast  of 
Station  8  beyond  the  plateau  on  which 
the  acoustic  experiment  was  conducted 
(see  Fig.  60).  The  sedimentary  cover 
over  basement  is  very  thin — only  308  m. 
In  fact,  the  short  drilling  distance 
through  the  sediments  to  the  basement 
was  one  of  the  reasons  for  choosing  the 
location  of  Site  136  [59]. 

The  shallowest  core  taken  in  the  Site 
136  drill  hole  and  the  only  one  in  the 
upper  215  m  is  in  the  130  to  139  m  in¬ 
terval.  In  this  core  no  evidence  of  tur¬ 
bidites  can  be  found.  Rjwever,  recent 
turbidites  at  the  site  cannot  be  ruled 
out  and  might  have  been  sampled  if  the 
sampling  density  in  the  drill  hole  were 
greater.  A  20  cm  thick  turbidite  with 
low  carbonate  content  was  found  4  m  down 
in  core  V27-162  taken  about  55  km  to  the 
west.  On  the  V-17  seismic-profiler  rec¬ 
ord  in  the  site  report,  there  appears 
to  be  no  topographic  barrier  which  would 
prevent  turbidite  layers  present  at  the 
site  of  core  V27-162  from  extending  to 
Site  136.  Whatever  turbidite  layers 
exist  at  Site  136  contribute  little  to 
the  sediment  thickness  and  consequently 
have  little  effect  on  the  acoustic  prop¬ 
erties  . 
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Despite  the  large  differences  in  water 
depth  and  sedimentary  thickness  at  Site 
136  and  Station  8,  the  composition  of 
the  sediments  at  the  two  locations  is 
probably  very  similar.  Strong  reflect¬ 
ing  horizons  in  the  WILKES  profiler 
record  in  Figure  62  may  have  been  creat¬ 
ed  by  changes  in  pelagic  sedimentation 
rates.  Note  that  these  horizons  can  be 
traced  from  the  base  of  the  plateau, 
past  the  obtrusion  at  Wf  and  up  the 
slope  of  the  plateau,  and  at  least  the 
uppermost  horizon  is  again  visible  ap¬ 
proaching  W* .  This  continuity  suggests 
that  the  same  deposition  mechanism  oc¬ 
curred  over  this  area. 

A  summary  of  the  sedimentary  descrip¬ 
tion  for  the  DSDP  Site  is  shown  on  the 
left  of  Figure  72.  The  upper  240  m  of 
sediments  at  the  site  have  moderate  to 
high  carbonate  content,  not  unlike  those 
samples  obtained  on  the  plateau  near 
Station  8.  Variations  in  sediment  thick¬ 
ness  occur  along  the  base  of  the  pla¬ 
teau.  Hayes  et  al.  [59]  suggest  that  the 
thin  sedimentary  cover  is  due  to  erosion 
and  nondeposition  of  sediment  (major 
hiatuses  are  represented  in  the  litholo¬ 
gic  sections).  Nevertheless,  the  evi¬ 
dence  points  to  similar  sediment  com¬ 
position  in  the  upper  200  m  at  Site  136 
and  Station  8. 

Sediment  Surface  Sound-Speed  and  Density 

Sound-speeds  were  measured  on  cores  SB-3 
and  SB-4.  The  sound-speed  profiles  shown 
in  Figure  64  suggest  a  large  gradient 
with  depth,  with  only  the  upper  few  cen¬ 
timeters  of  core  having  a  speed  less 
than  the  sound-speed  of  the  bottom 
water.  However,  SB-4  (taken  at  Station 
8)  is  too  short  to  safely  extrapolate  a 
gradient  downward  into  the  ocean  floor. 
To  estimate  the  average  sound-speed  of 
the  upper  5  m  of  sediments,  th‘e  assump¬ 
tion  was  made  that  the  lower  4  m  had  the 
same  sound-speed  as  the  average  speed  of 
the  lower  30  cm  of  the  core.  Using  these 
values,  the  estimated  adjusted  sound- 
speed  for  the  surface  sediments  is  1540 
m/sec,  with  a  sediment-to-bottom  water 
sound-speed  ratio  of  1.013. 


The  measured  values  of  wet  density  in 
the  four  samples  were  all  very  close  to 

1.60  g/cm^. 

Sound-Speed  Profile 

Based  on  the  similarity  in  sediments  at 
DSDP  Site  136  and  Station  8  as  discus¬ 
sed  in  the  geological  section,  the 
sound-speed  profile  at  Station  8  can  be 
represented  by  the  estimated  profile  at 
Site  136.  Due  to  the  differences  in 
water  depth  and  temperature,  the  sound- 
speeds  will  be  offset,  especially  near 
the  surface,  but  the  gradient  will  not 
change  significantly. 

The  sound  speed  determined  for  Site  136 
is  a  constant-gradient  model  with  g  = 
2.33  sec”l.  The  data  are  limited  and 
certain  assumptions  must  be  made  (see 

Appendix  E.2).  Four  other  methods  for 

estimating  the  sound-speed  profile  are 

presented  here  and  are  compared  to  sup¬ 
port  the  validity  of  using  the  g  =  2.33 
sec~l  profile.  The  curves  correspond¬ 

ing  to  all  the  sound-speed  prediction 
methods  are  shown  in  Figure  65.  These 
curves  are  based  on  in  situ  measurements 
and  laboratory  experiments. 

Hamilton  et  al .  [42]  have  recommended 

using  Equation  D  from  Houtz,  Ewing,  and 
Buhl  [60]  for  estimating  the  sound-speed 
profile  in  pelagic  sediments  wherever  in 
situ  data  are  not  available.  The  aver¬ 
age  sound  speed  for  this  profile  is  1789 
m/sec,  with  a  standard  deviation  of  160 
m/sec.  In  comparison,  the  average  speed 
at  Site  136  is  1867  m/sec,  or  only  78 
m/sec  higher.  Therefore,  Site  136  is  not 
anomalously  high  when  compared  to  some 
Pacific  areas  which  have  high-carbonate 
sediments.  Details  are  given  in  Appendix 
E.4. 

Houtz  and  Ewing  [39]  derived  an  equation 
from  data  covering  a  range  of  sediments 
which  includes  the  gradient  calculated 
for  Site  136.  The  details  used  in  devel¬ 
oping  the  curve  in  Figure  65  are  given 
in  Appendix  E.3. 

Laughton  [40]  made  laboratory  measure¬ 
ments  of  sound-speed  versus  overburden 
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pressure.  One  of  his  samples  taken  from 
the  North  Atlantic  in  a  water  depth  of 
2577  m  had  a  high  carbonate  content 
(54%).  His  values,  modified  in  the  upper 
135  m,  are  presented  in  Figure  65.  Fur¬ 
ther  details  and  discussion  are  in  Ap¬ 
pendix  E.l . 

The  lowest  sound-speed  profile  in  Figure 
65  is  based  on  surface  sediment  data 
from  Hamilton’s  [52,  82]  work  as  des¬ 

cribed  in  Appendix  E.5.  The  results  are 
extrapolated  to  buried  compacted  sedi¬ 
ments  by  making  adjustments  to  in  situ 
conditions  based  on  pressure  and  temper¬ 
ature  corrections.  Because  the  data  are 
from  surficial  sediment  samples,  the 
profile  is  suspect,  as  explained  in  the 
appendix.  The  profile  is  presented  with 
this  caveat  in  mind.  The  details  of  the 
adjustments  and  relationships  used  are 
also  discussed  in  the  Appendix. 

As  seen  in  Figure  65,  and  confirmed  by 
travel-time  calculations,  the  four  pro¬ 
files  derived  from  published  data  do  not 
have  sufficiently  high  sound-speeds  to 
account  for  the  time  observed  on  the 
Site  136  profiler  record.  However,  as 
discussed  above,  the  comparatively  high 
sound  speed  values  of  the  g=2.33  sec"1 
curve  are  not  unreasonable  in  light  of 
the  published  data.  Difficulty  in  relat¬ 
ing  laboratory  measurements  to  in  situ 
conditions  occurs  because  of  the  lack  of 
understanding  of  how  certain  processes 
affect  acoustic  characteristics  of  sedi¬ 
ments.  Therefore,  at  present,  techniques 
to  determine  in  situ  sound  speed  offer 
the  best  approach.  Insufficient  infor¬ 
mation  is  available  at  Site  136  to  mod¬ 
ify  the  constant-gradient  profile  deri¬ 
ved  using  the  CRIC  method.  Consistent 
with  the  general  trend  of  a  decreasing 
gradient  with  depth  observed  in  the 
other  profiles,  the  gradient  value  of 


2.33  sec~^  would  be  a  lower  limit  in 
the  uppermost  sediments  and  an  upper 
limit  in  the  deeper  sediments.  However, 
the  average  gradient  is  the  best  esti¬ 
mate  when  relating  Site  136  to  the  upper 
200  m  of  sediments  at  Station  8. 

Attenuation 

Grain  size  determined  on  samples  from 
both  cores  taken  at  Stations  7  and  8  are 
shown  in  Table  13.  Only  one  of  four  sam¬ 
ples  contained  a  substantial  percentage 
of  sand-sized  particles.  The  values  of 
attenuation  are  determined  directly  from 
Hamilton’s  [45]  graph  and  a  simple  aver¬ 
age  is  computed  to  determine  the  attenu¬ 
ation  constant  of  < =  0.25  dB/m/kHz.  The 
small  sampling  may  not  be  adequate,  but 
with  shallow  carbonate  deposits,  large 
grain  size  is  expected  with  a  corre¬ 
sponding  higher  value  of  k  .  Thus,  the 
value  of  k  is  the  minimum  that  would  be 
calculated  by  this  method  even  with  a 
greater  sampling  of  size  distribution. 

Station  9 

Station  9  is  located  on  the  lower  conti¬ 
nental  rise  west  of  North  Africa  on  the 
eastern  edge  of  the  Cape  Verde-Madeira 
Abyssal  Plain  in  a  water  depth  of  about 
4600  m.  The  almost  featureless  topog¬ 
raphy  provides  a  site  conducive  to 
modeling.  The  bomb  run  for  Station  9, 
shown  in  Figure  66,  commences  at  29°N, 
20°50’W.  The  location  provided  the 
opportunity  to  conduct  a  long-range  pro¬ 
pagation  experiment  as  indicated  in  the 
figure  on  the  track  to  the  northwest. 

Geological  Information 

Several  cores  have  been  taken  in  this 
area  by  the  Lamont -Doherty  Geological 
Observatory.  These  samples  are  indicated 


Table  13.  Sound  attenuation  from  cores  at  stations  7  and  8 


Sampling 

Mean  Grain 

kG 

Core 

Interval 

Sand-Silt-Clay 

Size  i 

db/m/kHz 

SB-3 

12-18  cm 

27-22-51% 

7.1 

0.08 

Station  7 

22-28  cm 

33-25-43% 

6.6 

0.10 

SB-4 

22-28  cm 

77-10-12% 

4.2 

0.72 

Station  8 

64-72  cm 

26-22-53% 

7.2 

0.09 
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in  Figure  66.  For  the  most  part,  cores 
in  the  area  have  carbonate  sediments. 
Within  these  cores  are  layers  charac¬ 
teristic  of  turbidites,  but  of  carbon¬ 
aceous  rather  than  mineral  composition. 
However,  mineral  (terrigenous)  turbi¬ 
dites  are  observed  in  many  cores  as  in¬ 
dicated  in  Figure  66  by  the  solid  sym¬ 
bol.  Because  Africa  is  so  highly  arid 
and  no  major  rivers  empty  into  the 
Atlantic  that  far  to  the  north,  the 
major  source  of  turbidites  at  the  sta¬ 
tion  is  probably  the  Canary  Islands, 
with  a  possible  contribution  from  the 
Madeira  Islands.  Windblown  quartz  and 
clay  minerals  from  Africa  are  also  an 
important  contribution  to  the  sediments. 

Sedimentation  in  this  region  is  typical 
of  abyssal  plains,  i.e.,  alternation  of 
pelagics  and  turbidites,  but  topograph¬ 
ically  it  lacks  the  flat  character  of 
abyssal  plains.  The  seismic  data,  and  on 
a  more  detailed  scale,  the  3.5  kHz  data, 
show  the  surface  to  consist  of  irregular 
depth  variations  of  5  m  amplitude.  Re¬ 
gionally,  the  sea  floor  surface  rises 
toward  the  east  with  a  prominent  change 
in  topography  occurring  near  the  end  of 
the  long-range  acoustic  run.  The  records 
are  shown  in  Figure  67. 

The  track  approaching  Station  9  from 
Station  8  and  the  track  just  beyond  the 
end  of  the  bomb  run  traverse  a  small 
double  channel  that  is  oriented  NNW. 
Another  channel  with  the  same  orienta¬ 
tion  is  observed  on  both  tracks;  the 
crossing  on  the  propagation  run  track  is 
at  19°07,W.  These  channels  are  indica¬ 
ted  by  dashed  lines  in  Figure  66.  All 
the  crossings  show  a  channel  relief  of 
5-10  m  except  the  latter  crossing  which 
has  a  relief  of  20  m.  The  channels  are 
about  0.3  km  wide  on  the  southern  KANE 
tr^ck  crossing,  but  are  not  observed  on 
the  WILKES  southern  track,  which  is  far¬ 
ther  south  and  up-slope  from  *the  KANE 
track.  Embley  [61]  has  found  evidence  of 
debris  flow  deposits  just  southeast  of 
the  bomb  run  for  this  station.  He  has 
identified  this  type  of  deposit  on  the 
3.5  kHz  depth  recorder  as  a  thin  lens 
lying  on  the  ocean  floor.  This  charac¬ 
teristic  signature  does  not  appear  on 


the  records  in  Figure  67.  However,  the 
NNW  flow  direction  indicated  by  the 
channels  coincides  with  the  flow  direc¬ 
tion  indicated  by  Embley. 

The  Canary  Islands  form  a  barrier  along 
the  coast  of  Africa,  restricting  the 
flow  of  turbidites  to  a  path  south  of 
the  islands  in  the  vicinity  of  Station 
9.  Uchipi  et  al.  [58]  suggest  that  the 
whole  region  along  the  long-range  prop¬ 
agation  run  obtained  its  major  constit¬ 
uents  of  sediments  prior  to  the  forma¬ 
tion  of  the  Canary  Islands  by  volcanic 
activity  in  the  Oligocene. 

Sediment  Surface  Sound-Speed  and  Density 

An  estimate  of  the  surficial  sediment 
sound-speed  is  based  on  the  sediment 
type  as  estimated  from  cores  in  the 
area.  The  cores  have  a  higher  carbonate 
content  than  those  at  the  other  deep 
stations,  but  less  than  those  at  Station 
8.  The  contribution  of  wind-blown  clay 
particles  is  greater  than  at  Station  8 
farther  north,  but  the  clay  content  is 
less  than  that  received  from  turbidity 
currents  in  the  Tagus  and  Iberian  Abyss¬ 
al  Plains.  No  sound-speed  measurements 
were  taken  on  cores  from  the  area.  By 
estimating  an  average  type  clay  size 
particle  distribution  with  a  porosity  of 
70%,  statistical  relationships,  such  as 
in  Anderson  [62],  can  be  used  to  predict 
average  acoustic  parameters.  In  this 
relationship,  the  carbonate  content  is 
not  a  factor.  The  predicted  sediment- 
water  sound-speed  ratio  and  density  for 
porosity  of  70%  are  0.984  and  1.52 
gm/ crn^  9  respectively . 

Sound-Speed  Profile 

Extensive  sonobuoy  data  measured  by 
Lamont-Doherty  personnel  in  this  area 
were  made  available  by  Robert  Houtz;  the 
locations  of  most  are  shown  in  Figure 
66.  The  high  density  of  stations  (al¬ 
though  no  stations  are  on  the  acoustic 
measurement  track)  provides  a  method  of 
estimating  a  sound-speed  profile.  This 
method  entails  the  use  of  a  statistical 
compilation  of  the  sonobuoy  data  to 
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derive  a  gradient  in  the  upper  200  m. 
The  data  and  details  of  the  method  and 
assumptions  are  in  Appendix  E.4.  The 
estimate  of  the  gradient  is  1.1  sec-^- 
with  a  standard  deviation  of  about  one- 
half  this  value. 

Attenuation 

In  order  to  estimate  the  attenuation  at 
Station  9,  the  mean  grain  size  distri¬ 
bution  must  first  be  estimated.  The  core 
data  suggest  a  greater  preponderance  of 
calcareous  sediments  than  are  found  in 
the  Tagus  Abyssal  Plain  at  Stations  4 
and  5.  Also,  the  sound-speed  gradient  at 
Station  9  is  greater  than  at  Stations  4 
and  5.  These  two  observations  are  self- 
consistent  and  further  suggest  a  larger 
mean  grain  size  at  Station  9  than  at  the 
other  abyssal  plain  stations,  but  not  as 
high  as  observed  at  Station  8.  Thus, 
without  actual  grain -size  measurements, 
an  estimated  value  of  about  0.10  dB/m/ 
kHz  based  on  Hamilton  [45]  is  proposed. 
Although  very  subjective,  this  value  is 
consistent  with  an  apparent  relationship 
between  estimates  of  the  gradient  and 
attenuation  constant  at  the  other  sta¬ 
tions  . 

Station  7 

The  site  for  Station  7  shown  in  Figure 
60  provides  a  rough  sea  floor  with  es¬ 
sentially  no  sediments.  Because  of  the 
high  acoustic  impedance  of  the  ocean 
bottom,  sound  propagation  paths  in  the 
bottom  are  not  a  consideration  in  ana¬ 
lyzing  the  data.  On  the  other  hand, 
reverberation  is  very  prominent.  This 
type  of  bottom  is  acoustically  repre¬ 
sentative  of  a  significant  portion  of 
the  ocean  floor,  which  consists  of  the 
crest  of  the  Mid-Oceanic  Ridge  or  any 
other  center  of  sea-floor  spreading.  The 
particular  site  was  chosen  because  its 
location  was  within  a  reasonable  vicin¬ 
ity  of  the  cruise  track  plan. 

Geological  Information 

The  elongated  basaltic  outcrop  which  the 
bomb  run  traverses  is  called  the  Dragon 
Seamount.  It  rises  to  a  minimum  depth  of 


1200  m.  The  most  detailed,  published 
geophysical  description  of  this  site  is 
reported  by  Laughton  et  al .  [63].  The 

bathymetry  of  the  Stations  7  and  8  vi¬ 
cinity,  as  illustrated  in  Figure  60,  is 
primarily  based  on  NAV0CEAN0  tracks, 
including  surveys  on  the  October  1973 
cruises  and  other  tracks  and  contours  by 
Hayes  et  al .  [59]  as  shown  in  the  DSDP 

Site  136  report. 


Because  of  the  roughness  of  the  track, 
3.5  kHz  records  were  not  well-defined. 
The  narrow  beam  12  kHz  record  is  shown 
in  Figure  63  at  one-half  the  scale  used 
for  the  3.5  kHz  records  of  other  sta¬ 
tions.  Along  the  bottom  of  the  figure  is 
a  line  drawing  of  the  bathymetry  used  in 
the  ray  trace  program  with  no  vertical 
exaggeration.  This  line  drawing  is  to 
assist  the  reader  in  visualizing  the 
effects  of  the  rough  topography  on  ray 
paths  with  the  true  angles  presented. 
The  seismic  profile  section  shown  in  the 
insert  in  Figure  63  is  one  crossing  of 
the  bomb  run  made  on  the  survey  by  the 
WILKES. 

Sediment  Sound-Speed  and  Density 


Core  SB-3  taken  in  the  outcrop  area  of 
Station  7  struck  hard  rock  after  pene¬ 
trating  less  than  a  meter  of  sediment. 
The  sediment  consists  of  predominantly 
foraminiferal  sandy  clay.  Sound-speed 
measurements  were  taken  on  this  core 
aboard  ship  and  are  shown  in  Figure  64. 
The  average  adjusted  sound-speed  in  the 
core  is  1530  m/sec,  giving  a  sediment - 
to-bottom  water  ratio  of  1.005.  Where 
this  sedimentary  cover  does  not  exist, 
however,  the  igneous  rock  is  in  direct 
contact  with  the  bottom  water.  Photo¬ 
graphs  taken  by  Laughton  et  al.  [63] 
show  steep  slopes  void  of  sediment.  Fox 
et  al.  [64]  measured  the  compressional 
sound-speeds  of  basalt  samples  taken 
from  DSDP  Site  136,  80  km  south  of  the 
station.  Measurements  at  0.25  kb  (equiv¬ 
alent  to  about  2500  m  water  depth)  were 
3880  m/sec,  and  the  density  was  2.5 
g/cm^.  However,  this  sample  is  highly 
altered  basalt,  and  Fox  et  al.  [65]  sug¬ 
gest  that  these  values  are  considerably 
lower  than  relatively  fresh  and 
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unaltered  basalt  measured  by  others.  For 
example,  Christensen  [66]  measured  5100 
m/sec  at  0.20  kb  on  a  water-saturated 
basalt  sample  with  a  density  of  2.90 
g/cm^  from  the  Juan  de  Fuca  Ridge. 
This  measurement  corresponds  well  with 
the  average  sound  speed  for  the  oceanic 
crust  (Layer  2).  Laughton  et  al.  [63] 
report  olivine  basalt  as  the  main  rock 
dredged  on  Dragon  Seamount.  Using  either 
density-sound  speed  set  gives  a  very 
high  reflection  coefficient  at  the 
water-basalt  interface  of  0.7-0. 8  (nor¬ 
mal  incidence  bottom  loss  of  1.8-2. 8 
dB).  The  sediment  covering  is  too  thin 
to  significantly  affect  the  reflectivity 
at  low  frequencies  over  the  bomb  run 
track.  A  possible  model  may  be  a  meter 
of  sediment  covering  a  few  meters  of 
weathered  basalt  over  a  core  of  unalter¬ 
ed  basalt. 
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Appendix  D. 

Geological  Aspects  of  Acoustic  Parameters 


Seafloor  sediments  contain  features  that 
create  reflecting  horizons,  preferential 
acoustic  propagation  paths,  and  other 
spatial  variations  of  acoustic  param¬ 
eters,  all  of  which  constitute  the 
acoustical  stratigraphy  of  the  ocean 
floor.  The  ability  to  extrapolate  the 
acoustical  stratigraphy  beyond  areas  of 
direct  measurements  requires  understand¬ 
ing  of  the  physical  structure  which 
influences  it.  We  are  still  far  from  a 
thorough  understanding  because  the  cor¬ 
relation  of  the  geological  stratigraphy 
and  acoustical  stratigraphy  is  not  sim¬ 
ple.  The  objectives  of  this  section  are 
to:  (1)  emphasize  sedimentary  processes 
relevant  to  the  acoustical  stratigraphy 
at  each  station,  (2)  point  out  problems 
in  the  correlation  of  geological  and 
acoustical  stratigraphy,  and  (3)  de¬ 
scribe  the  acoustic  parameters  empha¬ 
sized  in  this  report  and  their  methods 
of  measurement. 

Sedimentary  Processes 

A  few  predominant  sedimentary  proces¬ 
ses  have  produced  the  primary  types  of 
sedimentary  layers  which  affect  the 
acoustic  stratigraphy  at  the  station 
locations  in  this  report.  Sediments 
emplaced  by  other  processes  may  exist, 
but  not  to  a  significant  degree.  One 
such  process  is  the  precipitation  of 
dissolved  minerals  in  the  water  onto 
particles  on  the  sea  floor.  An  example 
is  the  formation  of  manganese  nodules, 
which  can  be  extensive  enough  to  form 
manganese  pavements.  The  more  apparent 
sedimentary  processes  are  discussed  in 
the  following  paragraphs. 

Pelagic  Sediments 

Pelagic  sediments  are  formed  by  the  dep¬ 
osition  of  the  shells  (tests)  of  fauna 
and  flora  which  float  in  the  ocean 
(plankta).  The  predominant  test  material 


of  plankta  is  calcium  carbonate  or  cal- 
cite.  Coarse  calcareous  sediments  con¬ 
sist  predominantly  of  the  tests  of  fo- 
raminifera,  whereas  the  clay-size  par¬ 
ticles  are  tiny  calcareous  plates  called 
coccoliths.  Under  proper  environmental 
conditions,  the  calcareous  sediments 
will  partially  or  totally  solidify  to 
form  chalk  or  limestone.  Variation  in 
these  conditions  and  in  sedimentation 
rates  results  in  variations  in  the 
acoustic  impedance  of  the  pelagic  sedi¬ 
ments.  In  general,  surface  sediments 
have  sound  speeds  less  than  the  bottom 
water,  with  the  sound  speed  increasing 
with  depth  of  burial.  Hamilton  et  al . 
[42]  have  observed  that  the  sound  speed 
increases  more  rapidly  with  depth  of 
burial  in  sediments  with  high  carbonate 
content  than  sediments  with  predominant¬ 
ly  terrigenous  clay.  The  sediments  at 
Station  8  have  high  carbonate  content 
with  an  implied  higher  sound-speed  gra¬ 
dient  than  at  the  other  stations.  In 
water  depths  greater  than  4000  m  the 
carbonates  tend  to  dissolve,  inhibiting 
their  contribution  to  the  sediments.  The 
degree  of  dissolution  is  dependent  cn 
the  hydrostatic  pressure  and  other  pa¬ 
rameters  in  the  water  which  have  fluctu¬ 
ated  throughout  geologic  history,  as 
well  as  the  amount  of  carbonates  avail¬ 
able.  (The  reader  is  referred  to  Berger 
[67]  for  a  discussion  of  the  calcite 
compensation  depth.) 

In  addition  to  the  calcareous  sources, 
silicof lagellates  and  diatoms,  both 
algae,  contribute  silica  constituents  to 
the  sediment.  Since  silicates  do  not 
dissolve  in  deep  waters  as  readily  as 
carbonates,  the  ratio  of  silica  to  cal¬ 
cite  increases  with  water  depth.  Where 
siliceous  sediments  are  dominant,  they 
may  change  form  from  siliceous  oozes  to 
the  denser  form  of  chert  with  a  much 
higher  sound  speed  (about  3.6  km/sec). 
The  formation  of  chert  layers  is  not 
understood,  although  explanation  relat¬ 
ing  their  formation  to  volcanism  have 
been  suggested  [67,  68].  Chert  layers 

have  been  found  in  the  DSDP  sites,  gen¬ 
erally  buried  deeper  than  250  m  below 
the  sea  floor,  although  they  have  been 
cored  as  shallow  as  100  m. 
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Generally,  the  layers  have  been  deter¬ 
mined  to  be  older  than  Miocene  (26 
mybp) •  These  layers  are  thin  (about  5 
cm),  but  several  layers  may  occur  to¬ 
gether,  thus  creating  a  good  reflector 
of  sound  over  a  large  frequency  range. 
Horizon  A,  an  extensive,  strong  sub¬ 
bottom  reflector  observed  in  the  north¬ 
western  Atlantic  and  correlative  to 
horizons  in  the  eastern  Atlantic  [58], 
has  been  drilled  and  correlates  with 
several  thin  chert  layers  [69,  70].  The 
age  of  this  layer  is  Eocene  (50  mybp). 
Horizon  A  may  therefore  exist  beneath 
the  deep  abyssal  plains,  but  is  masked 
by  more  recent  sediments.  The  evidence 
indicates  that  chert  layers  are  not  nec¬ 
essarily  time-synchronous,  and  that  they 
may  be  more  common  than  previously 
thought.  The  existence  of  chert  as  a 
reflector  at  any  of  the  stations  can  not 
be  ruled  out. 

Eolion  Sediments 

Eolian  sediments  are  particles  initially 
deposited  on  the  ocean  surface  by  the 
winds.  One  such  type  is  ash  that  has 
been  distributed  over  a  broad  area  by 
wind  during  volcanism.  Such  ash  layers 
may  be  only  a  few  centimeters  thick,  or 
extremely  thick  as  are  found  in  the  Med¬ 
iterranean  Sea  basins.  The  layers  are 
independent  of  water  depth  and  type  of 
surrounding  sediment  and  create  good 
sound  reflectors.  Several  thin  ash  lay¬ 
ers  which  are  close  together  may  form 
one  reflecting  horizon  for  a  range  of 
wavelengths  of  the  impinging  sound  ener¬ 
gy.  Ash  layers  are  the  most  likely  cause 
of  some  of  the  reflectors  seen  near  sev¬ 
eral  of  the  stations.  In  particular,  a 
rather  shallow  reflector  at  Station  8  is 
likely  to  be  an  ash  layer. 

Aside  from  ash,  other  land-derived  de¬ 
bris  is  carried  by  the  winds  far  out  to 
sea  and  contributes  to  the  sediments. 
This  debris,  even  up  to  the  silt-size, 
may  be  a  significant  constituent  of  the 
sediments  west  of  the  Sahara  Desert, 
such  as  at  Station  9.  Windblown,  clay¬ 
sized  particles  contribute  in  some  de¬ 
gree  to  almost  all  the  sediments  in  the 
world’s  oceans. 


Gas  Hydrates 

A  gas  hydrate  is  a  crystalline  struc¬ 
ture  of  water  molecules  in  which  gas 
molecules  are  physically  trapped  [71]. 
This  formation  produces  a  substantial 
increase  in  sound  speed  which,  in  turn, 
produces  a  strong  reflecting  horizon. 
Conditions  under  which  this  phenomenon 
occurs  exist  in  the  near-surface  sedi¬ 
ments;  their  geographical  distribution 
is  unknown  although  their  presence  at 
several  locations  has  been  verified. 
These  hydrates  leave  no  measurable  evi¬ 
dence  on  bottom  samples  other  than  the 
expansion  observed  before  a  core  sample 
is  extruded  from  the  coring  tube.  Be¬ 
cause  of  the  suspected  widespread  dis¬ 
tribution,  this  feature  may  account  for 
some  of  the  reflectors  observed  at  the 
stations . 

Turbid  ites 

Underwater  flow  mechanisms  that  con¬ 
tribute  to  the  redistribution  of  sedi¬ 
ments  or  emplacement  of  land-derived 
(terrigenous)  deposits  are:  bottom  con¬ 
tour  currents,  the  nepheloid  layer,  de¬ 
bris  flows  (observed  near  Station  9), 
and  turbidity  currents.  The  last  mech¬ 
anism  is  the  most  obvious  contributor  of 
sediments  at  the  stations. 

The  deposition  of  turbidite  layers  ac¬ 
counts  for  the  thick  sediments  in  the 
deep  ocean  basins’  abyssal  plains.  These 
deposits  are  generally  land-derived  min¬ 
eral  material  (predominantly  quartz) 
mixed  with  shallow-water  marine  deposits 
of  high  carbonate  content. 

Turbidite  currents  are  not  restricted  to 
the  continental  shelf,  but  can  occur 
along  any  slope  (such  as  that  off  an  is¬ 
land).  Material  from  rivers  and  the  con¬ 
tinental  shelf  piles  up  on  the  shelf 
until  the  deposit  becomes  unstable.  A 
microearthquake  can  disturb  the  pile  and 
trigger  a  down-slope  flow  of  a  mixture 
of  the  material  and  water.  This  turbid 
water  builds  up  speed  under  the  influ¬ 
ence  of  gravity  due  to  its  excess  den¬ 
sity  and  low  friction,  and  carries  the 
shelf  material  down  submarine  canyons 
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along  the  continental  slope  into  the 
deep  ocean  basins.  Here,  the  material  is 
deposited  as  the  flow  velocity  decreases 
due  to  decrease  of  the  downhill  slope. 
The  size  of  particles  carried  by  the 
turbid  water  is  a  function  of  the  veloc¬ 
ity:  the  larger  particles  settle  out 
first,  and  the  finer  particles  later. 
The  result  is  that  during  a  turbidity 
flow,  the  first  deposited  particles  are 
the  largest.  Thus,  the  larger  particles 
are  deposited  in  greater  amounts  at 
locations  closer  to  the  source.  Turbid¬ 
ity  flows  tend  to  flatten  the  topography 
of  the  receiving  basin  thus  producing 
the  characteristic  plain.  The  flatness 
is  conducive  to  acoustic  modeling;  how¬ 
ever,  the  complex  layering  complicates 
the  models. 

In  a  sedimentary  column,  the  turbidite 
layer  represents  a  rapid  deposition 
which  is  interjected  into  the  stream  of 
continuous  deposition  of  pelagic  and 
eolian  sediments.  The  turbidite  layer 
may  be  represented  in  an  idealized  model 
by  its  grain-size  distribution.  Large 
grain-sized  material  is  found  at  the 
base  and  gradually  grades  into  smaller 
sized  material  at  the  top.  The  decrease 
in  grain-size  upward  results  from  the 
decreasing  velocity  of  the  turbidity 
flow  in  time.  The  basal  interface  is 
sharply  defined  by  the  sudden  increase 
in  grain-size  compared  to  that  of  the 
underlying  sediments.  At  the  top,  how¬ 
ever,  the  interface  with  the  overlying 
sediments  of  similar  grain-size  is  im¬ 
perceptible.  A  common  perturbation  of 
the  model  is  the  diffusing  of  the  sharp 
bottom  interface  by  animal  burrowing.  A 
more  complex  grain-size  profile  than 
described  here  is  observed  at  DSDP  Site 
118  northeast  of  Station  2  by  Laughton 
et  al.  [49].  During  periods  of  glacia¬ 
tion,  the  mean  sea  level  is  lowered  due 
to  the  removal  of  water  to  the  glaciers. 
This  causes  an  increase  in  erosion  which 
provides  the  source  of  turbidites.  The 
glacial  period  commenced  in  the  northern 
hemisphere  about  3  mybp  in  the  Pliocene. 
The  peak  of  glacial  intensity  occurred 
about  0.4  mybp  [72].  (The  Pleistocene, 
often  associated  with  the  glacial  per¬ 
iod,  began  at  1.8  mybp.)  Although  four 


or  five  major  glacial  advances  are  ob- 
served  on  the  continents,  the  deep— ocean 
paleontological  record  shows  many  fluc¬ 
tuations  in  the  glacial  cycle;  each 
would  have  caused  a  variation  in  the 
turbidite  activity.  During  periods  of 
intense  glaciation,  less  accumulation  of 
pelagic  sediments  occurs  between  turbid¬ 
ity  flows;  thus,  the  percentage  of 
large,  grain-sized  basal  portions  of 
turbidites  in  the  sedimentary  column 
represented  by  this  period  is  greater. 

Turbidity  flows  contributed  to  the  bulk 
of  the  sediments  at  Station  2  in  the 
Iberian  Abyssal  Plain  and  Stations  4  and 
5  in  the  Tagus  Abyssal  Plain.  The  verti¬ 
cal  sound-speed  profile  through  a  turbi¬ 
dite  layer  responds  to  the  variation  of 
density  and  mean  grain-size  in  the 
layer.  The  sound-speed  profile  in  the 
turbidites  at  Stations  4  and  5  is  de¬ 
scribed  in  Appendix  C.  The  number  of 
turbidite  layers  influences  the  average 
sound  speed  of  the  sediments.  The  ver¬ 
tical  spatial  frequency  can  be  inferred 
from  seismic-profiler  records  as  dis¬ 
cussed  in  the  following  section. 

Acoustic  Reflectors  and  Geological  Interfaces 

The  correlation  between  the  geological 
and  acoustical  stratigraphy  is  not  al¬ 
ways  apparent.  This  lack  of  correlation 
between  acoustic  and  geological  inter¬ 
faces  has  often  been  observed  in  DSDP 
sites.  Theoretical  reasons  for  this 
phenomenon  exist.  Even  in  the  simplest 
approach  relating  the  acoustical  stra¬ 
tigraphy  to  the  geological  stratigraphy, 
three  limits  must  be  considered.  These 
limits  involve  the  thickness  of  the 
interface  zone,  the  thickness  of  a  re¬ 
flecting  layer,  and  the  separation  of 
reflecting  layers. 

A  reflecting  interface  occurs  where  an 
impedance  contrast  exists  between  two 
adjacent  media.  Maximum  reflection  oc¬ 
curs  when  the  contrast  occurs  over  a 
zone  which  is  infinitesimal  relative  to 
the  acoustic  wavelength.  On  the  other 
hand,  no  reflection  occurs  when  the  zone 
is  infinite  in  length  compared  to  a 


125 


wave-length.  According  to  ray  theory, 
this  effect  causes  a  refracting  of  the 
ray  when  the  ray  is  oblique  (i.e.,  if  a 
speed  change  is  involved  rather  than 
just  a  density  change).  Several  condi¬ 
tions  may  exist  which  mask  potential 
reflecting  horizons.  For  example,  slow 
changes  in  the  sedimentary  deposition 
process  may  create  a  wide  interface  zone 
even  at  wavelengths  comparable  to  low 
frequency  energy;  at  high-frequency 
wavelengths  one  may  see  the  effects  of 
burrowing  by  benthic  creatures. 

At  normal  incidence,  peculiar  phenomena 
occur  on  the  reflecting  of  half-wave- 
length  and  quarter-wavelength  layers. 
Brekhovskikh  [38,  pg.  55-56]  offers  a 
concise  explanation  of  this  phenomenon 
of  a  layer  lying  between  media  1  and 
media  3  with  impedance 

"The  halfwave  layer  has  no  effect  on  the 
incident  wave  (as  if  the  layer  were  ab¬ 
sent),  and  the  reflection  coefficient  is 
equal  to  the  usual  coefficient  of  re¬ 
flection  from  the  boundary  between  media 
3  and  1  just  as  if  they  were  in  direct 
contact  with  one  another.  In  particular, 
if  media  3  and  1  have  the  same  proper¬ 
ties  (Zj  =  Z^),  the  reflection  coef¬ 

ficient  is  zero." 

This  explanation  is  also  true  for  an  in¬ 
tegral  number  of  half-wavelength  thick¬ 
nesses.  This  limit  may  apply  where  tur- 
bidite  layers  alternate  with  pelagic 
layers.  For  the  quarter-wavelength  thick 
layer  (1/4A): 

"Reflection  at  the  boundary  between  any 
two  media  can  be  eliminated  entirely  by 
inserting  a  quarterwave  layer  between 
the  two  media,  the  layer  having  an  impe¬ 
dance  equal  to  the  geometric  mean  of  the 
impedances  of  the  two  media.  A  similar 
effect  will  be  observed  when  the  layer 
thickness  is  equal  to  3/4A,  5/4 X,  etc." 

The  effect  of  the  separation  of  layers 
can  be  thought  of  in  terms  of  the  layer 
thickness  of  the  separating  medium  with 
the  implications  mentioned  in  the  pre¬ 
ceding  paragraph.  Another  problem  occurs 
with  resolution  if  the  spacing  is  so 


small  that  the  returned  pulse  trains 
overlap.  The  chert  or  sandstone  layers 
that  would  individually  be  transported 
to  long  acoustic  wavelengths  may  appear 
as  a  single  reflector  if  several  are 
grouped  together. 

With  additional  mathematical  terms,  the 
above  relationship  can  also  be  applied 
to  oblique  angles.  These  relationships 
are  limits,  rather  than  definite  values 
for  layer  thicknesses,  implying  grada¬ 
tional  effects. 

Two  kinds  of  seismic  systems  for  record¬ 
ing  travel  time  to  subbottom  layers  were 
used  consistently  on  the  cruise:  3.5  kHz 
and  the  sparker  (50-150  Hz).  These  sys¬ 
tems  display  signals  reflected  at  normal 
incidence  from  interfaces  of  contrasting 
acoustic  impedance  below  the  ocean  floor 
as  well  as  the  water-sediment  interface. 
The  3.5  kHz  recorder  is  useful  in  defin¬ 
ing  reflectors  in  the  upper  50  m  of 
sediment.  The  sparker  seismic  records, 
on  the  other  hand,  are  responsive  to 
interfaces  which  affect  low-frequency 
sound  propagation  and  thus  are  closely 
related  to  the  bottom-reflectivity  meas¬ 
urements.  Although  the  near-surface  de¬ 
tail  is  lost  due  to  the  lower  frequency 
and  it  has  a  more  compressed  time  scale 
than  the  3.5  kHz  recorder,  penetration 
is  much  greater — reaching  to  two  seconds 
of  two-way  travel  time,  or  about  2  km  of 
sediment  in  the  areas  reported  here. 
Since  the  sparker  signal  would  interfere 
with  the  bottom-reflectivity  experiment, 
it  could  not  be  operated  during  the  ex¬ 
periment.  The  3.5  kHz  system,  on  the 
other  hand,  operating  at  a  much  higher 
frequency  and  at  lower  power  with  down¬ 
ward  directivity,  could  be  operated  dur¬ 
ing  the  bottom-reflectivity  run  without 
interfering  with  the  experiment. 

Most  of  the  reflectors  seen  on  the  seis¬ 
mic  records  are  within  the  sedimentary 
column.  Where  penetration  with  the 
sparker  system  is  sufficient,  a  strong 
reflector  called  basement  is  observed 
which  generally  correlates  with  oceanic 
layer  "2".  Basement  has  a  compressional 
wave  speed  of  at  least  3  km/ sec  and  con¬ 
sists  of  a  high  percentage  of  basalt 
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[73].  This  layer  was  generated  at  the 
source  of  sea  floor  spreading,  and,  not 
being  subjected  to  the  erosional  forces 
that  exist  on  land,  has  roughness  equiv¬ 
alent  to  that  of  the  younger  rough  areas 
of  the  ocean  floor  such  as  the  Mid- 
Atlantic  Ridge.  An  example  of  this  rough 
layer  can  be  seen  on  some  of  the  seis¬ 
mic-profile  records  in  this  report.  The 
acoustic  basement  is  defined  as  the  bot¬ 
tom-most  layer  possible  to  observe  on 
the  profiler  records.  Acoustic  basement 
may  not  be  the  same  as  basement  in  the 
geologic  sense,  but  instead,  may  be  the 
acoustic  signature  of  a  layer  buried 
under  unconsolidated  sediments.  In  this 
case,  a  small  amount  of  energy  is  trans¬ 
mitted  to  reflect  off  deeper  interfaces 
which  would  have  less  impedance  con¬ 
trast;  thus,  the  reflected  energy  would 
be  imperceptible. 

Subbottom  depth  recordings  in  abyssal 
plains  characteristically  display  sev¬ 
eral  reflecting  horizons  above  the 
acoustic  basement.  One  parameter  which 
describes  an  aspect  of  the  acoustical 
behavior  of  geological  features  in  the 
sediments  is  the  number  of  reflectors 
per  meter.  This  parameter  has  been  used 
in  Table  3  to  describe  the  3.5  kHz  sub¬ 
bottom  records  at  each  station. 

The  number  of  reflecting  horizons  ob¬ 
served  in  3.5  kHz  records  and  in  the 
lower  frequency  seismic-profiler  rec¬ 
ords,  and  the  number  of  reflecting  in¬ 
terfaces  observed  in  deep-sea  cores 
differ.  In  a  study  of  the  sediments  at 
Stations  4  and  5  in  the  Tagus  Abyssal 
Plain,  Ruddiman  and  Glover  [53]  made  a 
comparison  of  the  number  of  reflectors 
per  meter  in  the  records  and  the  cores. 
For  resolution  considerations,  in  com¬ 
paring  reflectors  on  the  3.5  kHz  records 
to  reflective  layers  in  the  cores,  two 
adjustments  must  be  made  in  the  core 
data:  first,  turbidites  spaced  less  than 
a  quarter-wavelength  (12  cm)  apart  were 
combined  into  a  single  reflector;  sec¬ 
ond,  layers  appreciably  less  than  a 
quarter  wavelength  thick  were  ignored. 
An  average  of  roughly  0.75  reflectors/ 
meter  in  the  cores  taken  at  acoustic 
Station  4  and  5  was  calculated.  In  con¬ 


trast,  the  underway  3.5  kHz  data  con¬ 
sistently  recorded  0.25-0.40  reflectors 
per  meter.  Some  of  the  loss  may,  of 
course,  be  cause  by  ship’s  motion  and 
other  masking  effects  due  to  noise. 
Thus,  the  near-surface  Tagus  Abyssal 
Plain  sediments  probably  consist  of  two 
to  three  times  as  many  individual  major 
turbidites  than  are  generally  visible  in 
the  3.5  kHz  data. 

To  examine  the  resolution  in  seismic- 
profiler  records,  Ruddiman  and  Glover 
[53]  counted  the  number  of  reflectors 
per  depth  of  sediment  on  the  records. 
Averages  of  0.065-0.085  reflectors  per 
meter  were  calculated  in  the  upper  0.5- 
1.0  seconds  (two-way  travel  time),  or 
about  400-800  m  on  the  profiles.  This 
figure  is  roughly  an  order  of  magnitude 
fewer  than  the  number  of  well-defined, 
thick  turbidites  observed  in  the  cores, 
and  a  factor  of  four  fewer  reflectors 
per  meter  than  are  visible  in  the  3.5 
kHz  data.  Each  individual  reflector 
observed  on  the  sparker  records  must, 
therefore,  be  caused  by  groups  of  adja¬ 
cent  turbidites.  Using  the  rough  esti¬ 
mate  of  sparker  resolution,  as  many  as 
10  major  turbidite  layers  (and  5-10  more 
minor  layers)  may  contribute  to  one  in¬ 
dividual  seismic  reflector. 

Geoacoustic  Parameters 

Sediment  Surface  Density  and  Sound-Speed 

The  detection  of  horizons  on  the  seismic 
recorders  by  normally  incident  acoustic 
waves  is  a  function  of  the  contrast  in 
acoustic  impedance  of  the  medium  at  an 
interface.  The  impedance  is  a  product  of 
two  measurable  properties:  density  and 
compressional  wave  (sound)  speed.  The 
density  is  a  commonly  measured  property 
on  core  samples  and  is  the  wet  weight 
divided  by  the  volume.  Most  near-surface 
sediments  consist  of  clay-sized  parti¬ 
cles  and  have  densities  of  about  1.5 
g/cm^  but  can  range  from  1.2  to  1.6 
g/cm^.  Higher  densities  are  generally 
associated  with  silts  or  sands  reaching 
an  extreme  of  2.6  g/cm^  and  above  for 
nonporous  material. 
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Fry  and  Raitt  [74]  described  a  method 
(using  a  receiver  and  source  at  the  sea 
surface)  of  determining  the  sound-speed 
of  ocean  floor  sediments  where  the 
speeds  are  less  than  that  of  the  bottom 
water.  Since  then,  new  evidence  indi¬ 
cates  that  low-speed  sediments  at  the 
water-sediment  interface  are  the  rule 
rather  than  the  exception  in  the  deep 
ocean.  Urick  [75]  considers  a  mixture  of 
dispersed  clay  particles  in  sea  water 
and  describes  the  theoretical  basis  for 
a  low— speed,  clayey  sediment.  Because 
pelagic  deposition  is  so  common,  and 
because  the  final  settling  particles  of 
turbidity  currents  are  commonly  clay 
particles,  the  location  of  low  sound- 
speed  sediments  is  predictable.  Often, 
the  tops  of  cores  taken  in  the  deep 
ocean  sediments  are  soupy,  further  sup¬ 
porting  Urick1 s  theory.  The  finest  par¬ 
ticle-sized,  deep-water  surface  sedi¬ 
ments  the  oozes  and  clays — have  sound 
speeds  of  about  0.985  of  the  sound  speed 
of  the  ocean-bottom  water  at  the  water- 
sediment  interface . 

Sound-speed  measurements  have  been  made 
on  linered  cores  from  the  upper  10  m  of 
sediment  by  many  investigators.  A  common 
method  is  to  measure  the  time  delay  of  a 
400  kHz  signal  transmitted  across  the 
diameter  of  a  core  as  compared  to  that 
through  a  distilled  water  sample  in 
order  to  calculate  the  compressional 
wave  speed  of  the  sediment.  The  calcula¬ 
ted  speed  is  dependent  on  the  laboratory 
conditions  and  needs  to  be  corrected  to 
a  standard  temperature  for  comparison 
with  other  measurements  or  to  in  situ 
temperature  and  pressure  for  absolute 
speed  determinations.  In  this  report, 
all  sound-speeds  have  been  adjusted  to 
laboratory  conditions  of  one  atmosphere 
and  20°C  for  convenience  in  expressing 
the  actual  speeds  without  reference  to 
in  .situ  pressure  and  temperature.  Ac¬ 
cording  to  Wilson*  s  equation  [44],  sea¬ 
water  sound-speed  at  these  conditions  is 
1521.5-1522.5  m/sec  for  respective  sa¬ 
linities  of  34.68-35.46  parts  per  thou¬ 
sand.  These  salinity  values  are  well 
within  the  range  for  ocean-bottom  water 
in  the  northeast  Atlantic  [76].  There¬ 
fore,  the  adjusted  sound-speed  for  the 


ocean-bottom  water  will  be  considered  to 
be  1522  m/sec.  Clays  and  oozes  are  ap¬ 
proximately  1500  m/sec  under  the  same 
conditions.  Sediment  sound-speed  can 
also  be  expressed  as  the  ratio  of  the 
sound-speed  of  sediment  to  that  of  the 
bottom  water,  which  relieves  the  neces¬ 
sity  of  reference  to  a  standard  tempera¬ 
ture. 

The  extensive  collection  of  literature 
addressing  the  subject  of  relating  meas¬ 
urable  geologic  parameters  to  sound- 
speed  has  been  summarized  by  Hamilton 
[41].  Most  of  the  studies  are  based  on 
sediments  from  the  upper  few  meters  of 
the  ocean  floor;  therefore,  the  statis¬ 
tical  correlations  observed  apply  to 
only  surficial  sediments.  Hamilton  em¬ 
phasizes  the  need  to  use  as  much  infor¬ 
mation  as  possible  before  relying  on  a 
general  relationship.  One  commonly  used 
relationship  permits  an  estimate  or  pre¬ 
diction  of  sound— speed  if  the  porosity 
is  known.  This  is  because  the  sound- 
speed  is  largely  determined  by  the  com¬ 
pressibilities  of  the  pore-water  and  the 
solid  material  and  the  percentage  of 
these  two  constituents.  Hamilton  sepa¬ 
rates  sediments  from  abyssal  hills, 
abyssal  plains,  and  continental  terra¬ 
ces  and  slopes  in  deriving  the  sound- 
speed-porosity  relationships. 

Other  relationships  can  also  be  used  in 
the  prediction  but  with  somewhat  less 
reliability.  Anderson  [62]  derived  re¬ 
lationships  using  data  from  several 
oceans  and  showed  that  the  highest  de¬ 
gree  of  correlation  exists  with  sound- 
speed  and  porosity  as  compared  to  sound- 
speed  and  other  parameters  such  as  mean 
grain  size,  carbonate  content,  or  shear 
strength.  The  relationship  between  grain 
size  and  sound-speed  is  complicated  by 
factors  related  to  the  sediment  type  and 
its  effect  on  porosity.  Surficial  sedi¬ 
ments  of  fine  mean  grain  size  normally 
consist  of  well— sorted  clay  particles. 
These  sediments  have  high  porosities 
due,  in  part,  to  the  adsorbed  water  on 
the  large  surface  area  of  clay  minerals. 
On  the  other  hand,  sediments  of  larger 
mean  grain  size  (silt-sand)  are  not 
well-sorted  because  the  fine  particles 
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fill  the  voids  which  exist  in  the  pack¬ 
ing  of  larger  particles  of  uniform  size. 
This  results  in  a  reduction  in  porosity. 
It  therefore  follows  that  an  inverse 
relationship  between  grain  size  and 
porosity  exist. 

Calcareous  sediments,  however,  do  not 
follow  this  relationship.  Clay-sized 
carbonate  sediments  consist  of  cocco- 
liths,  which  are  solid  flat  plates. 
Larger-sized  calcite  particles  and  sil¬ 
icate  tests  of  radiolaria  consist  of 
hollow  porous  tests  which  permit  a  much 
higher  porosity  of  the  sediment  than 
solid  materials.  Schreiber  [77]  obser¬ 
ved  that  sediments  with  a  grain  size 
exceeding  0.008  mm  median  diameter*  and 
consisting  of  radiolarian  ooze  or  more 
than  50%  carbonate  content  have  much 
higher  porosities  than  other  sediments 
of  the  same  grain  size.  She  also  shows 
how  the  high  carbonate  and  radiolarian 
sediments  deviate  from  the  general  sound 
speed-median  diameter  relationship,  al¬ 
though  this  deviation  is  not  as  striking 
as  the  deviation  from  the  median  diam¬ 
eter-porosity  relationship. 

An  interesting  aspect  of  sound-speed  is 
its  relationship  to  density.  Theoreti¬ 
cally  sound-speed  decreases  with  in¬ 
creasing  density.  However,  empirically 
it  is  higher  in  sediments  of  greater 
density.  The  associated  decrease  of  the 
compressibility  with  an  increase  in  den¬ 
sity  is  the  overriding  factor  which 
results  in  an  actual  increase  in  sound- 
speed.  The  additional  factor  of  an  in¬ 
crease  of  rigidity  also  increases  the 
compressional  wave  velocity  and  adds  a 
shear  wave  component.  The  latter,  how¬ 
ever,  is  probably  negligible  near  the 
surface  and  is  often,  although  perhaps 
erroneously,  ignored  in  the  intermediate 
depth  sediments. 


*The  median  and  mean  grain  size  are 
equivalent  only  if  a  symmetrical  distri¬ 
bution  of  grain  sizes  based  on  weight 
exists.  Horn  et  al .  [78]  states  that  the 
mean  is  a  better  measure  of  the  distri¬ 
bution  of  sizes,  and  this  description  is 
now  generally  used  when  discussing  sound 
speed  of  sediments. 


Sound-Speed  Profiles 

A  positive  sound-speed  gradient  with 
depth  in  the  upper  sediments  was  first 
postulated  by  Hill  [4]  in  1952,  and 
based  on  in  situ  experimental  results. 
He  derived  a  gradient  of  2.5  sec“^  at 
a  site  in  the  Norwegian  Sea.  This  value 
is  high  compared  to  gradients  reported 
since  then,  but  is  comparable  to  our 
results.  A  gradient  is  also  inferred 
from  geological  factors  such  as  compac¬ 
tion,  pressure,  and  temperature  grad¬ 
ients  in  the  sediments.  Evidence  from 
later  work  suggests  that  the  gradient 
decreases  with  depth.  In  this  report, 
several  different  methods  are  used  to 
estimate  the  sound-speed  profile  in  the 
upper  200  m  of  sediment.  These  methods 
are  described  in  Appendix  E. 

Attenuation 

The  attenuation  constant  <  is  the  con¬ 
stant  of  proportionality  relating  atten¬ 
uation  a  to  frequency.  Attenuation  is 
defined  as  the  loss  of  energy  by  absorp¬ 
tion  (dissipation  due  to  heat)  or  other 
mechanisms,  such  as  scattering,  which 
are  not  attributed  to  geometric  spread¬ 
ing.  Hamilton  [45]  discusses  physical 
factors  which  contribute  to  attenuation, 
such  as  particle  angularity  and  composi¬ 
tion.  These  factors  are  closely  related 
to  particle  size  distribution.  Hamilton 
[45,  41]  has  summarized  attenuation 

values  in  sediments  determined  by  sev¬ 
eral  investigators.  From  this  compila¬ 
tion,  he  derived  an  empirical  curve 
relating  to  mean  grain  size.  This  rela¬ 
tionship  is  used  in  the  present  study  to 
estimate  the  values  of  k  based  on  the 
mean  grain  size  distribution  as  deter¬ 
mined  from  core  samples.  Where  core 
samples  are  not  available,  the  mean 
grain  size  is  estimated. 
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Appendix  E. 

Methods  of  Estimating  Sediment  Sound-Speed  Profiles 

Methods  of  estimating  sediment  sound- 
speed  profiles  are  presented  in  detail 
in  this  appendix.  The  profiles  were  used 
in  Appendix  C  for  the  acoustic  stations. 
Some  of  the  profiles  were  general  in 
that  they  are  derived  from  data  at  dis¬ 
tant  sites  where  a  similar  sediment  type 
exists.  On  the  other  hand,  other  pro¬ 
files  have  parameter  values  unique  to  a 
site  not  too  distant  from  the  acoustic 
stations.  These  profiles  have  been 
chosen  as  the  best  estimate  for  each 
station  and  are  indicated  by  an  arrow  in 
Figure  10.  The  best  estimate  profile  can 
be  visually  compared  with  the  general 
profiles  and  the  acoustically  derived 
profile. 

1 .  Laboratory  Methods 

In  the  reports  of  numerous  laboratory 
experiments  involving  sea  floor  sedi¬ 
ments,  sound-speed  measurements,  high 
pressure,  and  compaction,  only  two 
sound-speed  profiles  applicable  to  sedi¬ 
ments  have  been  described  (Laughton  [40] 
and  Cernock  [50]).  The  authors  of  both 
reports  have  conducted  laboratory  ex¬ 
periments  involving  sound-speed  vs. 
overburden  pressure  (as  opposed  to 
hydrostatic  pressure  alone)  on  uncon¬ 
solidated  sediment  samples. 

Laughton’s  [40]  measurements  included  a 
sample  with  high  carbonate  content  (54%) 
and  a  density  of  1.59  g/cm^  from  the 
North  Atlantic  taken  at  a  water  depth  of 
2577  m.  Because  of  the  high  carbonate 
content  of  the  sediments  at  Station  8, 
his  results  are  used  to  suggest  how  the 
profile  of  high-carbonate  sediments  may 
possibly  appear.  Hamilton  [79]  has  made 
the  conversions  from  pressure  to  depth 
in  sediment  for  Laughton’s  globigerina 
ooze.  Laughton  used  equipment  designed 
for  high  pressures  (up  to  1024  kgm/ 
cm^)  and  measured  only  four  values  of 


sound-speed  at  pressures  corresponding 
to  depths  of  burial  less  than  400  m. 

This  includes  the  surface  value  of  1.62 
m/sec  and  the  values  at  135,  260,  and 
360  m.  Actual  in  situ  values  are  expect¬ 
ed  to  be  higher  than  Laughton’s  compu¬ 
tations  because  he  had  made  no  correc¬ 
tion  for  a  temperature  gradient  in  the 
sediments.  However,  his  surface  value  is 
100  m/sec  greater  than  the  sediments  at 
Station  8,  thus,  the  actual  profile  is 
unrealistic  in  the  near-surface  area. 
Cernock  [50]  suggested  that  Laughton’s 
initial  values  are  too  high  because  of 
his  experimental  procedure.  Following 
this  suggestion,  Laughton’s  surface 
value  is  replaced  by  the  surface-sedi¬ 
ment  sound-speed  at  Station  8.  Laugh¬ 
ton’s  values,  with  this  modification, 
are  represented  by  the  curve  in  Figure 
65. 

Cernock  [50]  made  measurements  on  a  core 
sample  126  cm  deep  in  an  abyssal  plain 
(water  depth  3636  m)  in  the  Gulf  of 
Mexico.  The  sample  is  clay  which  has  a 
low  carbonate  content  (8%)  and  a  density 
of  1.31  g/cm^.  The  sound-speed  ratio 
of  his  sample  is  the  same  as  that  re¬ 
ported  for  Station  2  (0.974).  A  correc¬ 
tion  of  26  m/sec  is  added  to  each  of  his 
calculations  to  correct  for  the  differ¬ 
ence  in  bottom-water  sound-speed  between 
the  Gulf  of  Mexico  station  and  Station  2 
in  the  Iberian  Abyssal  Plain.  The  re¬ 
sulting  curve  is  shown  in  Figure  54. 

As  pointed  out  by  Cernock  [30],  results 
obtained  by  this  method  must  be  used 
with  caution  because  of  many  unknowns. 
Laboratory  experiments  are  performed 
over  a  short  time  period  without  the 
understanding  of  the  chemical  and  phys¬ 
ical  effects  which  occur  over  geologic 
time.  Also,  present  and  past  environ¬ 
mental  conditions  (pressure,  temper¬ 
ature,  and  salinity)  are  unknown.  In 
addition,  laboratory  measurements  are 
made  at  frequencies  on  the  order  of  100 
kHz,  whereas  in  situ  acoustic  measure¬ 
ments  were  taken  at  frequencies  on  the 
order  of  100  Hz. 
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2.  CRIC  Method 

The  CRIC  method  is  useful  in  estimating 
the  sound-speed  profile  in  a  section  of 
sediments  where  the  surface  sound-speed 
ratio,  CR,  and  the  average  or  interval 
sound  speed,  IC,  is  known.  In  partic¬ 
ular,  this  method  is  applicable  to  sur¬ 
face  sections  of  DSDP  sites  where  IC  can 
be  calculated  by  dividing  the  thickness 
of  a  sediment  section  by  the  travel  time 
through  the  sediment.  The  section  is 
bounded  by  acoustic  reflectors  which  can 
be  identified  by  the  drilling  process; 
e.g.,  the  water  interface  on  top  and  a 
sandstone  layer  or  basalt  at  the  bottom. 
The  thickness  is  measured  by  the  drill¬ 
ing  depths  at  these  interfaces,  and  the 
travel  time  is  obtained  from  seismic 
profiler  records.  The  sediment  surface 
sound-speed  cQ  is  often  indirectly 
determined  by  measuring  the  sound-speed 
of  a  core  sample  in  the  laboratory.  The 
ratio  of  the  sediment  sound-speed  to 
that  of  the  bottom  water  sound— speed  CR 
can  then  be  calculated  from  the  labor¬ 
atory  measurements.  The  bottom  water 
sound-speed  cB  can  be  established  from 
historical  or  in  situ  data.  Thus, 

cQ  =  CR  x  cB  (1) 

CR  may  be  used  from  sediments  at  other 
locations  if  they  are  similar  to  those 
at  the  site. 

The  value  of  sound-speed  as  a  function 
of  depth  z  can  be  determined  if  a  linear 
increase  in  sound— speed  with  depth  is 
assumed.  That  is,  a  constant  gradient  g 
is  assumed: 

C(z)  =  CG  +  gz  (2) 

The  two-way  travel  time  to  depth  z  as 
observed  on  the  profiler  record  is 

t  =  (2/g)  ln(l  +  gz/c0)  (3) 

By  reiteration,  the  value  of  g  which 
satisfies  Equation  3  can  be  calculated. 

This  method  was  used  for  determining  the 
values  of  g  at  DSDP  Sites  118  and  136. 
These  values  were  used  to  estimate  the 


sound-speed  profiles  at  Stations  2,  4, 

5,  and  8. 

The  method  is  illustrated  in  Figure  68 
where  h,  the  thickness  of  the  upper 
layer,  is  used  instead  of  z.  Other  lines 
of  evidence  discussed  in  following  sec¬ 
tions  indicate  that  the  true  sound-speed 
profile  is  not  a  constant  gradient,  but 
tends  to  have  a  decreasing  gradient  with 
depth.  Therefore,  the  value  of  the  con¬ 
stant  gradient  calculated  by  the  CRIC 
method  is  the  lower  limit  of  the  actual 
gradient  in  the  upper  part  of  the  sec¬ 
tion,  and  the  upper  limit  of  the  actual 
gradient  in  the  lower  part. 

DSDP  Site  1 18 

Site  118  [49]  is  located  in  the  Biscay 
Abyssal  Plain  (see  Fig.  52).  On  the  pro¬ 
filer  records  obtained  at  the  site,  a 
strong  reflecting  horizon  called  the 
mid-R  reflector  appears  at  0.48  seconds 
two-way  travel  time  below  the  ocean 
floor.  The  mid-R  reflector  was  found  to 
be  several  sandstone  layers  lithified 
from  coarse  turbidite  layers  in  the 
depth  interval  from  400  to  450  m. 

The  in  situ  bottom  temperature  and  sali¬ 
nity  at  Site  118  are  2.6°  and  34.9 
°/00,  respectively  [76].  The  corre¬ 
sponding  bottom  water  sound-speed  (cB) 
at  4901  m  is  1545  m/sec  [44].  The  value 
of  CR  is  not  available  in  the  site  re¬ 
port,  but  is  available  from  an  area  of 
similar  sediments.  Based  on  Ruddiman  and 
Glover's  study  [33],  the  Tagus  Abyssal 
Plain  is  such  an  area.  Using  the  ratio 
of  sediment  to  bottom-water  sound  speed 
of  0.995  observed  at  Stations  4  and  5  in 
the  Tagus  Abyssal  Plain,  the  calculated 
surface-sediment  sound-speed  at  Site  118 
is  1537  m/sec.  The  gradient  g  which 
satisfies  Equation  3  using  this  value  of 
cQ  and  the  observed  values  of  t  (0.48 
sec)  and  z  (400  m)  is  0.67  sec™^-.  The 
sound-speed  profile  defined  by  Equation 
2  using  this  value  of  g  is  used  as  the 
best  estimate  for  Stations  4  and  5. 

The  value  of  g  determined  by  the  CRIC 
method  at  Site  118  is  also  estimated  to 
be  the  value  at  Station  2  in  the  Iberian 
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Abyssal  Plain.  At  Station  2,  the  source 
of  sediments  is  more  distant  than  at 
Site  118  or  at  Stations  4  and  5.  This 
results  in  a  smaller  grain  size  in  the 
basal  turbidites  in  the  northern  Iberian 
Abyssal  Plain,  and  is  corroborated  by 
the  lower  ratio  of  sediment  to  bottom- 
water  sound-speed  ratio  (0.977  as 
compared  to  0.995  in  the  Tagus  Abyssal 
Plain).  The  smaller  grain  size  suggests 
a  respectively  slower  accumulation  of 
sediments  deposited  by  turbidites.  Thus, 
sediments  of  a  given  age  are  less  deeply 
buried  at  Station  2  than  at  Site  118. 
Hamilton  et  al.  [42]  measured  gradients 
at  different  locations  on  the  Bengal 
Fan,  with  different  depths  of  burial  to 
sediments  of  the  same  age.  He  observed 
that  the  gradients  in  the  Bengal  Fan  are 
greater  farther  away  from  the  source  of 
sediments  where  the  accumulation  rates 
are  less.  He  suggests  that  lithif ication 
increases  with  age  with  an  accompanying 
increase  in  sound-speed.  Therefore, 
sediments  of  a  given  age  being  less 
deeply  buried  at  the  distant  stations 
account  for  the  larger  gradients.  A 
similar  argument  might  be  used  here  to 
suggest  that  the  gradient  at  Station  2 
is  higher  than  that  at  Site  118. 

However,  the  information  is  insuf¬ 
ficient  to  support  an  adjustment  to  the 
value  of  0.67  sec”-'-. 

DSDP  Site  136 

The  CRIC  method  can  also  be  applied  to 
DSDP  Site  136.  The  sediments  at  this 

site  are  discussed  in  Appendix  C  for 
Station  8.  The  value  of  cQ  =  1530 

m/sec  (CR  =  1.00)  is  assumed  at  Site 

136.  The  basalt  interface  was  drilled  at 
308  m.  The  corresponding  two-way  travel 
time  observed  on  the  profiler  records 
from  the  DSDP  report  [59]  is  0.33  sec. 
The  value  of  g  satisfying  Equation  3  is 
2.33  sec"  . 

The  value  of  2.33  sec”1  is  consider¬ 
ably  higher  than  those  values  esti¬ 
mated  at  the  other  stations .  This  value 
is  a  function  of  the  assumed  surface 
sound-speed  CQ.  However,  the  value  of 
CQ  is  unlikely  to  be  more  than  30 


m/sec  in  error.  (This  would  allow  for  a 
low-speed  bottom,  CR=0.98,  or  a  high¬ 
speed  bottom  as  measured  at  Station  8, 
CR=1.013.)  This  amount  of  error  does  not 
significantly  affect  the  determination 
of  g:  for  CQ  varying  +  30  m/sec,  g 

varies  +  0.23  sec"-'-. 

Two  possible  sources  of  error  are  appar¬ 
ent  in  the  DSDP  data;  neither  is  taken 
into  account  in  this  analysis:  (a)  The 
measured  thickness  of  sediment,  deter¬ 
mined  by  the  length  of  the  drill  cable, 
could  be  in  error  due  to  the  uncertainty 
as  to  when  the  drill  bit  actually  enter¬ 
ed  the  ocean  floor.  This  error  probably 
is  no  greater  than  5  m.  (b)  The  seismic 
profile  may  not  have  been  taken  right  at 
the  core  site — the  topography  of  the 
basement  is  very  rough  and  this  could  be 
in  error.  Or,  the  time  read  on  the  seis¬ 
mic  profile  may  be  in  error.  In  consid¬ 
ering  this  error,  the  site  report  states 
0.33  sec  to  basement  [59].  The  pick  may 
be  a  subjective  or  average  estimate.  In 
measuring  the  time  scale  of  the  record 
in  Figure  3  of  reference  [59],  it  ap¬ 
pears  to  be  0.37  sec.  It  is  assumed, 
however,  that  the  authors  of  the  report 
had  sufficient  information  to  adequately 
define  the  travel  time  of  0.33  sec. 

3.  Seismic  Refraction  Method 

Houtz  and  Ewing  [39]  developed  a  sound- 
speed-depth  curve  for  unconsolidated 
sediments  based  on  wide-angle  reflec¬ 
tion  data.  The  data  were  obtained  from 
the  results  of  60  seismic  refraction 
profiles  taken  predominantly  along  the 
western  North  Atlantic  continental  rise 
and  the  Bermuda  Rise.  The  geometry  used 
in  this  technique  is  equivalent  to  that 
used  in  the  present-day  sonobuoy  method. 
In  five  profiles,  observed  arrivals  were 
interpreted  as  paths  refracted  in  the 
sediments  due  to  sound-speed  gradients 
(i.e.,  turnaround  rays).  The  average  of 
the  computed  gradients  is  1.85  sec”^-. 
In  the  same  profiles,  arrivals  were  ob¬ 
served  which  reflected  from  the  base  of 
the  first  sedimentary  layer  (validated 
by  the  T^  -  analysis).  The  cal¬ 

culated  interval  velocity  implied  a 
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decrease  in  the  gradient  from  that  cal¬ 
culated  in  the  near-surface  refracted 
arrival . 

Houtz  and  Ewing  [39]  compiled  33  deep¬ 
water  sedimentary  sound-speeds  calcu¬ 
lated  from  the  subbottom  reflection 
curves  on  the  refraction  profiles •  These 
values  were  plotted  against  the  travel 
time  in  the  sediment  to  the  midpoint  of 
the  layer,  where  the  midpoint  time  is 
one-half  the  interval  time  of  the  layer. 
By  choosing  a  sound-speed  function  which 
provides  a  decreasing  gradient  with 
depth,  sound-speed  was  derived  as  a 
function  of  time  t  in  seconds: 

C  =  1520  (1  +  8.75)1/5  m/sec  (4) 

or  expressed  as  a  function  of  depth  z  in 
km: 

C  =  1520  (1  +  6.91Z)1/6  m/sec  (5) 

Equation  5  is  plotted  in  Figure  65. 
Equation  5  less  5  m/sec  for  a  surface 
sound-speed  correction  is  plotted  in 
Figure  54. 

4.  Wide-Angle  Reflection  Method 

Extensive  wide-angle  reflection  pro¬ 
files  using  the  so nobuoy  technique  [80] 
provide  average  sound  speeds  for  sed¬ 
imentary  layers  in  many  ocean  areas.  We 
have  used  the  results  of  this  method  in 
predicting  possible  sound-speed  pro¬ 
files  at  several  of  the  stations. 

Data  for  Turbidite  Sediments 

Hamilton  et  al.  [51]  have  compiled  data 
from  17  ocean  areas  with  sediments 
formed  mainly  by  turbidite  deposition. 
The  average  gradient  g  is  given  by 
equation  6  where  CQ  is  the  sediment 
surface  sound  speed,  and  IC  is  the 
interval  sound  speed  over  the  layer  of 
thickness  h: 

g  =  2  (IC  -  CG)/h  (6) 

Equation  6  is  an  approximation  because 
the  value  of  IC  is  assumed  to  be  iden¬ 
tical  to  the  value  of  sound  speed  at  the 


midpoint  of  the  layer.  Hamilton  et  al. 
used  the  directly  measurable  parameter 
of  travel  time  rather  than  the  calcul¬ 
ated  thickness  to  compile  average  grad¬ 
ients.  Their  objective  was  to  produce  a 
simple  means  to  estimate  thickness  from 
the  travel  time  observed  on  a  profiler 
record.  They  plotted  the  average  gra¬ 
dient  versus  one-way  travel  time  at  in¬ 
crements  of  0.1  sec.  Results  indicate  a 
constant  decrease  of  gradient  with  one¬ 
way  travel  time,  t,  as  expressed  by 

g  =  1.316  -  1.117t  (7) 

The  sound  speed  is  given  by 

C  =  CQ  +  gh  (8) 

where 

h  =  CQ  (eSt  -  l)/g  (9) 

and  using  the  value  of  g  corresponding 
to  t  in  Equation  7 •  By  assuming  a  sed¬ 
iment  surface  sound  speed  of  CQ 
1515.0  m/sec,  the  instantaneous  sound- 
speed  c  has  been  plotted  to  produce  the 
profile  presented  in  Figure  54  for  the 
Iberian  Abyssal  Plain.  The  profile  is 
also  shown  in  Figure  10  for  Stations  4 
and  5  in  the  Tagus  Abyssal  Plain. 

Data  for  Highly  Calcareous  Sediments 

Houtz,  Ewing  and  Buhl  [60]  have  com¬ 
piled  sonobuoy  data  collected  in  dif¬ 
ferent  areas  of  the  Pacific  Ocean.  They 
plotted  the  interval  sound  speed  of  a 
layer  against  one-way  travel  time  in 
seconds  to  the  midpoint  of  the  layer  for 
all  the  valid  data  in  each  area.  They 
derived  an  equation  for  each  area  by 
making  a  least  squares  fit  of  a  poly¬ 
nomial  for  each  set  of  data.  One  set  of 
data  is  from  pelagic  sediments  in  the 
equatorial  zone.  Houtz  et  al .  [60]  noted 
the  high  sound-speed  gradient  for  this 
set  of  data  representing  sediments  of 
high  carbonate  content.  Hamilton  et  al . 
[42]  have  recommended  using  the  equation 
for  this  set  of  data  for  pelagic  sedi¬ 
ments  anywhere  in  the  world  as  a  best 
estimate  of  sound-speed  if  there  is  no 
in  situ  data  available.  The  sediments 
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at  Station  8  are  pelagic,  with  a  high 
carbonate  content.  In  following  Hamil¬ 
ton's  recommendation,  this  equation  has 
been  used  in  deriving  a  profile  for  Sta¬ 
tion  8. 

The  Houtz,  Ewing,  and  Buhl  [60]  equa¬ 
tion  for  sound-speed  in  the  Pacific  equ¬ 
atorial  zone  is: 

C  =  1521  +  277 6t  +  2931t2  m/sec 

(10) 

where  t  is  the  one-way  travel  time  in 
seconds.  (Their  expression  of  sound 
speed  versus  depth  [in  their  Table  2]  is 
the  best  fit  for  data,  but  gives  an  in¬ 
terface  sound-speed  of  1600  m/sec.) 

In  order  to  use  this  data  for  the  near¬ 
surface  sound-speed  profile,  we  have 
used  the  value  of  sound-speed  from  Equa¬ 
tion  10  and  the  corresponding  depth  z 
from  Equation  11  for  the  same  value  of 
one-way  travel  time.  Equation  11  is  an 
exact  expression  for  depth  z  in  terms  of 
t  derived  by  integration  of  Equation  10: 

z  =  /  C  dt  =  1521t  +  1388t2  +  977t3m 

(U) 

Data  for  Station  9  Area 

Many  sonobuoy  stations  were  taken  on 
Lamont -Doherty  cruises  in  the  area  of 
Station  9,  as  shown  in  Figure  66.  The 
large  variety  of  layer  thicknesses  with 
respective  interval  sound-speeds  pro¬ 
vides  a  data  base  to  estimate  an  average 
sound-speed  profile  for  the  region.  This 
type  of  information  has  been  derived 
from  sonobuoy  data  by  Hamilton  et  al. 
[42]  for  several  areas. 

Twenty-three  of  the  stations  have  been 
separated  geographically  into  two  groups 
as  indicated  in  Figure  66:  Group  W 
(west)  as  denoted  by  triangles,  and 
Group  E  (east)  with  circles.  In  addi¬ 
tion,  the  four  stations  nearest  Station 
9  are  grouped  separately  and  are  indi¬ 
cated  by  squares;  all  others  are  indi¬ 
cated  by  crosses.  The  layer  thickness 


and  interval  speeds  as  well  as  deeper 
refraction  speeds  are  shown  for  all  the 
stations  in  Figure  69. 

Groups  E  and  W  data  are  considered  as 
two  separate  data  bases  to  estimate  a 
sound  speed  profile  near  Station  9.  For 
all  the  layers  in  a  group,  the  layer 
speed  is  plotted  against  the  depth  below 
seafloor  to  the  midpoint  of  the  layer. 
These  plots  are  shown  for  Group  W  and 
Group  E  in  Figures  70  and  71.  A  least- 
squares  second-degree  curve  for  speed  as 
a  function  of,  depth  of  burial  was  com¬ 
puted  and  drawn  on  the  figures . 

In  addition  to  the  data  base  points, 
points  for  the  layers  determined  from 
the  sonobuoy  stations  near  Station  9  are 
shown  on  both  figures  as  squares.  The 
curve  for  Group  E  appears  to  better  rep¬ 
resent  the  data  points  near  Station  9 
than  does  the  curve  for  Group  W.  There¬ 
fore,  the  equation  for  this  curve  is 
used  to  define  the  function  of  speed 
with  depth  for  this  area: 

C(z)  =  1505  +  1448z  -  305z^  m/  sec 

(12) 

where  z  is  in  kilometers  below  sea 
floor,  and  C  is  sound-speed.  The  aver¬ 
age  gradient  for  the  upper  200  m  is 
g  =  AC/Az  =  1.39  sec“-*-. 

The  in  situ  bottom-water  sound-speed 
range  over  the  range  of  water  depth  of 
all  the  sonobuoys  is  about  1515-1535 
m/sec.  This  variation  (due  to  differ¬ 
ences  in  hydrostatic  pressure)  will  af¬ 
fect  the  sound-speed  in  the  unconsoli¬ 
dated  sediments  by  the  same  amount.  How¬ 
ever,  the  variation  is  rather  insignifi¬ 
cant  on  the  scale  used  in  Figures  70  and 
71.  The  value  of  the  intercept  for  Group 
W  is  1435  m/sec,  which  is  obviously  too 
low,  whereas  1505  m/sec  for  Group  E  is 
within  the  range  expected  for  clays. 

The  least-squares  equation  was  recalcu¬ 
lated  using  the  data  weighted  according 
to  the  inverse  of  the  standard  error  of 
each  data  point,  and  adding  a  heavily 
weighted  (x20)  surface  sound-speed  of 
1520  m/sec.  The  effect  of  the  weighted 
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data  for  either  group  was  to  decrease 
the  near-surface  gradient  by  less  than 
0.2  sec“l.  (The  difference  due  to 
adding  the  heavily  weighted  surface 
sound  speed  is  less  than  0.1  sec"^ .) 


A  danger  exists  in  extrapolating  to  the 
surface  the  sound  speed  versus  depth 
curve  determined  predominantly  by  deeper 
layers.  Laboratory  experiments,  such  as 
by  Laughton  [40]  (Fig.  65),  show  that 
the  greatest  increase  of  speed  with 
overburden  pressure  due  to  depth  of  bur¬ 
ial  is  in  the  near-surface  sediments. 
Houtz  and  Ewing  fs  [39]  in  situ  experi¬ 
ments  support  this  concept.  Their  data 
consist  only  of  midlayer  points  to  a 
maximum  depth  of  0.8  km.  Below  this 
depth  a  discontinuous  increase  in  speed 
exists  at  the  interface  of  unconsoli¬ 
dated  sediments.  As  the  method  used  to 
derive  Equation  12  includes  high-speed 
deeper  layers  in  a  model  with  a  contin¬ 
uous  sound-speed  function,  the  gradient 
near  the  surface  will  be  in  error.  If 
one  could  separate  sound-speed  deter¬ 
minations  of  the  unconsolidated  sedi¬ 
ments  in  the  sonobuoy  data,  the  estimate 
would  be  more  accurate. 


unconsolidated  sediments.  The  calculated 
average  gradients  using  layers  with 
midpoint  depths  less  than  0.5  and  1.0 
km  are  1.1  and  1.2  sec"-*-,  respec¬ 

tively.  In  the  second  approach,  layers 
with  speeds  less  than  a  specified  sound 
speed  are  assumed  to  consist  of  only 
unconsolidated  sediments.  The  calcu¬ 
lated  average  gradients  using  layers 
with  speeds  less  than  2000  and  2300 
m/sec  are  0.8  and  1.2  sec"^ ,  respect¬ 
ively.  The  standard  deviations  of  each 
of  these  four  determinations  are  high, 
ranging  from  0.6  to  0.9  sec-1,  which 
shows  the  variability  of  the  sediment 
sound-speed  profiles  in  this  area. 
Throughout  the  DSDP  reports,  the  report¬ 
ed  sound-speeds  for  the  unconsolidated 
sediments  overlap.  A  similar  difficulty 
arises  in  attempting  to  determine  a  cri¬ 
tical  depth  which  separates  the  uncon¬ 
solidated  and  consolidated  sediments. 
Nevertheless,  the  average  of  the  four 
calculated  gradients  is  about  1.1 
sec*-1,  which  is  less  than  one-half  a 
standard  deviation  from  each  gradient 
value . 

5.  Surf icia I  Sediment  Data  Method 


An  alternate  method  of  calculating  the 
gradient  is  used  when  considering  only 
near-surface  layers.  Rather  than  esti¬ 
mating  a  sound-speed  profile  represen¬ 
ted  by  a  least-squares  curve,  only  two 
points  are  determined.  The  upper  point 
is  the  average  water-sediment  sound- 
speed  value  (1525  m/sec).  The  lower 
point  has  a  sound-speed  which  is  the 
average  speed  of  all  the  layers  consid¬ 
ered.  The  depth  at  this  point  is  the 
average  of  the  midpoints  of  these 
layers.  The  gradient  is  then  the  dif¬ 
ference  in  the  two  sound-speed  values 
divided  by  the  depth.  This  value  of  the 
gradient  is  about  0.9  times  the  value 
which  satisfies  the  constant-gradient 
model  with  the  same  speed  at  the  mid¬ 
point  depth. 

Two  approaches  to  determine  a  gradient 
in  the  above  manner  are  used  here.  In 
the  first  approach,  layers  with  a  mid¬ 
point  depth  less  than  a  specific  depth 
are  assumed  to  consist  only  of 


A  sequence  of  empirical  and  theoretical 
relationships  are  used  to  develop  a 
sound-speed  profile  at  the  DSDP  Site 
136.  Commencing  with  a  semi-theoretical 
porosity-depth  profile  based  on  surfi- 
cial  sediment  type,  and  applying  an  em¬ 
pirical  porosity-sound-speed  relation¬ 
ship  also  based  on  surficial  sediment 
type,  a  sound-speed  profile  is  derived 
which  is  then  corrected  for  in  situ 
pressure  and  temperature.  The  pressure 
and  temperature  corrections  applied  to 
sediments  are  those  established  for 
water.  However,  in  order  to  estimate  the 
temperature,  an  empirical  relationship 
between  temperature  and  porosity  is 
used.  Within  this  sequence  of  relation¬ 
ships  are  many  implied  assumptions  which 
must  be  kept  in  mind  when  assessing  the 
final  results. 

The  DSDP  report  [59]  provides  some  poro¬ 
sity  data  for  Site  136.  Porosity  was 
measured  by  two  methods  on  board  the 
GLOMAR  CHALLENGER  on  all  the  cores  from 
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the  drill  hole:  GRAPE  (Gamma  Ray  Attenu¬ 
ation  Porosity  Evaluator)  and  a  volu¬ 
metric  determination.  The  GRAPE  method 
uses  an  average  of  800  determinations 
per  section  to  compute  a  single  value 
for  porosity,  whereas  the  volumetric 
method  determines  one  value  for  one 
sample  of  a  section.  The  range  of  poro¬ 
sity  values  of  all  sections  for  each 
core  by  the  these  two  methods  are  indi¬ 
cated  in  Figure  72.  Because  these  meas¬ 
urements  are  made  after  the  sample  has 
suffered  a  release  of  overburden  pres¬ 
sure  from  the  in  situ  environment,  a 
correction  should  be  made  to  in  situ 
conditions.  This  correction  is  shown  by 
the  length  of  the  arrows  from  the  mean 
of  the  GRAPE  porosity  range  for  each 
core.  The  correction  adjusts  the  poro¬ 
sity  to  the  in  situ  porosity  based  on 
depth  of  burial  using  the  values  deter¬ 
mined  by  Hamilton  [81]. 

Three  theoretical  curves  by  Hamilton 
[81]  are  presented  in  Figure  72  that 
show  the  relationship  between  porosity 
and  depth  of  burial  for  carbonates,  tur- 
bidites,  or  pelagic  clays.  The  actual 
porosity  measurements  of  the  Site  136 
cores  vary  considerably;  therefore, 
these  measurements  cannot  be  used  to 
justify  the  application  of  a  particular 
curve.  Because  of  the  predominantly  cal¬ 
careous  sediment  at  Site  136,  the  "car¬ 
bonates'*  curve  is  used  in  determining  a 
smooth  porosity  profile  at  the  site.  A 
sound-speed  porosity  relationship  is 
used  to  transform  this  curve  into  a 
sound-speed-versus-depth  curve. 

Hamilton  [52]  derived  curves  represent¬ 
ing  sound-speed  as  a  function  of  poros¬ 
ity  for  different  types  of  sediments. 
His  curve '"T"  is  based  on  data  derived 
from  high  carbonate  content  sediments 
found  on  the  continental  terrace.  This 
set  of  data  encompasses  porosities  of 
35-80%.  The  sediments  at  Site  136  are 
probably  similar  in  porosity  and  car¬ 
bonate  content  to  the  sediments  in 
Hamilton’s  data  base.  Hayes  et  al.  [59] 
calculated  a  very  high  average  sound 
speed  of  1867  m/sec  from  the  two-way 
travel  time  through  the  308  m  sediment 
column  at  Site  136.  The  high  values  are 


consistent  with  Hamilton’s  "T"  curve. 
This  curve  is  among  the  sound-speed 
porosity  curves  summarized  by  Anderson 
[62]  and  has  the  highest  sound-speed 
values  for  the  porosity  range  of  50-70%, 
i.e.,  the  range  expected  at  Site  136. 
Using  the  "T"  curve,  the  porosity  at 
each  depth  is  used  to  estimate  the 
sound^speed  at  that  depth.  However,  the 
application  of  porosity  relationships 
from  surficial  sediments  to  buried  sedi¬ 
ments  presents  a  problem  which  is  dis¬ 
cussed  later. 

In  determining  the  in  situ  sound-speeds 
for  buried  sediments,  increases  in  temp¬ 
erature  and  hydrostatic  pressure  below 
the  ocean  floor  must  be  considered.  A 
common  practice  is  to  assume  that  the 
sound-speed  of  sediment  varies  with 
pressure  and  temperature  to  the  same 
degree  as  does  sea  water,  as  defined  by 
Wilson's  equation  given  in  Reference 
[44].  Considering  only  the  hydrostatic 
pressure  increase  from  the  ocean  floor 
(4169  m)  to  308  m  into  the  sediment 
(4477  m),  the  speed  increases  5.5  m/sec. 

To  determine  the  in  situ  temperature 
below  the  sea  floor,  heat  flow  infor¬ 
mation  must  be  considered.  Bullard  [82] 
suggests  a  general  relationship  between 
the  thermal  conductivity  k  of  sediment 
and  water  content  w: 

k  -  (168  +  678  w)"1  cal/cm/sec/°C 

(13) 

where  w  can  be  expressed  in  terms  of  the 
porosity  n,  water  density  p  w,  and 
sediment  grain  density  Ps: 

w  =  [ps/pw  (1/n-l)  +1]-1 

(14) 

The  temperature  at  the  lower  boundary  of 
the  depth  interval  zn  is 

Tn  =  -QAzn/k  +  Tn-1 

(15) 
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where  Tn_-^  is  the  temperature  at  the 
upper  boundary  of  zn,  and  Q  is  the 
heat  flow.  The  average  heat  flow  in 
this  area  is  -1.06xl0-6  cal/cm2/sec 
[83],  representing  upward  flow  in  the 
negative  z  direction,  and  the  bottom 
temperature  is  2.3°C  [76].  Using 

^  s  =  2*7  and  P  w  =  1.05  g/cm-2  in 

Equation  14  and  the  porosities  obtained 
from  the  curve  in  Figure  72,  the  calcu¬ 
lated  in  situ  temperature  increases  with 
depth  from  2.3°C  at  the  ocean  floor  to 
15.9°C  at  308  m  into  the  bottom.  The 
resulting  increase  in  sound-speed  is 
50.7  m/sec  over  the  308  m  interval. 

The  above  pressure  and  temperature 
effects  have  been  applied  in  deriving 
the  in  situ  profile.  A  problem  exists, 
however,  in  using  sound— speed— porosity 
relationships  determined  from  sediments 
on  the  sea  floor.  Hamilton's  [52]  sam¬ 
ples  are  all  from  the  upper  30  cm  of  the 
sea  floor.  Indeed,  most  valid  sound- 
speed  measurements  on  ocean  sediments 
are  taken  on  surficial  sediments.  (In 
deep  cores ,  such  as  the  DSDP  cores ,  the 
effect  of  the  release  of  the  overburden 
pressure  presents  a  disturbance  which 
has  unknown  effects  in  the  sound-speed 
determinations.)  Ocean— bottom  clay  sam¬ 
ples  which  have  not  been  previously  bur¬ 
ied  have  high  porosities.  When  buried, 
the  resulting  compaction  reduces  the 
porosity.  An  uncompacted  sample  at  the 
water-sediment  interface  with  the  same 
porosity  will  have  larger  grain  size, 
and  therefore  will  likely  consist  of 
carbonaceous  or  siliceous  tests  or 
quartz  grain  instead  of  clay  minerals. 
The  structural  matrix  is  different  as 
the  mineralogy  is  different;  therefore, 
the  sound-speed  of  the  surface  sample 
would  not  necessarily  be  the  same  as 
that  of  the  compacted  clay  of  the  same 
porosity.  Thus,  a  problem  exists  in  ap¬ 
plying  the  sound  speed— porosity  rela¬ 
tionships  derived  from  surface  sediment 
samples  to  buried  sediments. 
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Appendix  F. 

Water-Column  Sound-Speed  Structure  at  Each  Station 


The  sound-speed  profile  of  the  water 
column  at  each  acoustic  station  was 
determined  from  data  obtained  with  an 
SVSTD.  The  measurements  were  taken  at 
each  site  immediately  before  or  after 
the  reflectivity  experiment.  The  pro¬ 
files  for  each  station  are  shown  in 
Figure  73.  The  variation  of  water  depth 
along  the  bomb  run  is  also  noted. 

The  depths  of  the  source  and  receiver 
for  each  station  are  shown  on  the  res¬ 
pective  profiles  in  Figure  73.  The 
authors  chose  these  depths  prior  to  the 
cruise  to  obtain  the  optimum  geometry 
(discussed  in  Sect.  Ill)  without  con¬ 
sideration  of  the  sound-speed  profile. 
Near  midwater  receiver  depths  were  used 
at  all  stations,  except  Station  7,  where 
the  receiver  depth  was  the  same  as  the 
shot  depth. 

Fenner  and  Bucca  [84,  85]  have  describ¬ 
ed  the  water  mass  circulation  and  sound- 
speed  profiles  in  the  North  Atlantic. 
The  predominant  structure  of  each  pro¬ 
file  is  related  to  three  oceanographic 
features:  the  surface  layer,  the  inter¬ 
mediate  water,  and  the  deep  water.  The 
warm,  saline  water  of  the  Gulf  Stream/ 
North  Atlantic  Current  system  produces 
the  high  sound-speed  in  the  surface 
layer  above  a  strong  negative  gradient 
common  to  all  the  stations.  Below  2000 
m,  the  North  Atlantic  Deep  Water  pro¬ 
duces  a  characteristic  positive  gradient 
which  varies  only  slightly  between  deep 
ocean  basins.  Interjected  between  these 
two  water  masses  is  the  Mediterranean 
Intermediate  Water  (MIW)  which  flows  out 
of  the  Straits  of  Gibraltar.  This  body 
of  warm,  highly  saline  water  creates  a 
sound-speed  maximum  which  varies  between 
900  and  1200  m  at  our  stations,  and  also 
causes  the  pronounced  bichannel  struc¬ 
ture  in  the  profiles.  Consequently  the 
deep  sound-channel  axis  is  depressed 


from  about  1200  m  to  2000  m  in  the  pres¬ 
ence  of  MIW. 

The  MIW  flows  preferentially  along  sev¬ 
eral  routes  as  shown  in  the  location 
chart  in  Figure  73.  The  extremely  high 
sound  speed  at  the  maximum  for  Station  5 
(about  1518  m/sec)  indicates  that  this 
station  lies  on  a  preferential  MIW  flow- 
path.  The  sound-speed  maximum  for  Sta¬ 
tion  4  is  about  4  m/sec  less  than  that 
for  Station  5,  reflecting  a  greater 
distance  off  the  MIW  flow-path.  Sound- 
speed  maxima  for  Stations  7  and  8  (about 
1510  m/sec)  indicate  even  further  dilu¬ 
tion  of  the  MIW  to  the  southwest.  At 
Station  2,  the  sound-speed  maximum  of 
about  1510  m/sec  indicates  dilution  of 
MIW  to  the  northwest.  Here,  the  depth 
and  sound-speed  at  the  upper  channel 
axis  are  the  minima  encountered  in  our 
data  due  to  the  lower  temperature  at  the 
northern  location.  At  Station  9,  the  bi¬ 
channel  sound-speed  structure  has  been 
compressed  to  an  irregular  structure 
around  the  depth  of  the  deep  sound-chan¬ 
nel  axis,  as  expected,  in  the  region 
southwest  of  Madiera  due  to  further  di¬ 
lution  and  sinking  of  the  MIW.  The  45% 
MIW  concentration  shown  in  Figure  73 
southwest  of  Station  9  represents  an 
anomalous  condition  (D.  F.  Fenner,  pers. 
c  omm  * ) • 
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